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Abstract 


Connected-pipe,  subsonic  combustion  ramjet  and  ducted  rocket  performance  determination  procedures  used  by  the  NATO  coun¬ 
tries  have  been  reviewed  and  evaluated, 

A  working  document  has  been  produced  which  provides  recommended  methods  for  reporting  test  results  and  delineates  the  para¬ 
meters  that  are  required  to  be  measured, 

Explanations  and  detailed  numerical  examples  are  presented  covering  the  determination  of  both  theoretical  and  experimental  per¬ 
formances,  the  use  of  air  heaters  and  uncertainty  and  error  analyses. 


Abrege 


Les  m^thodes  de  determination  des  performances  des  statordacteurs  et  des  statofus6es  au  banc  d’essais  en  conduite  forc6e, 
utilis6es  au  sein  de  la  communaut6  de  I'OTAN,  ont  6te  examinees  et  evaluees, 

Un  document  de  travail  a  et6  6tabor6  afin  de  foumirdes  recommandations  concemant  la  presentation  des  r6sultats  d’essais  et  de 
pr6ciser  les  parametres  indispensable  s  S  mesurer. 

Des  explications  sont  donnfees  et  des  exemples  numeriques  detailles  sont  prfesentes  afin  de  determiner  les  performances 
theoriques  et  expert mentales,  incluant  I’utilisation  de  foyers  de  pr6chauffage  de  Tan  ainsi  que  I’emploi  de  procedures  d’analyse 
des  erreurs  et  des  incertitudes. 
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constant  or  speed  of  sound,  m/s 
critical  speed  of  sound  (M  =  1),  m/s 
geometric  area, 
velocity,  m/s 
discharge  coefficient 
thrust  coefficient 

specific  heat  at  constant  pressure,  J/kg/K 
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mass,  kg 

mass  flow  rate,  kg/s 
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universal  gas  constant,  kj /kmole/ K 
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propellant  burning  surface  area,  m^ 
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efficiency  parameter  based  on  the  ratio  of  experimental-to-theoretical  values  for  a  specific  performance 
parameter 

efficiency  parameter  based  on  the  ratio  of  experimental-to-theoretical  characteristic  velocity  (c*)  where 

the  experimental  c*  can  be  determined  from  either  measured  thrust  or  pressure 

efficiency  parameter  based  on  the  ratio  of  experimental-to-theoretical  vacuum  specific  impulse  (fj^c) 
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Carbon 

effective 

experimental 

final,  fuel  or  propellant  grain,  frozen 

ducted  rocket  gas  generator 

geometric 

Hydrogen 

water 

initial 

ideal  air 

injector 

load  cell 

measured 

Oxygen 

make-up  oxygen  (i.e.,  vitiator  oxidizer  plus  replenishment  oxygen) 
process 

local  or  shifting  equilibrium 
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total 

tare 

theoretical 
temperature  rise 
vacuum 
vitiator  fuel 
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vehicle  station  identifications  (see  Fig.  3.1  -  3.3) 

process  representation  between  two  stations  (e.g.  station  4  to  station  5) 
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Superscripts 

*  secondary  vehicle  inlet  designation  (see  Fig.  3.2)  or  Af  =  1  location 

molar  basis 


Acronyms 

AAAM  Advanced  Air-to-Air  Missile  (U.S.) 

ANS  Anti-Navire  Supersonique,.  anti-ship  supersonic  missile  (France) 

ASMP  Air-Sol  Moyenne  Portee,  air-to-surface  medium  range  missile  (France) 
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CARS  Coherent  anti  Stokes  Raman  spectroscopy 

CFD  Computational  fluid  dynamics 

DR  Ducted  rocket 

EFA  Experimental  Feststau  Antrieb,  ducted  rocket  (Germany) 

IR  Infrared 

LFRJ  Liquid  fuel  ramjet 

LIF  Laser  induced  fluorescence 

MMH  Monomethylhydrazine 

MPSR  Missile  Probatoire  Statofusee  “ Rustique" ,  “Riistijae”. ducted  rocket  (France) 
NASP  National  Aero  Space'  Plane  (U.S.) 

NG  Natural  Gas 

SA4  Surface-to-air  missile  type  no.  4  (USSR) 

SA6  Surface-to-air  missile  type  no.  6  (USSR) 

SFRJ  Solid  fuel  ramjet 

SI  Systeme  international  d’unites  (International  system  of  units) 

SLAT  Supersonic  Low  Altitude  Target  missile  (U.S.) 

UDMH  Unsymmetrical  Dimethylhydrazine 

VFDR  Variable  Flow  Ducted  Rocket  (U.S.) 


Definitions 

Heat  of  formation  (AhJ)  —  Increment  in  enthalpy  associated  with  the  reaction  of  forming  the  given  compound 
from  its  elements,  with  each  substance  in  its  thermodynamic  standard  state  at  298. 15K.  (Also  referred  to  as  standard 
enthalpy  of  formation). 

Make-up  Oxygen  —  The  sum  of  oxygen  flow  rates,  vitiator  oxidizer  and  replenishment  oxygen,  that  must  be  added 
to  maintain  the  “mole  or  mass  fraction  of  oxygen  in  air”  in  the  vitiated  air  stream  supplied  to  the  propulsion  system. 
Stream  Thrust  —  Fs  =  me  +  pA 
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1  Introduction 


Many  NATO  nations  are  now  conducting  research  and  de¬ 
velopment  of  ramjets  for  supersonic,  extended  range  mis¬ 
siles  and  projectiles  In  the  context  of  this  report  the 
ramjet  is  taken  to  be  a  generic  class  of  propulsion  devices 
which  comprise  the  liquid  fuel  ramjet,  the  solid  fuel  ramjet 
and  the  ducted  rocket  (sometimes  referred  to  as  the  ram- 
rocket).  Accurate  data  are  needed  for  trade-off  studies  (es¬ 
pecially  in  the  medium-range  area  where  solid  propellant 
motors  can  effectively  compete)  in  which  thrust-time  char¬ 
acteristics  are  input  into  mission  analysis  codes  It  is  also 
necessary  to  standardize  the  techniques  as  much  as  possi¬ 
ble  so  that  performance  data  reported  by  one  NATO  na¬ 
tion  (or  facility)  can  be  effectively  and  fairly  compared  to 
performance  data  reported  by  other  organizations.  Each 
facility  may  employ  different  types  of  air  heaters,  chemical 
equilibrium  codes,  instrumentation,  calibration  techniques 
and  testing  methods.  Identical  motors  tested  in  different 
facilities  can  result  in  different  reported  thrust  levels  and 
combustion  efficiencies.  These  differences  become  even 
more  critical  with  the  introduction  of  metallized  fuels  and 
with  increased  flight  Mach  number. 

In  recognition  of  these  needs,  the  AGARD  Propulsion  and 
Energetics  Panel  established  Working  Group  22  to  gen¬ 
erate  a  working  document  which  delineates  the  recom¬ 
mended  methods  for  the  determination  of  connected-pipe 
ramjet  and  ducted  rocket  performance.  To  accomplish  this 
goal  the  Working  Group  collected,  reviewed  and  evaluated 
the  methods  and  techniques  used  in  the  NATO  commu¬ 
nity. 

The  scope  of  this  document  restricts  itself  to  experi¬ 
mental  and  analytical  methods  for  the  determination  of 
connected-pipe,  subsonic  combustion  ramjet  (with  solid 
and  liquid  fuels)  and  ducted  rocket  internal  performance. 
In  addition  to  an  overview  of  ramjet  and  ducted  rocket 
propulsion  devices  the  following  six  major  areas  of  inter¬ 
est  are  addressed: 

1.  Recommended  methods  for  reporting  test  results,  in¬ 
cluding  the  methodology  for  uncertainty  analysis. 

2.  Explanation  of  the  methods  used  for  the  prediction 
of  theoretical  thermodynamic  and  performance  para¬ 
meters  (codes  employed,  values  used,  input  data  re¬ 
quirements,  etc.). 

3.  Delineation  of  the  parameters  required  to  be  mea¬ 
sured. 

4.  Explanation  of  the  calculation  methods  for  experi¬ 
mental  performance  parameters. 


.5.  Sample  performance  calculations  for  each  of  the 
propulsion  devices  and  determination  of  the  sensi¬ 
tivity  of  the  various  performance  parameters  to  the 
measured  variables.  Also  included  is  an  example  of 
uncertainty  analysis. 

6.  Techniques  for  the  utilization  of  air  heaters  and  the 
effects  of  air  heaters  on  theoretical  and  experimental 
performance. 
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2  Overview 


After  the  end  of  the  second  world  war,  major  research  ef¬ 
forts  were  undertaken  in  several  countries  on  supersonic, 
airbreathing  propulsion.  This  led  to  numerous  experi¬ 
mental  missile  and  aircraft  flight  tests,  for  example  the 
French  GRIFFON  aircraft,  and  to  a  few  operational  sys¬ 
tems.  Among  the  latter  can  be  cited  a  first  generation 
of  airbreathing  missiles,  such  as  the  American  BOMARC 
and  TALOS,  the  British  BLOODHOUND  SEA  DART 
(Fig.  2.1)  and  the  Soviet  SA^.  In  the  seventies  a  second 
generation  began  to  appear  with  new  technologies,  princi¬ 
pally  the  ducted  rocket  and  the  integral  booster.  The  So¬ 
viet  SA6,  a  ducted  rocket  with  an  integral  booster,  showed 
the  effectiveness  of  the  new  design  in  the  field,  in  Middle 
East  conflicts.  In  France,  the  ASMP  missile,  a  liquid  fuel 
ramjet,  has  been  deployed  by  the  aircraft  of  the  strategic 
forces,  since  1986  (Fig.  2.2). 

Today,  the  ramjet  is  drawing  attention  again,  throughout 
the  world,  for  potential  military  (tactical  and  strategic) 
and  civilian  applications; 

•  For  missiles; 

France  and  Germany  are  preparing  the  ANS  super¬ 
sonic  antiship  missile  (Fig.  2.3).  The  USA  is  working 
on  air-to-air  advanced  developments  such  as  VFDR 
and  AAAM,  etc.  There  are  also  other  cooperative 
efforts  between  several  NATO  countries. 

•  For  hypersonic  or  orbital  aircrafts  and  space 
launchers: 

Research  on  ramjet  propulsion,  with  subsonic  or  su¬ 
personic  combustion,  is  making  a  strong  comeback 
in  a  number  of  countries,  mainly  in  the  USA  and  in 
Germany  (Figs.  2.4  and  2.5),  but  the  operational  ap¬ 
plications  are  expressed  only  in  the  long  term. 


2.1  Ramjet  Configurations 

The  basic  advantages  of  all  ramjet  configurations  over  con¬ 
ventional  rocket  propulsion  systems  are  twofold.  Firstly 
they  have  the  potential  to  achieve  an  increased  range  and 
secondly  having  ”power-on-to-target”  and/or  higher  ter¬ 
minal  velocity  they  offer  increased  effectiveness  against 
manoeuvring  targets.  Either  or  both  of  these  advantages 
can  be  sufficient  to  justify  the  use  of  a  ramjet  over  a  con¬ 
ventional  rocket  motor  in  certain  applications. 

Among  the  various  configurations  which  have  been  stud¬ 


ied,  a  classification  can  be  established  according  to  the  na¬ 
ture  of  the  fuel,  either  liquid  with  its  high  performance,  or 
solid  with  its  operational  simplicity  and  potentially  lower 
cost. 


2.1.1  Ramjets  Using  Liquid  Fuel  (LFRJ) 

The  liquid  fuel  ramjet  can  use  classical  kerosene,  high  den¬ 
sity  or  slurry  fuels.  The  liquid  fuel  ramjet  (Fig.  2.6  and 
2.7)  dominates  the  operational  applications,  mainly  be¬ 
cause  of  its  high  throttleability  and  excellent  performance. 


2.1.2  Ramjets  Using  Solid  Fuel  (SFRJ) 

It  is  possible  to  use  special  solid  fuels  in  a  solid  fuel  ramjet 
in  order  to  obtain  conditions  of  maintenance  and  storage 
similar  to  that  of  ordinary  ammunitions  or  classical  solid 
propellant  rocket  missiles  (Fig.  2.8).  The  engine  uses  only 
one  chamber,  resulting  in  a  very  simple  construction.  The 
pure  solid  fuel,  that  is  without  oxidizer,  covers  the  wall 
of  the  combustor.  By  ablation  in  the  hot  air  flow,  it  is 
transformed  into  gases  which  burn  in  a  combustion  cham¬ 
ber.  It  is  particularly  well-suited  for  the  high-acceleration 
environment  of  projectiles. 


2.1.3  Ducted  Rockets  (DR) 

The  ducted  rocket  contains  a  solid  propellant  with  only  a 
small  portion  of  oxidizer  (fuel  rich  solid  propellant).  Just 
the  quantity  of  oxidizer  is  used  which  is  necessary  to  pro¬ 
duce  gases  through  pyrolyzing  and/or  burning  reactions. 
The  ducted  rocket,  like  the  solid  fuel  ramjet,  has  the  main¬ 
tenance  and  storage  characteristics  of  a  solid  rocket  mo¬ 
tor.  Like  the  liquid  fuel  ramjet,  the  ducted  rocket  may 
also  have  a  throttling  ability. 

Two  variants  of  the  ducted  rocket  exist; 

1.  Ducted  rocket  with  separate  gas  generator 
(Figs.  2.9,  2.10  and  2.11) 

The  fuel  is  stored  in  a  separate  container,  or  gas  gen¬ 
erator,  which  works  like  a  rocket.  The  gases  pro¬ 
duced,  relatively  low  in  temperature,  can  be  injected 
into  the  combustion  chamber  through  a  control  valve. 


Figure  2.2;  ASMP  missile 


lir  intake  fkiel  it^ector  and 


oombustkm  nooie 
chamber 


Figure  2.6;  Sketch  of  a  liquid  fuel  ramjet 
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Figure  2,7;  Si/lTmissile 
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Figure  2.8:  Sketch  of  a  solid  fuel  ramjet 
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Figure  2.9:  Sketch  of  a  ducted  rocket  (with  separate  gas  generator) 


Figure  2.10:  Projectile  with  fuel  rich  solid  propellant  ducted  rocket  propulsion 


As  the  burning  rate  of  the  fuel  is  influenced  by  pres¬ 
sure,  it  is  possible  to  regulate  the  gas  flow. 


2.  Ducted  rocket  with  integrated  gas  generator 
(Fig.  2.12,  2.13) 

An  example  is  the  original  French  design,  called  “Rvs- 
tique”.  It  has  a  single  chamber,  fuel  rich  solid  propel¬ 
lant  in  direct  contact  with  the  combustion  chamber, 
and  wide  altitude  variation  capability  because  of  self¬ 
regulation  qualities  (as  the  burn  rate  is  pressure  de¬ 
pendent). 


2.1.4  Scramjets  (Supersonic  Combustion) 

The  ramjet  is  without  any  doubt  the  most  suitable  air- 
breathing  propulsion  system  for  hypersonic  flight  in  the 
atmosphere.  Efficient  operation  of  a  ramjet  is  reached  by 
subsonic  combustion  up  to  Mach  6  or  7,  and  supersonic 
combustion  beyond;,  in  the  latter  case,  the  engine  is  des¬ 
ignated  as  a  scramjet. 

Theoretically,  on  the  basis  of  its  power  performance,  it 
could  reach  orbital  velocities.  In  practice,  it  will  be  diffi¬ 
cult  to  go  beyond  a  hypersonic  speed  of  about  Mach  10  to 


7 


Figure  2.11:  BFA  ducted  rocket  missile 


Figure  2.12.  Sketch  of  “Rustique"  ducted  rocket  missile 


12,  because  of  the  sensitivity  of  the  engine  thrust  to  small 
disturbances  at  higher  Mach  numbers.  However,  at  the 
present  time,  the  complete  operation  envelope  remains  to 
be  explored  at  the  cost  of  considerable  research  and  devel¬ 
opment  efforts.  Several  countries  have  begun  to  consider 
the  potential  use  for  the  scramjet  and  have  launched  am¬ 
bitious  programmes  to  develop  demonstrators. 


2.2  Ramjet  Performance  Deter¬ 
mination 

It  is  more  difficult  to  evaluate  the  performances  of  air- 
breathing  engines  than  those  of  rockets,  because  they  vary 
strongly  with  the  flight  conditions  (Mach  number,  alti¬ 
tude,  atmospheric  conditions,  angle  of  attack,  etc). 

The  main  difficulties  are  discussed  below, 

1,  Required  precision  of  performance  determina- 
tioii 

Ramjet  engine  net  thrust,  available  for  vehicle  propul- 
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Figure  2.13:  “Rusiiqut”  ducted  rocket  missile 


sion  is  equal  to  the  thrust  generated  bj'  the  nozzle 
minus  the  ram  drag,  This  ram  drag  is  equal  to  the 
momentum  flux  of  the  incoming  air  flow,  and  is  also 
called  inlet  momentum  drag. 

At  high  flight  speeds  the  nozzle  thrust  and  the  inlet 
momentum  drag  will  have  the  same  order  of  magni¬ 
tude.  Thus,  an  error  in  the  nozzle  thrust  will  prop¬ 
agate  into  an  error  in  the  net  thrust.  This  increases 
with  flight  Mach  number  as  shown  in  Fig.  2,14.  For 
instance,  at  a  flight  speed  of  Mach  4,  a  1  percent  error 
in  nozzle  thrust  may  lead  to  a  3  percent  error  in  net 
thrust  or  3  percent  in  range. 


2.  Engine  airframe  integration 

Especially  for  high  Mach  numbers,  the  different  com¬ 
ponents  of  the  engine,,  mainly  the  air  intake  and  the 
nozzle,  are  integrated  in  the  aerodynamic  configura¬ 
tion  of  the  vehicle.  It  is  therefore  difficult  to  separate 
the  thrust  and  drag  terms. 


3.  Specific  problems  for  ramjets  using  solid  fu¬ 
els/propellants 

It  is  more  difficult  to  know  the  combustion  efficiency 
of  a  solid  fuel  ramjet  or  a  ducted  rocket  than  that 
of  a  liquid  fuel  ramjet,  because  of  the  difficulty  of 
mass  flow  measurement,  and  sometimes  because  of 
the  presence  of  condensed  material  on  the  nozzle 
and/or  in  the  exhaust.. 


2.3  Main  Stages  of  Ramjet  Devel¬ 
opment 

As  for  any  propulsion  system,  developing  a  ramjet  engine 
goes  first  through  successive  development  phases,  then 
through  detailed  debugging  and  demonstration,  and  fi¬ 
nally,  acceptance  testing  under  all  flight  conditions.  This 
demands  a  great  deal  of  experimental  research  and  devel¬ 
opment. 

For  example,  to  debug  one  current  operational  ramjet 
powered  missile  and  to  qualify  it  with  its  equipment  under 
all  flight  conditions,  600  test  runs  were  required  each  year 
for  seven  years  (90%  of  blowdown  tests  lasted  30  seconds 
and  10%  lasted  longer),  using  some  80,000  kg  of  liquid 
fuel! 

In  the  same  way  as  for  aircraft,  the  current  trend  is  to 
qualify  the  missile  on  the  ground,  in  the  most  realistic 
environment  possible,  so  that  the  flight  tests,  always  very 
costly,  have  a  high  probability  of  success. 

During  a  ramjet  development,  the  following  successive  ex¬ 
perimental  steps  are  needed. 

1.  Design  test  on  components 

This  includes  tests  of  air  intakes  in  a  vvind  tunnel, 
optimisation  of  the  combustor  with  the  help  of  flow 
visualisation  techniques,  etc. 

2.  Connected-pipe  tests 
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Flight  Mach  Number 


Figure  2.14;  Sensitivity  of  net  thrust  to  a  1%  change 
in  nozzle  thrust 


This  is  a  very  important  step  of  a  development,  and 
the  heart  of  this  report,  The  simulation  of  flight  con¬ 
ditions  {velocity,  altitude)  is  obtained  with  the  help 
of  different  devices,  such  as  the  air  heater  (heat  accu¬ 
mulation  or  fuel  combustion).  The  engine  is  supplied 
with  subsonic  hot  air,  with  simulation  of  the  aero¬ 
dynamic  conditions  at  the  end  of  the  inlet  diffuser. 
An  important  problem  is  the  quality  of  the  air  enter¬ 
ing  the  ranyet,  for  instance  the  percentage  of  water 
vapour.  With  the  use  of  special  devices,  to  eliminate 
or  to  compensate  for  the  inlet  momentum,  it  is  pos¬ 
sible  to  determine  the  ramjet  nozzle  thrust. 

3.  Semi-free  jet  tests 

The  engine,  including  air  intakes,  is  supplied  with  su¬ 
personic  air  coming  from  nozzles  just  in  front  of  each 
inlet.  The  air  mass  flow  required  is  roughly  50%  to 
100%  higher  than  in  connected-pipe  tests,  due  to  ex¬ 
ternal  air  flow. 

4.  Free  jet  tests 

This  is  the  best  simulation,  because  the  entire  vehicle 
forebody  is  surrounded  with  supersonic  air  flow,  as  in 
flight.  However,  a  free  jet  test  installation  needs  to 
be  very  powerful  and  its  cost  is  very  high.  Therefore, 
free  jet  testing  is  not  always  employed,  resulting  in 
higher  risks  during  flight  tests. 

Facility  size  can  limit  or  even  preclude  free  jet  test¬ 
ing  with  a  full  scale  forebody/inlet  installed.  Spe¬ 
cial  test  techniques  are  occasionally  utilized  to  cir¬ 
cumvent  free  jet  size  or  .flow  rate  limitations.  Two 
prevalent  techniques  are  the  forebody  simulator  and 
jet  stretcher,  The  forebody  simulator  (a  contoured 
forebody  which  may  be  half  the  length  of  the  com¬ 
plete  forebody)  provides  the  same  inlet  flow  field  in 
the  free  jet  environment  as  the  complete  forebody  in 
the  flight  environment.  The  jet  stretcher  (an  aero- 
dynamically  shaped  surface  which  simulates  a  free  jet 
streamline.  Fig.  2.15)  extends  the  test  rhombus  of  the 
free  jet  nozzle  by  precluding  extraneous  shock  or  ex¬ 


Figure  2.15:  Jet  stretcher  concept 


pansion  waves  from  the  nozzle  or  jet  boundary  from 
being  reflected  into  the  flow  upstream  of  the  inlet, 

5.  Flight  tests 

This  is  the  final  objective  after  several  years  of  ground 
tests. 

2.4  Need  for  Standardization  in 
NATO  Nations 

In  summary,  it  is  seen  that: 

•  the  ramjet  configurations  can  be  various, 

•  the  determination  of  performance,  with  good  preci¬ 
sion,  is  not  easy, 

•  the  ramjet  test  facilities  are  complex  and  costly,  es¬ 
pecially  to  simulate  high  Mach  numbers,  and  can  use 
many  different  techniques. 

Many  exploratory  development  programmes  are  carried 
only  through  connected-pipe  testing.  Connected-pipe 
tests  also  constitute  the  initial  and  a  major  portion  of  any 
ramjet  missile  propulsion  development  programme.  The 
data  obtained  from  these  tests  can  be  used  to  validate 
the  designs  of  components  such  as  inlets,  combustors  and 
exhaust  nozzles  and  can  also  be  used  for  some  system  in¬ 
tegration  (combustor-inlet  coupling,  etc.).  The  resulting 
data  can  also  be  used  for  preliminary  trajectory  analysis 
for  the  missile  system.  Thus,  there  is  a  need  for  compar¬ 
ison  and  recommended  procedures  for  the  connected-pipe 
methods  used  by  the  various  NATO  nations.  This  is  the 
aim  of  this  report. 


to 


3  Methods  for  Reporting  Test  Results 


Established- reporting  standards  used  by  the  participat¬ 
ing  nations  were  reviewed.  These  included  national  stan¬ 
dards  and  those  used  by  individual  organizations.  These 
standards  were  used  as  a  basis  for  the  recommendations 
presented  in  this  chapter.  These  include: 

1.  Identification  of  vehicle  stations 

2.  General  test  information  to  be  reported 

3.  Geometric  data  to  be  reported 

4.  Description  of  equipment  and  instrumentation 

5.  Test  data  to  be  reported,  as  appropriate 

6.  Performance  data  to  be  reported 

7.  Error  analysis 

The  application  of  Sl-units  is  required  for  reporting. 

3.1  Identification  of  Vehicle  Sta¬ 
tions 

Station  locations  for  flight  conditions  were  specified  as  fol¬ 
lows: 

oo  Freestream. 

0  Flow  field  immediately  upstream  of  the  inlet  shock 
system 

1  Inlet  lip  cross-section  or  capture  station;  beginning  of 
the  internal  flow  system 

2  End  point  of  the  inlet  compression  process  or  end  of 
the  inlet  diffuser 

3  An  appropriate  position  at  the  upstream  end  of  com¬ 
bustor  section 

4  Downstream  end  of  combustor  section 

5  Exhaust  nozzle  georhetric  throat 
Nozzle  exit 


Figure  3.1:  LFRJ  configuration  station  numbers 


Figure  3.2:  SFRJ  configuration  station  numbers 


Figure  3.3:  Ducted  rocket  configuration  station  numbers 


For  connected-pipe  test  installations  these  station  num¬ 
bers  should  be  adopted.  Schematic  diagrams  are  shown 
in  Figures  3.1  to  3.3  which  illustrate  the  application  of 
station  numbers  to  configurations  which  are  typical  of  liq¬ 
uid  fuel  ramjets,  solid  fuel  ramjets  and  ducted  rockets. 
Subnumbers  can  be  used  to  describe  intermediate  stations 
downstream  of  the  primary  stations  defined  above.  For 
example,  stations  along  the  combustor  would  be  identi¬ 
fied  as  3.1,  3.2,  etc.  Secondary  inlets,  such  as  the  bypass 
inlet  shown  on  the  solid  fuel  ramjet  in  Figure  3.2,  should 
be  identified  with  an  asterisk  (*),  with  further  description 
as  deemed  necessary  by  the  author. 
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o  Determining  combustion  efficiency 
o  Determining  pressure  losses 
o  Determining  fuel/propellant  mass  flow  rate 
o  Time  averaging 
o  Area  averaging 

•  Special  measurements,  such  as  gas  sampling,  spec¬ 
troscopy,  etc,  should  be  described  as  welt, 


3.3  Geometric  Data  to  be  Re 
ported 


Figure  3.4:  Example  of  station  letters  for  multiple  inlets 

For  multiple  inlets  at  a  common  combustor  station  or  mul¬ 
tiple  nozzle  configurations,  station  numbers  should  be  fol¬ 
lowed  by  capital  letters  which  distinguish  the  individua! 
components,  beginning  with  "A”  at  the  first  component 
to  be  identified  from  the  top  of  the  vehicle  cross-section, 
proceeding  in  a  clockwise  direction  when  looking  upstream 
from  the  exhaust  nozzle  (Fig,  3.4), 

3.2  General  Test  Information  to 
be  Reported 

The  information  listed  beneath  each  topic  below  should 
be  reported  and  described,  as  appropriate. 


•  A  schematic  representation  to  describe  the  geometry 
of  the  experimental  hardware,  as  required,  including 
appropriate  station  numbers 

•  Geometric  flow  area  or  dimensions  characterizing  the 
flow  area  for  the  applicable  stations,  including  area 
variations  during  the  test,  if  any 

•  Corrective  coefficients  for  flow  area,  if  used  for  per¬ 
formance  calculations,  such  as  the  nozzle  discharge 
coefficient 

•  Air  injection  configuration  (coaxial,  multiple  inlets, 
side  dump,  etc,) 

•  Air  injection  angle  relative  to  the  ramcombustor  cen¬ 
ter  axis 

•  Fuel  injection  configuration 

•  Fuel  injection  angle,  relative  to  the  ramcombustor 
center  axis 


•  Fuel/Propellant: 

o  Chemical  formula(s)  representing  the  composi- 
■  tion 

o  Heat  of  formation  at  298,15 /\ 
o  Enthalpy  at  initial  temperatures 
o  Stoichiometric  fuel/air  ratio 
o  Density  at  298,15  A' 


Exhaust  nozzle  internal  geometry 

Afterbody  and  external  nozzle  geometry,  as  appropri- 
ate 

Axial  dimensions  as  required 

Type  of  ramburner  construction  and  materials,  in- 
eluding  type  and  extent  of  thermal  protection 


•  Air: 

o  Chemical  formula  representing  air 
o  If  vitiated  air  is  used,  also  include: 

*  Fuel  used  for  the  vitiator  (fuel  defined  as 
above) 

*  Oxygen  make-up  (vitiator  oxidizer  and  re¬ 
plenishment) 

*  Composition  of  vitiated  air 

•  Corrections  and  assumptions  made  to  arrive  at  the" 
final  test  results  may  include  methods  for: 


3.4  Description  of  Equipment  and 
Instrumentation 

While  it  is  generally  accepted  that  information  about  the 
related  test  equipment  and  instrumentation  is  important 
to  properly  document  an  experimental  investigation,  it  is 
not  essential  for  the  purpose  of  determining  ramjet  perfor¬ 
mance,  It  is,  however,  considered  essential  to  report  the 
type  and  location  of  all  measurements  required  for  perfor¬ 
mance  calculations.  The  locations  should  be  indicated  at 


12 


each  relevant  station  or  substation  by  a  sketch  or  written 
description. 

It  is  assumed  that  the  instrumentation  and  the  associated 
treatment  of  the  data  comply  with  the  state  of  the  art  for 
obtaining  the  required  physical  properties.  For  example,  it 
is  assumed  that  thermocouple  readings  are  corrected  for 
recovery  factor  and  radiation  effects,  if  necessary  for  an 
accurate  temperature  measurement. 

3.5  Test  Data  to  be  Reported,  as 
Appropriate 

•  Mass  flow  rates  for  air,  vitiator  fuel,  make-up  oxygen, 
and  ramcombustor  fuel  or  propellant 

•  Actual  compositions,  such  as  air,  vitiated  air  and  fu¬ 
els,  if  measured 

•  Inlet  profiles,  such  as  total  pressures,  velocities  and 
temperatures,  if  known 

•  Equivalence  ratio 

•  Burning  time 

•  Fuel/propellant  initial  temperature 

•  Fuel  regression  rate 

•  Force  measurement,  corrected  for  external  effects 

•  Measured  nozzle  throat  area  variations  during  test 

•  Mach  numbers  at  Stations  2,  3  and  4,  including  a 
description  of  the  method(s)  of  calculation 

•  Static  and  total  pressures  and  temperatures  at  Sta¬ 
tions  2,  3  and  4,  Indicate  if  measured  or  describe 
method  of  calculation 

•  Specific  impulse  f*  or  actual  characteristic  velocity  c* 

at  Station  5  ,, 

•  Measured  hardware  temperature,  if  used  in  perfor¬ 
mance  calculations 

•  Measured  values  used  in  heat  loss  calculations,  if  ac¬ 
complished 

•  Amplitudes  and  frequencies  of  pressure  oscillations, 
if  recorded  Identify  measurement  locations  wuth 
station  number  and  include  modes  of  oscillation,  if 
known 

•  Exhaust  gas  composition,  if  known 

•  Other  measurements,  as  appropriate 

•  General  comments,  to  include  descriptions  of  items 
unique  to  the  experiment  which  may  not  be  well 
known  to  others 


3.6  Performance  Data  to  be  Re¬ 
ported 

•  Combustion  efficiencies 

•  Expulsion  efficiency 

•  Total  pressure  losses 

•  Nozzle  efficiency 

•  Isentropic  exponent  7p  (refer  to  appendix  B) 

•  Rich  and  lean  flammability  limits,  if  obtained 

•  Stability  limits,  if  encountered 

•  General  remarks  regarding  items  which  could  affect 
performance 

3.7  Uncertainty  Analysis  Meth¬ 
odology 

All  test  data  have  errors  or  inaccuracies.  A  means  for 
quantifying  these  inaccuraweies  is  identified  in  this  section 
and  discussed  in  Appendix  A.  The  accepted  practice  in 
the  technical  community  is  to  express  such  measurement 
inaccuracies  as  an  “uncertainty”  which  is  obtained  by  an 
uncertainty  or  error  analysis.  Error  analysis  quantifies  the 
uncertainty  for  test  data  and  serves  cis  an  invaluable  en¬ 
gineering  tool  in  the  tasks  of  designing  measurement  sys¬ 
tems,  ensuring  compliance  to  data  accuracy  requirements, 
and  interpreting  test  results  (e.g.,  AGARD-PEP  Work¬ 
ing  Group  15,  Uniform  Engine  Test  Programme  ([Ij  and 
[2))).  The  methodology  was  developed  by  R.  B.  Aber- 
nethy  and  J.  W.  Thompson  [3]  and  is  used  by  the  In¬ 
ternational  Standards  Organization  (4),  by  the  American 
National  Standard  Institute  and  American  Society  of  Me¬ 
chanical  Engineers  [5],  and  by  the  Instrument  Society  of 
America  [6). 

The  methodology  used  to  determine  test  data  uncertainty 
is  based  on  quantifying  elemental  errors  in  the  measur¬ 
ing  system  for  each  Basic  Measurement  (e.g.,  pressure, 
temperature,  force,  length  and  time),  classifying  the  er¬ 
rors  into  two  categories,  either  as  precision  (random  or 
scatter  error)  or  bias  (fixed  or  offset  error),  and  propagat¬ 
ing  the  errors  by  using  Influence  Coefficients  (determined 
from  Sensitivity  Analysis)  into  an  estimate  of  uncertainty 
for  Performance  Parameters, 


3.8  Error  Analysis 

An  error  analysis  provides  a  detailed  layout  for  error  book¬ 
keeping  and  auditing  which  yields  uncertainties  for  the  Ba- 


Pressure; 
Temperature: 
Fuel  flow; 

Area: 

Force  (load  cell): 


±0.2%  to  ±0.5%  of  reading 
±0.2%  to  ±0.6%  of  reading 
±0.4%  to  ±1.0%  of  reading 
±0,4%  to  ±0.2%  of  reading 
±0.2%  to  ±0.5%  of  reading 


Table  3.1:  Typical  measurement  uncertainties  in  ground 
test  facilities 


Thrust; 

Airflow: 

Total  Pressure  Loss: 
Combustion  Efficiency: 


±0.2%  to  ±1.0%  of  calc.  val. 
±0.4%  to  ±1.0%  of  calc,  val, 
±0.5%  to  ±2.0%  of  calc,  val. 
±2.0%  to  ±5.0%  of  calc.  val. 


Table  3.2:  Typical  performance  parameter  uncertainties 


sic  Measurements  and  Performance  Parameters,  and  iden¬ 
tifies  the  contribution  of  each  error  source  to  the  total 
uncertainty  level,  Typical  measurement  uncertainties  in 
ground  test  facilities  are  given  in  Table  3.1. 

Representative  detailed  layouts  (or  audits)  of  error  sources 
for  the  Basic  Measurements  are  presented  in  Appendix  A 
(Tables  A.l  through  A. 4,  respectively).  These  error  ana¬ 
lyses  results  and  Influence  Coefficients  (determined  by  a 
Sensitivity  Analysis  using  equations  which  calculate  the 
Performance  Parameters)  provide  an  error  propagation 
basis  (i.e.,  bookkeeping  layout)  which  yields  the  uncer¬ 
tainties  of  Performance  Parameters.  A  representative  lay¬ 
out  for  any  specific  Performance  Parameter  at  a  selected 
test  condition  is  presented  in  Appendix  A  (Table  A. 5). 
Typical  Performance  Parameter  uncertainties  are  given  in 
Table  3.2. 

Measurement  uncertainties  are  strongly  influenced  by 
ramjet  operating  conditions  (i.e.,  altitude,  Mach  number 
and  power  setting)  and  by  test  goals,  resources  and  sched¬ 
ules.  It  is  important  to  note  that  the  above  uncertainty 
values  represent  a  range  from  the  best  uncertainties  that 
can  be  achieved  to  values  that  can  be  obtained  using  every¬ 
day  measurement  practices.  Of  course,  much  larger  values 
of  uncertainty  will  result  if  any  part  of  the  measurement 
process  is  carelessly  executed. 
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4  Theoretical  Performance  Determination 


The  theoretical  performance  determination  of  ramjets  is 
based  on  an  ideal  combustion  process  which  assumes 
chemical  equilibrium.  Computer  codes  that  were  origi¬ 
nally  developed  to  evaluate  the  rocket  motor  combustion 
process  are  used  to  model  the  ramjet  combustion  process. 
The  codes  however  are  based  on  a  flow  velocity  of  zero 
in  the  combustion  chamber,  since  in  a  (nozzled)  rocket 
engine  the  Mach  number  is  small.  However,  in  a  ramjet 
combustor  substantial  flow  velocities  may  be  encountered 
(up  to  about  Mach  0.8).  This  requires  the  use  of  stagna¬ 
tion  flow  properties  as  input  to  the  codes  instead  of  static 
flow  conditions.  Some  features  of  these  codes  and  their  ap¬ 
plication  to  the  ramjet  combustion  and  vitiated  air  heater 
processes  are  discussed  in  this  chapter. 

Two  codes  are  most  frequently  used  for  performance  deter¬ 
mination.  Codes  based  on  NASA  CET89  [7,  8, 9]  are  in  use 
by  organisations  in  all  countries  while  the  Naval  Air  War¬ 
fare  Center  Weapons  Division  (formerly  Naval  Weapons 
Center)  Propellant  Evaluation  Program  (PEP)  is  a  code 
[10]  used  by  many  organisations  in  the  United  States. 

Test  cases  were  developed  to  compare  thermochernical 
properties  of  ramjet  reacting  flows  predicted  by  these 
•codes  as  applied  by  different  organisations.  The  differ¬ 
ences  in  code  outputs  are  discussed. 

Combustor  performance  can  be  affected  by  combustion  in¬ 
stability  and  heat  loss  and  may  also  be  driven  by  plume 
signature  requirements.  Several  references  are  given  in  the 
bibliography  [11]-[14]  that  are  used  to  predict  combustion 
instability  and  plume  signature.  Two-phase  flow  is  ad¬ 
dressed  by  some  of  these  codes. 


4.1  List  of  Theoretically  Deter¬ 
mined  Parameters  for  Perfor¬ 
mance  Calculations 

Many  parameters  were  identified  that  may  be  required  in 
order  to  calculate  performance.  Some  of  these  parameters 
are  difficult  to  measure  or  determine  experimentally  and, 
therefore,  are  theoretically  calculated.  They  are  listed  and 
described  below: 


Station  2; 

•  Molecular  Weight  (M2) 

•  Mole  Fraction  of  Species  (xz) 

•  Isentropic  Exponent  (72) 

Station  4; 

•  Stagnation  Temperature  (T14) 

•  Molecular  Weight  (Ma) 

•  Mole  Fraction  of  Species  (^4) 

•  Isentropic  Exponent  (74)  may  be  defined  three  ways: 

o  7^4  =  Cp(T4)/cv(T4)  is  the  frozen  isentropic  ex¬ 
ponent 

o  7,4  =  — (51np/5lnu)|,4  or  — 7/(51nu/5lnp)|T4 
is  for  local  equilibrium  (Appendix  B) 

o  7p,4-5  is  a  process  isentropic  exponent  from  sta¬ 
tion  4  to  station  5  expressed  as 
ln(p4/P5)/(ln(p4/P5)  -  In(r4/T5)) 

Station  5: 

•  Static  Pressure  (ps) 

•  Static  Temperature  (T5) 

•  Characteristic  Velocity  (cj) 

•  Stagnation  Temperature  (Tts) 

•  Molecular  Weight  (Ads) 

•  Mole  Fraction  of  Species  (xs) 

•  Isentropic  Exponent  (7/5,  7,5  or  7p,4-5) 

Station  6: 

•  Static  Pressure  (ps) 

•  Static  Temperature  (Ts) 

•  Specific  Impulse  (j,p) 

•  Molecular  Weight  (Adg) 
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•  Mole  Fraction  of  Species  (xe) 

•  Isentropic  Exponent  (7^6,  7,6  or  7p,s-6) 

•  Specific  Heat  (cps) 

•  Thrust  Coefficient  (cj?) 

4.2  Description  of  Aerothermo- 
chemical  Equilibrium  Codes 

The  aerothermochemical  equilibrium  codes  which  are  in 
use  are  listed  below.  Contacts  and  addresses  for  obtaining 
these  codes  are  given  in  Table  4.1, 

1.  The  NASA  CET89  [7,  8,  9]  computer  program  is  used 
for  calculations  involving  chemical  equilibria  in  com¬ 
plex  systems  The  method  applied  is  based  on  mini¬ 
mization  of  the  Gibbs  free  energy,  The  program  per¬ 
mits  calculations  such  as: 

•  Chemical  equilibrium  for  assigned  thermody¬ 
namic  states 

o  Temperature  and  pressure 
o  Enthalpy  and  pressure 
o  Entropy  and  pressure 
o  Temperature  and  volume  or  density 
o  Internal  energy  and  volume  or  density 
o  Entropy  and  volume  or  density 

•  Theoretical  rocket  performance 

o  Frozen  flow 
o  Equilibrium  flow 

•  Chapman-Jouguet  detonations 

•  Shock  tube  parameter  calculation 

The  program  considers  condensed  species  as  well  as 
gaseous  species.  Condensed  and  gaseous  species  are 
supposed  to  have  the  same  velocity  and  temperature. 

2.  PEP  code  [10]  was  developed  for  the  calculation  of 
high- temperature  thermodynamic  properties  and  per¬ 
formance  characteristics  of  propellant  systems.  De¬ 
termination  of  chemical  equilibrium  is  accomplished 
by  a  combination  of  two  methods  [15]-[20],  An  op¬ 
timized  basis,  which  is  a  subset  of  molecular  species, 
is  chosen.  The  chemical  systeih  is  then  divided  into 
a  number  of  subsystems,  each  relating  a  nonbasis 
species  to  the  basis.  The  subsystem  with  the  greatest 
discrepancy  in  its  equilibrium  relationship  is  corrected 
stoichiometrically  until  convergence  is  obtained,  The 
program  permits  calculations  such  as; 

•  Chemical  equilibrium  for  assigned  thermody¬ 
namic  states 

o  Temperature  and  pressure 


o  Enthalpy  and  pressure 
o  Entropy  and  pressure 
b  Temperature  and  volume  or  density 
o  Internal  energy  and  volume  or  density 
o  Entropy  and  volume  or  density 

•  Theoretical  rocket  performance 

o  Frozen  flow 
o  Equilibrium  flow 

This  code  was  revised  to  run  on  a  personal  computer 
(IBM  compatible)  in  1989, 

Other  codes  that  may  be  used  are; 

1.  The  COPPELIA  code  [21]  is  based  on  the  NASA 
CET89  code  (NASA  CEC  71  version).  It  is  more 
complete  and  extended  than  the  NASA  CEC  71  code. 
It  has  limited  distribution, 

2.  The  STAN  JAN  [22]  program  is  used  to  calculate 
chemical  equilibrium  in  a  complex  system,  including 
several  phases,  The  calculation  technique  is  based 
on  the  method  of  element  potentials.  The  method 
of  element  potential  uses  theory  to  relate  the  mole 
fractions  of  each  species  to  quantities  called  element 
potentials.  There  is  one  element  potential  for  each 
independent  atom  in  the  system,  and  these  element 
potentials,  plus  the  total,  number  of  moles  in  each 
phase,  are  the  only  variables  that  must  be  adjusted 
for  the  solution.  In  large  problems,  that  is  in  cases 
with  many  species,  this  number  of  element  poten¬ 
tials  is  a  much  smaller  number  than  the  number  of 
species,  and  hence  far  fewer  variables  need  to  be  ad¬ 
justed.  The  program  assumes  that  the  gas  phase  is  a 
mixture  of  ideal  gases  and  that  condensed  phases  are 
ideal  solutions,  The  program,  called  STANJAN  be¬ 
cause  of  its  roots  at  Stanford  and  its  connection  with 
the  JANAF  thermochemical  data  tables,  is  an  inter¬ 
active  program  designed  for  use  with  either  personal 
or  mainframe  computers.  Thermodynamic  cycle  an¬ 
alysis  is  easily  executed  with  STANJAN,  because  the 
user  may  specify  the  state  parameters  in  a  variety  of 
ways  including; 

•  Temperature  and  pressure 

•  Pressure  and  entropy 

•  Enthalpy  and  pressure  same  as  last  run 

•  Volume  and  entropy  same  as  last  run 

3.  The  Aerotherm  Chemical  Equilibrium  (ACE)  [23] 
computer  program  is  a  versatile  code  for  calculating 
quantities  of  importance  for  many  thermochemical 
processes.  It  is  well  adapted  to  study  ablative  pro¬ 
cesses.  Closed  and  open  systems  (i.e,  constant  vol¬ 
ume  and  constant  pressure)  can  be  handled.  The  rel¬ 
ative  amount  of  each  chemical  element  in  the  system 
is  specified  for  closed  systems.  The  relative  amounts 
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of  chemical  elements  depend  on  various  mass  trans¬ 
fer  rates  for  open  systems.  Systems  may  be  treated 
in  chemical  equilibrium  or  certain  reactions  may  be 
kinetically  controlled.  It  has  limited  distribution. 

4.3  Application  and  Procedure 
for  Theoretical  Performance 
Determination 

Aerothermo chemical  equilibrium  codes  may  be  used  either 
interactively  with  a  data  reduction  program  or  to  gener¬ 
ate  data  tables  of  thermochemical  equilibrium  properties. 
The  tables  should  be  of  sufficient  fineness  such  that  inac¬ 
curacies  do  not  result  due  to  the  interpolation  procedure 
used.  The  codes  are  usually  used  interactively,  since  the 
computational  time  is  no  longer  a  major  concern. 

In  order  to  calculate  theoretical  performance  one  must 
define  the  inputs  to  the  aerothermochemical  equilibrium 
codes,  select  the  most  appropriate  calculation  procedure 
and  choose  the  desired  outputs.  Assumptions  must  be 
made  in  determining  theoretical  performance.  For  exam¬ 
ple,  the  gas  velocity  in  the  combustion  chamber  may  be 
assumed  to  be  zero  and  heat  losses  through  the  combus¬ 
tor  wall  as  well  as  equilibrium  (velocity  and  temperature) 
between  condensed  and  gaseous  phases  may  be  neglected. 
The  output  data  may  be  needed  for  calculation  of  other 
theoretical  performance  parameters. 

4.4  Inputs  to  the  Aerothermo¬ 
chemical  Equilibrium  Codes 

Aerothermochemical  equilibrium  codes,  such  as  NASA 
CET89  and  PEP,  require  certain  inputs.  These  inputs  con¬ 
sist  of  species  data,  pressures,  temperatures  or  enthalpies, 
mass  flow  rates  or  mass  fractions,  geometric  areas  and  in¬ 
gredient  properties. 

4.4.1  Species  Thermochemical  Data 

Species  data  are  obtained  primarily  from  the  JANAF 
Thermo  chemical  Tables  [24,  25].  The  Thermophysical 
Properties  Research  Center  (TPRC)  provides  mainte¬ 
nance  of  the  tables  [26j.  The  JANAF  tables  have  also 
been  supplemented  from  other  sources  (27]-(29]  The  tab¬ 
ulated  data  have  been  put  into  polynomial  form  for  use  by 
the  codes.  The  accuracy  of  both  the  tabulated  data  and 
the  curve  fits  to  the  data  can  significantly  affect  the  code 
output.  Later  versions  of  the  codes  contain  improved  data 
and/or  curve  fits.  Therefore,  care  must  be  taken  if  older 
versions  of  the  codes  are  utilized. 


4.4.2  Pressure  Inputs 

Requirements  are  the  stagnation  pressure  at  station  4  (M^ 
is  assumed  equal  to  zero)  and  the  pressure  ratio  between 
stations  4  and  6.  The  ambient  air  pressure  is  also  required, 
if  absolute  values  are  used.  The  stagnation  pressure  at 
station  4  is  calculated  from  the  static  properties  at  station 
4  and  the  geometric  areas  at  stations  4  and  5. 


4.4.3  Mass  Flowrates  and/or  Mass  Frac¬ 
tions 

Mass  flow  rates  are  required  to  determine  mixture  ratios 
and  for  solution  of  the  continuity  equation  at  different  sta¬ 
tions.  The  total  flowrate  for  vitiated  air  (or  flowrates  for 
air,  vitiator  fuel  and  make-up  oxygen)  and  the  flowrate  for 
ramjet  combustor  fuel  or. ducted  rocket  propellant  must  be 
measured  or  deduced  from  analysis. 


4.4.4  Geometric  Areas 

The  expansion  ratio  (As/Aj)  is  a  program  input  and  may 
be  corrected  for  boundary  layer  thickness  when  appropri¬ 
ate  (see  Section  5.2.1). 


4.4.5  Compositions  and  Temperatures  or 
Enthalpies  of  Constituents 

The  compositions  and  temperatures  or  enthalpies  of  all 
constituents  at  station  3  may  be  needed.  These  are  de¬ 
termined  from  the  constituent  temperatures  and  flowrates 
and  the  inlet  air  temperature.  Heat  loss  between  the  air 
heater  and  combustor  inlet  is  accounted  for  as  described 
in  Section  4. 4. 5. 3. 


4. 4. 5.1  Fuel/Propellant 

Fuel  properties  are  determined  from  the  NASA  CET89  in¬ 
gredient  file,  the  PEP  ingredient  file,  military  standards, 
laboratory  analyses,  handbooks  (27]-[31]  and'  manufac¬ 
turer  data  sheets.  Occasionally,  a  fuel  consists  of  sev¬ 
eral  ingredients.  Mass  weighted  calculations  of  ingredient 
properties  may  be  used  to  determine  the  properties  for  the 
fuel. 


4. 4.5.2  Ideal  Air 

The  composition  of  the  air  supplied  to  the  air  heater  is 
usually  assumed  to  be  that  of  ideal  air.  However,  for  ideal 
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code 

remarks 

address 

NASA  CET89 

NASA  CET89  is  disseminated  under 
the  sponsorship  of  the  National  Aero¬ 
nautics  and  Space  Administration  by 
the  Computer  Software  Management 
and  Information  Center  (COSMIC) 

COSMIC 

Software  Information  Services 

The  University  of  Georgia, 

Computer  Services  Annex 

Athens,  GA  30602  USA 

Program  Number:  LEW-15113 

Program  Name:  9  A  CET89 

PEP 

PEP  code  is  available,  but  distribution 
is  limited 

Naval  Air  Warfare  Center  Weapons  Division 
Code  C2776  (3276) 

China  Lake,  CA  93555 

USA 

COPPELIA 

COPPELIA  code  is  available,  but  dis¬ 
tribution  is  limited 

ONERA 

29  Avenue  de  la  Division  Leclerc 

92322  Chatillon  sous  Bagneux 

France 

STANJAN 

_ _ J 

Prof.  W.C.  Reynolds 

Dept,  of  Mechanical  Engineering 

Stanford  University 

Stanford,  CA  94305-3030 

USA 

Table  4.1:  Contacts  and  addresses  for  obtaining  aerothermochemical  equilibrium  codes 


air,  several  compositions  are  in  use  by  different  organisa¬ 
tions.  A  list  is  presented  in  Table  4.2.  The  enthalpy  of 
air  as  a  function  of  temperature  is  given  in  Table  4.3  and 
used  to  T  <  2500 A'. 

An  evaluation  of  the  effects  of  air  composition  on  theo¬ 
retical  combustor  performance  was  conducted.  Different 
ideal  airs  were  combusted  with  JP-IO  hydrocarbon  fuel. 
The  predicted  mole  fractions  for  some  species  are  listed  in 
Table  4.4.  Calculated  values  show  only  small  variations. 
Also,  the  theoretical  performance  parameters  (i.e.  T,4,  74 
and  M\)  do  not  show  significant  differences  for  the  various 
air  compositions. 


Method  1  . 

The  vitiated  air  is  assumed  to  be  ideal  air  with  a  tempera¬ 
ture  or  enthalpy  corresponding  to  the  inlet  air  temperature 
(Table  4.3)  Vitiated  air  properties  are; 

•  Flowrate  of  the  vitiated  air. 

•  Composition  of  ideal  ait. 

•  Temperature  or  enthalpy  for  ideal  ait. 

Method  2 


Ideal  air  and  vitiator  combustion  products  (assumed  com¬ 
plete)  are  input  as  oxidant  reactants  to  the  aerothermo¬ 
chemical  equilibrium  code  at  the  measured  T12  and  pt2- 
Vitiated  air  properties  are: 

4.4. 5. 3  Vitiated  Air 


The  composition  and  enthalpy  of  vitiated  ait  can  be  ob- 
tained'using  one  of  several  different  approaches  described 
below  and  depicted  in  Figures  4.2  and  4.3  ([32]).  All  meth¬ 
ods  account  for  heat  losses  that  always  occur  between  the 
vitiator  and  the  combustor  inlet  as  indicated  in  Figure 
4.1.  Chapter  8  provides  a  detailed  discussion  of  vitiated 
air  heaters.  To  is  assumed  to  be  equal  to  To,  or  the  tem¬ 
perature  at  some  intermediate  station,  since  the  la.tter  can 
be  more  easily  measured. 


•  Flowrates  of  ideal  air  and  vitiator  combustion  prod¬ 
ucts. 


•  Compositions  of  ideal  air  and  vitiator  combustion 
products. 


•  Vitiated  air  temperature  or  enthalpy,  where  enthalpy 
is  given  by  Eq.  4.1. 


fits  — 


E.  ni./ii(To  =  Tn) 
E. 


(4.1) 
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Species  | 

Composition  1 

XN 

xo  ; 

Xc 

XAt 

Ah^f,  [kJ/kg] 

N2,  O2  and  Ar 

iV83s0224^'‘S 

0.78477 

0,21053  1 

0 

0.00470 

0 

N2,  O2  and  CO2 

A54,6230i4  675C0.OIO 

0  78812 

0.21174 

0.00014  1 

0 

-4,187 

^2,  O2,  CO2  and  Ar 

^186.2041  seAro  93400.0314 

0.78443 

0.21072 

0.00016 

0.0049 

-4.187 

Table  4.2;  Composition  and  heat  of  formation  for  ideal  airs 


Temperature 

Enthalpy 

Temperature 

Enthalpy 

Temperature 

Enthalpy 

[kJ/kg] 

■K] 

[kJ/kg] 

[K] 

[kJ/kg] 

298 

0 

1050 

805 

1800 

1704 

350 

52,4 

1100  ■ 

863 

(  1850 

1766 

400 

103 

1150 

921 

1  1900 

1828 

'  ~¥o 

154 

1200 

979 

'  1950 

1890 

500 

205 

1250 

1038 

2000 

1953 

550 

257 

1300 

109? 

2050 

2015 

600 

309 

1350 

1157 

2100 

2078 

650 

362 

1400 

1217 

2150 

2141 

700 

415  ' 

1450 

1277 

2200 

2204 

750 

469  I 

1500 

1337  1 

2250 

2267 

800 

524  j 

1550 

1398 

2300 

2330 

850 

578 

1600 

■  1459 

2350 

2393 

900 

635 

1650 

1520 

•  2400 

2457 

950 

691 

1700 

1581 

2450 

2520 

1000 

748 

1750 

1643 

2500 

2584 

Table  4,3;  Air  (Aa3s0224Ar5)  enthalpy  as  function  of  air  temperature  (sea  level)  [32] 


Formulation 

Tt2,  [A'] 

Xrv2 

ACO3  ( 

XH,o 

Xco 

XAr 

XO2 

Xho 

Ns3s022iArs 

298.15 

0.729 

0.126  i 

0.108 

0.014 

0,009 

0,006  1 

0,003 

1500 

0.703 

0.084 

0.091 

0,052 

0.008 

0,020  1 

0.014 

7V54. 623^14  STsC'o  DID 

298.15 

0.736 

0.127  1 

0.108 

0,013 

0  0 

0.007  1 

0.003 

1500 

0.710 

0.085 

0.091 

0,050 

0.0 

0,021  1 

0.014 

7Vi56  204i.96Aro,934Co.031'l 

298.15 

0.729 

0.126 

0.108 

0,014 

0.009 

0.006  1 

0.003 

1  Isoo"; 

0.703 

0.084 

0.091 

^^052 

0.008 

1  0.020 

0.014 

A/r°  jp  iQ  =  -n^kj/kg,  f/a  =  0.0704,  Pi  =  21SkPa 


Table  4.4;  Some  species  mole  fractions  for  combustion  of  JP-10  with  ideal  airs. 


Ideal  Air 
Oxygen 
Vitiator  Fuel 


Vitiator 


Vitiated 

Air 


Combustor 


Tt2 


Figure  4.1:  Schematic  of  vitiator  test  setup 
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•  Flowrates  and  composi¬ 
tions  of; 


0  Ideal  air 

0  Make-uD  oxveen  ^ 

Equilibrium 

0  Vitiator  fuel 

Code 

•  T(3 


•  Composition 
of  vitiated  air 


•  P(3 


Figure  4.2:  Determination  of  flow  properties  at  station  3  assuming  equilibrium  vitiator  combustion  for  method  3 


Method  3 

Input  a  specified  temperature  (Tjs)  and  pressure  (po) 
into'  an  aero  thermochemical  equilibrium  code  and  calcu¬ 
late  the  equilibrium  flow,  given  the  mass  fractions  for  the 
air,  make-up  oxygen  and  the  vitiator  fuel  (Figure  4  2). 


Method  4 

If  the  T  —  p  approach  of  Method  3  cannot  be  done  by  the 
code  then  an  iterative  approach  can  be  used  that  gives 
identical  results.  Use  an  aerothermochemical  equilibrium 
code  to  calculate  the  composition  and  temperature  of  the 
vitiated  air  for  a  specified  pressure  (po))  given  the  mass 
fractions  for  the  air,  make-up  oxygen  and  vitiator  fuel. 
The  heat  losses  are  considered  by  adjusting  the  enthalpies 
(via  the  heat  of  formation)  of  the  vitiator  ingredients  until 
the  calculated  temperature  matches  the  measured  temper¬ 
ature  at  station  2  (Figure  4.3). 

Test  cases  were  used  to  evaluate  the  four  vitiated  air 
methods.  Air  was  burned  with  vitiator  fuels  (hydrogen 
and  methane)  to  obtain  inlet  air  temperatures  of  700  and 
1000/\.  A  lOOA'  temperature  loss  was  assumed  to  occur 
between  the  vitiator  and  the  combustor  inlet. 

The  results  of  the  vitiator  enthalpy  calculation  for  each 
method  are  tabulated  (Table  4.5).  Method  1  overpredicts 
the  value  of  enthalpy  when  compared  to  the  more  rigorous 
approaches  of  Methods  2-4.  This  was  due  to  the  failure 
to  account  for  vitiator  product  species.  The  values  of  vi¬ 
tiator  enthalpy  calculated  from  Methods  2-4  were  nearly 
identical.  It  is  recommended  that  one  of  these  methods  be 
used  for  most  applications  and  that  the  use  of  Method  1  be 
restricted  to  conditions  with  a  low  vitiator  temperature, 


4.5  Results 

from  the  Aerothermochemi¬ 
cal  Equilibrium  Codes 

The  numerical  vMues  of  some  thc’oretical  performance  pa¬ 
rameters  may  be  different  depending  on  which  equilibrium 
option  is  selected.  The  code  outputs  are  also  dependent 
on  assumptions  made  in  determining  the  inputs  or  on  the 
calculation  procedure  itself.  Each  of  these  items  is  consid¬ 
ered  further  in  the  following  subsections. 


4.5.1  Equilibrium  Option 


A  user  of  any  aerothermochemical  equilibrium  code  must 
decide  whether  frozen  or  local  equilibrium  flow  is  more 
appropriate.  Generally,  local  equilibrium  is  appropriate 
from  station  4  to  station  6.  Condensation  and  recombi¬ 
nation  processes  can  be  important  between  stations  5  and 
6.  Therefore,  from  station  5  to  station  6  one  must  decide 
between  local  equilibrium,  frozen  or  a  combination  of  local 
equilibrium  followed  by  frozen  flow,' 

However,  it  is  recommended  to  run  the  code  first  with  both 
the  local  equilibrium  and  the  frozen  flow  options  between 
stations  4  and  6  to  determine  the  difference  in  performance 
between  the  two  modes.  If  a  difference  of  more  than  5% 
is  observed  in  specific  impulse  it  is  recommended  that  the 
code  be  re-run,  assuming  local  equilibrium  flow  between 
stations  4  and  5  and  frozen  flow  between  stations  5  and  6. 
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•  Flowrates,  composition 


•  Composition 
of  vitiated  air 

•  hts 

•  Adjusted  Ah” 
of  vitiator  in¬ 
gredients 


Figure  4.3:  Determination  of  flow  properties  at  station  3  assuming  equilibrium  vitiator  combustion  for  method  4 


1  Product  enthalpies  [kJ/hg] 

Vitiator  Fuel 

Vitiator  Temperature  [/iT] 

Method  1 

Method  2 

Method  3  iz  4 

Hydrogen 

|. 

700 

415.0 

-112.8 

-112.8 

1000 

747.6 

-123.2 

-123.1 

.  Methane 

_ i 

700 

1  415.0 

-169.4 

-169.3 

1000 

1  747.6 

-216.9 

-217.1 

Table  4.5:  Summary  of  product  enthalpies  (Air  N835  0224Ars  at  pt2  =  650i(:Pa) 


4.5.2  Effects  of  Input  Parameters  on  The¬ 
oretical  Performance 

The  effects  of  combustor.  Mach  number,  condensed  species 
and  the  aerothermochemical  equilibrium  code  used  on  the¬ 
oretical  performance  parameters  were  evaluated.  The  ef¬ 
fects  of  ideal  air  versus  vitiated  air  vvere  discussed  in  Sec¬ 
tion  4. 4. 5. 3,  but  are  addressed  in  more  detail  in  Chapter 
7  and  8.  The  inputs  used  in  this  investigation  ate  listed 
below. 

T,2  =  625A'  : 

Akj  =  QJ jkg  for  Aig3502242lr5  at  298/\ 

Ah”  -  335  kJ/kg  for  at  625A" 

Ah)  =  -77UJ/kg  for  CioHis  at  298  A’ 

Ah”  =  OJ/kg  for  Boron  at  298/1 
Ah”  =  OJ/kg  for  Magnesium  at  298/<' 

Pi  =  2l8kPa 
Ps  =  lOOi  Pa 


f/a  =  0.0704  for  CioAie 

f/a  =0.1  for  (C'io//i6)50%(5)40?t(Af?)io% 

Two  ramjet  fuels  were  used  in  order  to  address  the  is¬ 
sue  of  condensed  species.  One  was  a  liquid  hydrocarbon 
fuel  (Cio/Zis)  and  the  other  a  metallized  fuel  (a  blend  of 
CioMie  (50%),  B  (40%)  and  Mg  (10%)  by  mass).  Ideal  air 
(A83s0224Ar5)  rather  than  vitiated  air  (ideal  air,  make-up 
oxygen  and  vitiator  fuel)  was  used  to  simplify  the  calcu¬ 
lation  procedure.  The  results  are  presented  in  Table  4.6. 

The  codes  require  as  an  input,  but  P4  is  usually  the 
only  available  measurement.  If  it  can  be  assumed  that 
Mi  =  0  then,  p(4  =  p^.  However,  if  this  assumption  is  not 
valid,  then  it  becomes  important  to  calculate  P(4  and  use 
this  value  as  the  input  to  the  code.  Specific  impulse  is  very 
dependent  on  the  value  used  for  pa  as  seen  from  cases  1 
and  4  of  Table  4.6.  Combustion  temperatures  were  some¬ 
what  affected  by  pressure  and  other  parameters  were  only 
slightly  affected.  The  validity  of  assuming  Pi4  =  p4  is  de¬ 
pendent  on  the  parameters  of  interest  and  the  combustor 
Mach  number. 
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Case 

Fuel 

M4 

Pt4 

kPa 

T,4 

K 

Ts 

^Sp 

Ns/kg 

Code 

1 

Liquid 

0 

218 

2449 

2271 

1022 

PEP 

2 

Liquid 

0 

218 

2450 

2269 

1022 

NASA  CET89 

3 

Liquid 

0.41 

240 

2453 

2272 

1072 

PEP 

4 

Liquid 

0.86 

330 

2465 

2280 

1247 

PEP 

!  5 

Metallized 

0 

218 

2709 

2540 

1052 

PEP 

6 

Metallized 

0.41 

240 

2714 

2544 

1112 

PEP 

7 

Metallized 

0.86 

330 

2733 

2556 

1287 

PEP 

8 

Metallized 

0.86 

330 

2736 

2558 

1288 

NASA  CET89 

Tkble  4.6:  VaJues  of  theoretical  performance  parameters 


There  were  no  significant  differences  (compare  cases  1  and 
2  or  cases  7  and  8  of  Table  4.6)  in  the  calculated  values  ob¬ 
tained  from  the  two  equilibrium  codes  used  (NASA  CET89 
and  PEP). 


method  will  give  a  result  identical  to  method  2,  if  the 
actual  mole  fraction  of  oxygen  in  the  vitiated  air  is  0.2095. 

The  NASA  CET89  code  also  provides  the  equivalence  ratio 
from  which  the  stoichiometric  fuel/air  ratio  can  be  deter¬ 
mined. 


4.6  Determination  of  the  Stoi¬ 
chiometric  Fuel/Air  Ratio 


The  stoichiometric  fuel/air  ratio  is  required  for  calculat¬ 
ing  equivalence  ratio  (Equation  4.2)  and  is  generally  de¬ 
fined  by  the  stoichiometric  equation  for  ideal  air  and  the 
fuel/propellant  assuming  complete  combustion,  regardless 
of  whether  the  actual  air  is  vitiated  or  not. 

0  _  (//‘^)meaj 
(/ /<2)}loicfc 


4.7  other  Aspects  of  Theoretical 
Performance  Prediction 


Several  computer  codes  are  available  for  the  prediction 
of  combustion  instabilities  [11]-[13],  plume  signature  and 
multi-phase  flow  losses  [14].  Any  modern  finite  element 
code  will  calculatepotential  acoustic  frequencies  that  may 
be  excited.  There  are  no  codes  for  the  prediction  of  fre¬ 
quencies  that  will  be  excited. 


Other  methods  are  sometimes  used  to  determine  this  ratio, 
however,  since  the  air  is  almost  always  vitiated.  Method 
2  considers  only  the  ideal  air  portion  of  the  vitiated  air 
to  determine  a  stoichiometric  fuel/air  ratio.  Combustion 
products  of  the  vitiator  are  not  included. 

Method  3  requires  the  determination  of  a  theoretical 
equivalence  ratio  for  a  hydrocarbon  ramjet  fuel  (the  ra¬ 
tio  of  oxidizer  required  for  complete  combustion  of  the 
vitiator  and  ramjet  fuels  to  that  available  for  combustion) 
according  to  Equation  4.3.  The  numerator  consists  of  the 
oxidizer  used  in  the  vitiator  combustor  process  subtracted 
from  that  required  for  complete  combustion  of  all  fuels  at 
station  4.  The  denominator  is  the  oxidizer  available  for 
combustion  of  the  ramjet  fuel.  The  oxidizer  needed  for 
complete  combustion  of  the  vitiator  fuel  is  appropriately 
accounted  for  by  this  method. 

(iVtf/2.0d-Arc  X  2.0)4 -(NH/2.0-bNc  x  2.0); 

“  {No)2  -  (Nh/^  O  +  Nc  X  2.0)2 

(4.3) 

The  stoichiometric  fuel/air  ratio  may  then  be  calculated 
using  Equation  4.2  and  the  measured  fuel/air  ratio.  This 
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5  Required  Measured  Parameters 


The  main  objective  of  this  chapter  is  to  specify  the  mini¬ 
mum  set  of  parameters  that  should  be  measured  to  deter¬ 
mine  the  performance  of  a  ramjet,  engine.  Typical  ways 
of  measuring  these  parameters  are  indicated.  Some  ad¬ 
ditional  nonessential  (but  useful)  data  that  can  also  be 
obtained  are  briefly  mentioned. 

5.1  Required  Measured  Parame¬ 
ters  to  Evaluate  Ramjet  En¬ 
gine  Performance 

The  measured  parameters  listed  below  use  the  station 
numbers  and  nomenclature  specified  in  Chapter  3  and  are 
valid  for  the  liquid  fuel  ramjet,  ducted  rocket  and  solid  fuel 
ramjet.  The  parameters  are  listed  according  to  whether 
they  are  measured  before  and/or  after  a  test  or  during  a 
test. 


5.1.1  Measurements  Taken  Before  and/or 
After  a  Test 

These  measurements  are  the  following; 

A2  geometric  area  at  station  2 

geometric  area  at  station  2*  for  solid  fuel  ramjet 
with  by-pass  inlets 

A^  geometric  area  at  station  4 

As  geometric  area  at  station  5 

yle  geometric  area  at  station  6;  if  used 

Ai  nozzle  base  area  (Figure  5.1) 

Ag  ducted  rocket  gas  generator  throat  area 

m;,,-  initial  mass  of  fuel/propellant  (SFRJ  and  DR  only) 

mj  j  final  mass  of  fuel/propellant  (SFRJ  and  DR  only) 

7}  ,■  initial  fuel/propellant  temperature 


Figure  5.1:  Base  area  for  a  convergent-divergent  nozzle 

5.1.2  Measurements  Taken  During  a  Test 

Each  parameter  is  measured  as  a  function  of  time.  It  is 
assumed  here  that  inlet  conditions  are  measured  at  station 
2  (and  possibly  station  2*  for  by-pass  inlets).  If  a  loca¬ 
tion  beyond  station  2  (or  2*)  is  utilized  then  appropriate 
subscripts  should  be  used. 

rhz  mass  flow  rate  at  station  2. 

If  vitiated  air  is  used 

^2  =  mo,r  +  rhvf  +  mo,  (5.1) 

rhair  non  vitiated  air  mass  flow  rate 
rn^/  vitiator  fuel  mass  flow  rate 

wo,  oxygen  make-up  mass  flow  rate  (vitiator 
oxidizer  and  replenishment) 

P2  or  P<2  static  or  total  pressure  at  station  2 

7t2  total  temperature  at  station  2 

For  by-pass  inlets  with  SFRJ  (Figure  3.2)  rh2*,  P2-  or  p!2* 
and  T(2»  are  also  needed. 

rh/  fuel  or  propellant  mass  flow  rate 

Pi  static  pressure  at  station  4 

Pi  base  pressure 

Pamh  local  ambient  pressure 
Fi,c  load  cell  force 
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If  heat  losses  or  thermal  induced  area  changes  are  to  be 
taken  into  account,  additional  measurements  are  required. 
(For  example,  if  a  water  cooled  device  is  applied,  mz/jOi 
Tn:,o,in  and  Tn^o.out  are  required.) 


5.2  Typical  Methods  for  Measur¬ 
ing  Parameters 


It  is  assumed  that  mass  flow  rate,  pressure,  temperature 
and  force  are  measured  using  conventional  techniques.  It 
is  further  assumed  that  pressures  and  temperatures  should 
represent  appropriate  averaged  values  across  the  flow  area. 
Special  ramjet  measurements  are  delineated  below. 


5.2.1  Nozzle  Discharge  Coefficient 


SFRJ; 


with 


m/  = 


h' 


(5.3) 


mj  i—  TUfj  total  burnt  mass  of  fuel 

tb  burn  time 


•  DR  with  choked  gas  generator,  constant  injection 
throat  area  and  for  a  constant  c*: 

.  Jo 


with 

Pt^g  gas  generator  total  pressure 

(Practically,  the  gais  generator  static  pressure, 
Pg,  is  measured  and  used  in  the  instantaneous 
fuel  mass  flow  rate  calculation  Eq.  5.4.) 


The  discharge  coefficient  is  traditionally  a  streamline  cur¬ 
vature  correction  to  yield  the  effective  flow  area  for  one¬ 
dimensional  isentropic  flow. 

For  sonic  flow  the  following  equation  can  be  used  to  esti¬ 
mate  cd3  experimentally  (using  pre-  and/or  post-test  air 
flow): 

(:^)  (5.2) 

with 

Rs  =  R/Ms 


•  DR  with  unchoked  gas  generator  or,  choked  gas  gen¬ 
erator  with  variable  injection  throat  area: 

=  ejSi{t)rt(t)  (5.5) 

with 

■  Qf  propellant  density 

Si  propellant  burning  surface  area 

Ti  '  burning  rate 

The  burning  rate  for  solid  propellants  can  be  mea¬ 
sured  directly  or  calculated  with  a  burning  law,  gen¬ 
erally  expressed  by: 


TZ  universal  gas  constant 

Mh  molecular  weight  at  station  5 

75  =  7/5  or  =  7,5  or  =  7p,4-5  (Appendix  B) 

However,  it  is  generally  very  difficult  to  obtain  an  accurate 
value  of  C£)5  because  of  Inaccuracies  in  the  measurements 
of  mass  flow  rate,  pressure  and  temperature  and  the  effect 
of  molecular  weight  and  specific  heat  ratio. 

If  accurate  values  cannot  be  determined  experimentally 
then  CD  5  may  be  estimated  from  reference  [33]. 


5.2.2  Fuel  Mass  Flow  Rates  for  SFRJ  arid 
DR 


rt  =  apj 

(5.6) 

with 

a 

a  parameter  which  depends  primarily  upon 
propellant  temperature 

n 

pressure  exponent 

The  parameters  a  and  n  can  be  a  function  of  pressure 
and  initial  propellant  temperature. 

For  solid  propellant  configurations  with  varying  burn¬ 
ing  area,  there  is  a  relation  between  burning  grain 
area  and  burnt  thickness  such  that  Si(t)  = 
where  this  Izist  value  is  obtained  by 

Ei(t}=  I  n{t)di  (5.7) 

Jo 


In  the  cases  of  the  solid  fuel  ramjet  and  ducted  rocket, 
the  mass  flow  rate  is  not  measured  directly,  but  can  be 
approximated  by  the  methods  given  below. 


The  method  also  requires  knowledge  of  the  relation 
between  the  burning  rate  and  the  gas  generator  pres¬ 
sure. 
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5.3  Useful  Data  not  Essential  for 
Performance  Calculations 


For  a  better  understanding  of  ramjet  engine  behaviour, 
other  useful  data,  although  not  essential,  are  frequently 
taken.  For  instance: 


•  Combustor  surface  temperature  (e  g  thermocouple, 
pyrometer,  IR  thermography). 

•  Direct  instantaneous  measurements  of  DR  and  SFRJ 
fuel  regression  rate  (e.g  ultrasonic  or  X  ray  methods) 

•  Exhaust  plume  signature  (e.g  temperature,  IR,  parti¬ 
cle  size). 

•  Local  flowfield  (e  g  velocity,  temperature,  distortion, 
turbulence).  . 


5.4  Pressure  Oscillations  and 
Combustion  Instabilities 

The  analysis  of  pressure  oscillations  and  combustion  insta¬ 
bilities  is  a  difficult  problem  and  requires  a  specific  detailed 
description.  Nevertheless,  it  is  generally  assumed  that  un¬ 
steady  pressure  measurements  are  made  (e.g  piezoelectric 
transducers)  at  locations  on  the  combustor  wall  (water 
cooled  transducers)  and/or  the  inlet  ducts. 

From  the  signal-time  histories  that  are  obtained,  fre¬ 
quency,  amplitude  and  phase  data  can  be  determined. 


25 


6  Experimental  Performance  Evaluation 


The  performance  evaluation  of  ramjet  motors  is  accom- 
plished  with  the  help  of  characteristic  parameters.  Usu¬ 
ally  these  parameters  are  determined  by  a  combination 
of  measurement  and  analytical  calculation.  Typical  ram- 
combustor  performance  parameters  are: 

•  Combustion  efficiency 

•  Total  pressure  loss 

•  Nozzle  expansion  efficiency 

•  Expulsion  efficiency  of  the  fuel  tank  or  gas  generator 

In  addition,  other,  parameters  characterizing  the  opera¬ 
tional  function  or  limits  of  the  engine  may  be  used.  They 
include  ignition  limits,  blow-off  limits  and  combustion  star 
bility.  The  latter  parameters  are  not  addressed  in  this 
report, 


6.1  Assumptions  and  Procedures 

Required  measured  parameters  are  listed  in  Chapter  5, 
These  parameters  are  used  as  the  data  base  for  the  cal¬ 
culations  which  can  be  performed  using  the  equations  in 
this  chapter. 

Furthermore,  for  the  calculation  procedure  several  as¬ 
sumptions  are  necessary: 

•  One  dimensional  flow. 

•  Mass  flow  at  station  4  is  the  same  as  at  station  5: 
rTi4  =  ms.  That  means  that  the  complete  mass  flow 
of  the  combustor  is  expanded  through  the  nozzle.  In 
the  following  sections  either  rh4  or  rhs  are  used,  de¬ 
pending  on  which  is  most  appropriate  to  the  context. 

•  Total  pressure  losses  between  stations  4  and  5  are 
neglected:  p(4  =  p,5 

•  Heat  losses  between  stations  4  and  5  are  neglected: 

—  T(5 

•  Between  stations  4  and  5  isentropic,  equilibrium  flow 
is  assumed. 

•  Normally,  the  Mach  number  at  station  5  is  Ms  =  1. 
The  nozzle  throat  is  choked. 


The  evaluation  process  necessitates  knowledge  of  the  com¬ 
position  of  the  combustion  gases.  In  principle,  a  chemical 
analysis  of  the  combustion  products  is  needed.  Because 
of  the  high  effort  that  this  would  entail,  and  in  order  not 
to  introduce  additional  sources  of  error,  this  procedure  is 
avoided  in  most  cases.  Instead,  it  is  assumed  that  the  com¬ 
position  of  the  combustion  gases  is  the  same  as  in  the  case 
of  chemical  equilibrium.  Concerning  the  physical  model, 
upon  which  the  evaluation  process  is  based,  this  means 
that  the  energy  losses  of  incomplete  reactions  are  substi¬ 
tuted  by  simple  heat  losses.  Only  in  the  case  of  highly 
incomplete  reaction  within  the  combustion  chamber  will 
this  assumption  lead  to  noticeable  errors  of  evaluation:'  the 
better  the  efficiency,  the  truer  the  assumption. 

Referring  to  the  calculation  procedures,  there  are  two  ap¬ 
proaches  possible. 

One  is  to  relate  stagnation  to  static  properties  (Appendix 
B).  Any  use  of  y-values  is  to  some  extent  inaccurate, 
but  the  inaccuracies  are  small.  It  is  recommended  in  this 
procedure  to  use  the  between  the  chamber  stagnation 
and  static  throat  conditions 

The  more  correct  procedure  is  to  extract  from  the  aero- 
thermochemical  equilibrium  calculation  the  direct  rela¬ 
tionship  between  the  several  parameters.  This  does  not  re¬ 
quire  the  determination  of  y-values.  The  normal  approach 
is  to  apply  the  theoretical  values  of  adiabatic  combustion. 
A  refinement  by  using  a  plausible  average  combustion  ef¬ 
ficiency  in  the  theoretical  calculation,  may  be  reasonable 
in  some  cases.  The  effort  of  a  true  iteration  is  usually  not 
worthwhile,. 


6.2  Combustion  Efficiency 

In  general,  an  efficiency  definition  compares  a  measured 
performance  value  with  a  theoretically  evaluated  one.  The 
resulting  value  gives  information  about  the  quality  of  the 
examined  process.  Hence,  the  general  formula  is: 

experimentally  determined  value 
theoretically  determined  value 

where  the  numerator  and  denominator  consist  of  characte¬ 
ristic  performance  parameters.  Due  to  practical  or  tradi¬ 
tional  reasons,  several  definitions  of  combustion  efficiency 
are  in  use.  Presented  herein  are  combustion  efficiencies 
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based  on  characteristic  velocity  ('?£-),  vacuum  specific  im¬ 
pulse  (ijJoc),  temperature  rise  (iJat)  and  equivalence  ratio 

Incomplete  combustion  of  ramjet  fuel  (with  some  fuel  un¬ 
burned)  in  fuel-rich  situations  results  in  a  nonequilibrium 
combustion  temperature  that  is  higher  than  the  theoret¬ 
ical  equilibrium  combustion  temperature.  Thus,  combus¬ 
tion  efficiencies  greater  than  unity  may  be  calculated  when 

There  are  two  basic  ways  to  determine  combustion  effi¬ 
ciency,  either  by  measured  static  chamber  pressure  or  by 
measured  thrust.  Both  ways  can  be  executed  utilizing 
isentropic  exponents  (Appendix  B)  .or,  as  recommended 
in  6.1,  by  utilizing  parameters  derived  from  the  aerother- 
mochemical  equilibrium  code,  not  using  7-values.  The  re¬ 
sulting  four  methods  are  outlined  below. 

Derivations  of  the  equations  utilized  in  this  chapter  are 
given  in  Appendix  C. 


6.2.1  Efficiency  Based  on  Characteristic 
Velocity 


Then 


Tp,  r 

Pm  =  Pvxp  (^1  +  (6.4) 

2.  Based  on  measured  static  pressure  without  using  y: 
Ph  can  be  determined  using  the  following  expression 


where  {pm/p4)m  is  obtained  from  an  aerothermo- 
chemical  equilibrium  code  for  a  specified  area  ratio 
Ai/(AsCDs)- 


3.  Based  on  measured  thrust  utilizing  7.’ 

In  the  case  of  thrust  measurement  with  a  convergent 
nozzle  the  total  pressure,  at  station  4  (assumed  =  p(s) 
can  be  determined  from  the  following  equation; 


_  -^5  +  PamhAj  f  7p,i  +  1 

(^  +  yp.i<^Ds)As  \  2  ) 


(66) 


Fs  is  the  load  cell  thrust,  corrected  for  base  pres¬ 
sure  force  and  preloads  on  the  thrust  stand 


The  c*  efficiency  compares  the  characteristic  velocities  of 
the  exhaust  jet  derived  from  experiment  and  from  theory. 


'-trp 

(6,1) 

given  by; 

.*  _  PM^sCDS 

(6.2) 

where 

4.  Based  on  measured  thru.st  without  utilizing  7: 
Pm  {=  Pis)  is  given  by: 


( ^5  d-PomtAs 
AsCdS 


(6.7) 


6.2.2  Efficiency  Based  on  Vacuum  Spe¬ 
cific  Impulse 


c*(,  is  calculated  with  the  aerothermochemical  equilib¬ 
rium  code 

Cos  is  the  discharge  coefficient  at  station  5 
rh4  is  given  by: 

ITI4  =  ma.r  -f  rh^it  +  m/ 


The  methods  to  determine  the  total  pressure  pn  are  as 
follows; 


describes  the  relation  of  an  experimental  vacuum  spe¬ 
cific  impulse  to  the  theoretical  value: 


*vdc,<gp 

7* 

ac^ih 


(6.8) 


is  defined  as  the  thrust  per  unit  mass  flow  of  a  con¬ 
vergent  nozzle  discharging  into  a  vacuum.  This  value  is 
identical  with  the  stream  thrust  per  unit  mass  flow  in  the 
sonic  throat: 


Fss 

rhs 


(69) 


1.  Based  on  measured  static  pressure  utilizing  7; 
In  this  case,  M4  must  first  be  calculated. 


CDS  As 
A4 


( 


2 

Tp,^  +  1 


-I- 


7p,5  -  1 

7p,.  +  1 


tp,  I  +1 


(6.3) 


**flc  1/1  be  obtained  using  an  aerothermochemical  equi¬ 
librium  code  and  the  following  equation: 


.  Vae.lfi 


rilscs  +  P5A5 


(6.10) 


i*ac,crp  6e  derived  from  one  of  the  following  tech¬ 
niques: 
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1.  Based  on  measured  static  pressure  using  7: 


1,  pressure  measurement  at  station  4  using  7; 


Up.  +  lj  (1  +  7p,.CD5) 

(6.11) 

where  pm  is  obtained  from  Eqs.  6.3  and  6.4. 

2.  Based  on  measured  pressure  without  using  7; 

■  (6.12) 

where  c*^j,  is  calculated  with  Eqs.  6.2  and  6.5. 

3.  Based  on  measured  thrust  using  7. 

In  this  particular  case  method  3  does  not  exist  be¬ 
cause  7  does  not  appear  in  the  resulting  expression. 


(6.16) 


or; 


lt4,6Xp 


‘’tacicrp 


2(7p,J  H"  1)  ^4, exp 


(6.17) 


For  practical  reasons,  is  replaced  by  Ri^th.  ob¬ 

tained  from  an  aerothermochemical  equilibrium  code. 
c*^p  is  obtained  from  Eqs.  6. 2-6. 4  and  from 

Eq  6.11. 


2.  pressure  measurement  at  station  4  without  using  7: 


4.  Based  on  measured  thrust  without  using  7; 

J,  ^5+Pamk^5 


\jat,txp 


m5 


(6.13) 


or, 


There  exists  another  method  using  measured  static  pres¬ 
sure  in  the  sonic  throat; 


Cc,«p=,(Tp,,cd5  +  1)^  (6.14) 

This  approach  is  not  often  used  due  to  the  difficulties  of 
measuring  ps, 

The  efficiencies  and  rjf  are  directly  equivalent. 


The  same  comment  on  /?4,eip  as  given  above  applies 
here  also.  c*j.p  is  obtained  from  Eqs.  6.2  and  6.5  and 
*J<ic,«p  from  Eq.  6.12. 

3.  thrust  measurement  using  7; 

Use  Eq.  6.16  together  with  Eq.  6  2  and  Eq.  6.6  to 
obtain  jp 


Note  that  the  use  of  either  efficiency  depends  on  the  type 
of  measurement.  If  the  test  arrangement  has  a  convergent 
nozzle  and  thrust  is  measured,  the  evaluation  based  on 
is  more  straight  forward.  If  only  pressure  is  measured,  it 
is  more  appropriate  to  use  c*. 


6.2.3  Efficiency  Based  on  Temperature 
Rise 


The  value  shows  the  experimental  stagnation  tem¬ 
perature  rise  in  comparison  to  the  theoretical  temperature 
rise  in  the  combustion  chamber: 


Tt4,exp  —  7)2 

~  Ti2 


(6.15) 


4,  thrust  measurement  without  using  7;  ■ 

Use  Eq.  6.19  together  with  Eq.  6.13  to  obtain  Tm.sip- 


6.2.4  Efficiency  Based  on  Equivalence  Ra¬ 
tio 


The  efficiency  of  a  combustion  process  can  be  character¬ 
ized  by  comparison  of  the  experimentally  injected  equiv¬ 
alence  ratio  against  the  theoretical  value  which  is 
necessary  to  gain  the  experimentally  determined  perfor¬ 
mance  (c*,  J. 


Vi  — 


<f>i 


(6.20) 


where 

Tt^^th  can  be  obtained  from  an  aerothermochemical  equi¬ 
librium  code 

In  principle  Tti^cxp  can  be  measured  directly  (total  tem¬ 
perature  probes,  calorimeters  etc.)  but  good  results  are 
difficult  to  achieve.  In  most  cases  the  total  temperature 
at  station  4  is  determined  directly  from  the  experimen¬ 
tally  obtained  characteristic  velocity  and/or  vacuum  spe¬ 
cific  impulse,  using  one  of  the  following  methods: 


with; 

^  (6.21) 

For  example,  an  aerothermochemical' equilibrium  code 
is  used  to  generate  the  theoretical  relationship  between  c* 
and  4)  (Fig.  6.1).  The  figure  is  entered  with  to  obtain 
<j)b'  Normally,  this  technique  is  used  for  experimentally 
injected  equivalence  ratios  <  1.  c*  varies  insignificantly 
with  small  changes  in  pressure,  Therefore,  pn  obtained 
using  any  of  the  above  mentioned  methods  is  suitable  for 
input  to  the  code. 
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c*  as  function  of  <j> 


this,  they  become  bulky.  Therefore,'  there  is  a  risk  that 
the  test  results  may  be  affected  by  blockage  of  the  cham¬ 
ber  cross-section.  If  metallized  propellants  are  used  the 
severe  problem  arises  to  protect  the  orifices  from  being 
clogged  by  the  combustion  products.  Therefore,  in  most 
cases  the  total  pressures  are  derived  from  the  measure¬ 
ments  of  static  wall  pressures  or  thrust. 

In  the  following  sections  the  equations  for  determination 
of  the  total  pressures  are  presented. 


6.3.1. 1  Evaluation  of  P( 2 
Pt2  may  be  determined  using 

Pt2=P2[i  +  ^^^Miy"'  (6.24) 


Figure  6.1:  Principle  of  determination  of  from  c*  =  ^2  is  found  as  the  subsonic  solution  of  the  following  rela- 

f{4>) 


6.3  Additional  Performance  Pa¬ 
rameters 


P2A2 


72 


(6.25) 


6.3.1  Pressure  Losses 


Usually,  pressure  losses  are  characterized  by  the 
total  pressures: 

PtA 
Pi  2 

or  by  the  relative  difference  of  total  pressures: 


ratio  of 
(6.22) 


Referring  to  section  6.1  it  is  best  to  use  the  process  7  be¬ 
tween  stagnation  and  static  conditions  at  station  2.  Low 
Mach  number  and  moderate  temperatures  at  station  2  per¬ 
mit  the  process  7  to  be  replaced  by  7  values  obtained  from 
air  tables  at  the  measured  temperature  Tn- 


6.3. 1.2  Evaluation  of  pt 4 


P<2  —  Pm 
Pt2 


(6.23) 


Pressure  losses  occur  at  different  places  in  the  ramjet  en¬ 
gine  The  above  mentioned  pressure  ratios  give  only  global 
values.  The  following  pressure  losses  may  be  found  in  a 
ramjet  engine: 


The  total  pressure  at  station  4  can  be  derived  from  mea¬ 
surements  either  of  static  pressure  p4  or  of  the  thrust. 
In  addition,  tne  evaluation  of  pt4  can  be  done  with  or 
without  utilising  isentropic  exponents,  using  the  formulas 
presented  in  Section  6  2.1. 


6.3.2  Expulsion  Efficiency 


•  aerodynamic  stagnation  pressure  losses  caused  by 

o  sudden  expansion  (dump)  from  station  2  up  to 
station  4  (Carnot  diffusor) 

o  wall  friction 
o  flow  turning 

o  fuel  injector  and  flameholder  drag 

•  stagnation  pressure  losses  caused  by  combustion 


The  expulsion  efficiency  characterizes  the  completeness  of 
the  fuel/propellant  utilization  and  has  the  following  defi¬ 
nition: 

_  _  _  ,  Am/ 

f/ex  —  — 

mj^i 

TTif^i  is  the  total  stored  propellant.  Am/  is  the  residue  of 
the  fuel/propellant  in  the  tank  system,  gas  generator  or 
combustor  (in  the  case  of  the  solid  fuel  ramjet)  after  use. 


(6.26) 


Direct  measurement  of  the  total  pressures  with  rakes 
of  pitot-tubes  mostly  leads  to  complicated  test  arrange¬ 
ments.  The  pitot-tubes  must  be  cooled  and,  because  of 


The  liquid  fuel  ramjet  has  a  rather  high  expulsion  effi¬ 
ciency  of  about  0.98.  Fuel  residues  are  due  to  wetting  of 
the  bladder,  remainder  in  pipes  and  pumps,  etc. 
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High  expulsion  efficiencies  are  not  easily  obtained  for  solid  significance,  and  have  to  be  calculated  or  at  least  esti- 
propellant  ducted  rocket  gas  generators.  However,  with  mated  by  appropriate  methods  before  the  expansion  effi- 
a  proper  design  of  the  propellant  and  gas  generator  an  ciency  is  applied, 
expulsion  efficiency  of  above  0.96  is  possible. 

The  wall  heat  transfer  is  of  noticeable  influence  only  in 
The  expulsion  problems  are  even  more  critical  in  the  case  the  case  of  heat-sink  or  active  cooling  and  then  has  to  be 
of  the  solid  fuel  ramjet,  primarily  due  to  non-uniform  considered, 
burning  of  the  fuel  grain.  Generally,  the  expulsion  effi¬ 
ciency  will  remain  below  0.95.  The  two-phase  flow  effects  which  significantly  reduce  the 

performance  of  solid  propellant  rocket  motors  are  of  mi¬ 
nor  importance  in  the  case  of  ramjets.  Due  to  the  high 

6.3.3  Nozzle  expansion  efficiency  f solid  mass  fraction  of  ramjet  ex¬ 

haust  IS  always  small,  even  if  metals  are  applied  as  a  fuel 
compound.  Nevertheless,  the  negative  influence  has  to  be 
In  section  6.1  the  evaluation  of  ramcombustor  performance  considered  and  it  is  recommended  to  reduce  the  expan- 
was  based  on  the  assumption  that  the  flow  within  the  con-  sion  efficiency  by  about  one  percentage  point  if  a  highly 
vergent  part  of  the  nozzle  was  isentropic.  This  idealization  metallized  propellant  is  used, 
is  admissible  since  the  losses  within  a  well-shaped  conver¬ 
gent  nozzle  are  small.  Moreover,  the  idealization  of  the 
flow  process  in  the  subsonic  part  of  the  nozzle  does  not 
suppress  these  small  losses,  but  only  shifts  them  to  the 
performance  balance  of  the  combustor. 

The  losses  within  the  supersonic  part  of  the  thrust  nozzle 
can  be  higher.  Losses  are  caused  by: 


1.  wall,  friction 

2.  divergence  and  local  shocks 

3.  heat  transfer  to  the  wall 

4.  two-phase  flow  effects 

5.  incomplete  recombination  (or  dissociation) 


The  first  two  types  of  losses  are  approximately  propor¬ 
tional  to  the  exit  stream  thrust  and  are  usually  considered 
by  the  expansion  efficiency  being  defined  as  follows 
(assuming  no  base  drag); 


(P6^6  +  m6C6)ejp 
(pe-de  +  rTi6C6)t)( 

Te  T  Pamb-^S 

[P6d6(l  +  T6jW|)]tA 


(6.27) 


where  Fg  equals  the  load  cell  thrust  corrected  for  any 
preloads  on  the  thrust  stand  and  base  pressure  forces. 


In  the  regime  of  moderate  supersonic  speeds  up  to  about 
Mach  4,  the  aerodynamic  types  of  losses  (”!”  and  ”2”)  are 
dominant,  and  the  expansion  efficiency  is  a  useful  tool  to 
determine  the  effective  thrust  from  the  theoretical  values 
calculated  with  the  aerothermochemical  equilibrium  code. 
With  well  shaped  thrust  nozzles,  expansion  efficiencies  of 
about  0.98  are  possible. 


In  the  hypersonic  flight  regime,  the  nozzle  performance 
losses  due  to  incomplete  recombination  have  increeising 
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7  Sample  Calculations 


Sample  calculations  are  provided  to  assist  the  reader  with 
applications  of  methodologies  established  in  this  report. 
Ail  aerothermochemical  equilibrium  calculations  were  per¬ 
formed  with  the  help  of  the  NASA  CET89  code  [9)  accord¬ 
ing  to  the  assumptions  as  outlined  in  Chapter  6.  Results 
from  four  cases  are  presented.  Case  1  is  for  a  LFRJ  per¬ 
formance  without  vitiation  heating  of  the  test  air  medium 
(ideal  air).  For  this  case,  calculation  results  are  provided 
for  combustion  efficiency  based  on  characteristic  velocity, 
vacuum  specific  impulse,  temperature  rise  and  equivalence 
ratio.  Each  efficiency  parameter  is  computed  both  from 
measured  static  combustor  pressure  and  measured  sonic 
thrust.  For  each  type  of  measurement  calculations  were 
made  with  and  without  7p,,.  Therefore,  results  of  the  four 
computational  methods  are  presented.  Case  2  is  a  sam¬ 
ple  calculation  for  LFRJ  performance  when  a  vitiated  air 
heater  is  employed.  Cases  3  and  4  are  sample  calculations 
only  for  mass  flow  rates  for  the  DR  and  SFRJ,  respec¬ 
tively,  because  the  remaining  calculations  are  identical  to 
those  of  the  LFRJ.  Tables  7.1  and  7.2  present  a  summary 
of  the  equations  necessary  for  the  calculation  of  the  per¬ 
formance  parameters  using  the  four  methods  presented  in 
Chapter  6.  :■ 


7.1. 1.1  General  test  information 

•  Fuel:  CiQ/f20  kerosene 

•  Air;  Afi56.204i.96-^r’o. 934^0.0314 

7.1. 1.2  Geometric  data 

•  Flow  areas 

A2  =  0.010325 
A4  =  0  022698 

•  Exhaust  nozzle  (Fig.  5.1  and  C.l) 

A5  =  0.012668 

cd5  =  0.996 
Ai  =  0.004304 


Combustor  and  nozzle  heat  losses  are  neglected  in  the  the¬ 
oretical  calculations  of  the  following  examples.  For  all  cal¬ 
culations  the  actual  input  and  output  files  for  the  NASA 
CET89  are  given  in  Appendix  D. 


7.1  LFRJ  Performance  with  Ideal 
Air  (Case  1) 


7, 1.1. 3  Measured  test  data 
•  Mass  flows 


m„/  =  0.0  kg/s 

moj  =  0.0  kg/s 
rhair  =  6.692  kg/s 
mj  =  0.311  kg/s 


•  Thrust  stand  forces 


This  case  is  for  a  kerosene-fueled  LFRJ  tested  in  an  ideal 
air  test  medium  (no  vitiator).  Detailed  procedures  are 
presented  only  for  the  method  “measured  static  combustor 
pressure  without  using  7" .  However,  numerical  results  for 
all  four  methods  are  summarized  in  Table  7.3. 


Fic  =  13400  N 
Ftart  =  5000  A" 


•  Pressures  and  temperatures 


7.1.1  Calculation  Procedure 

The  following  describes  a  methodology  that  can  be  fol¬ 
lowed  to  calculate  the  efficiency  parameters  and  is  orga¬ 
nized  in  the  order  the  information  is  required,  and  the 
order  of  Chapter  3,  Methods  for  Reporting  Test  Results. 


Pumh 

101300  Pa 

P2 

= 

650200  Pa 

P4 

= 

568800  Pa 

P5 

=1 

341203  Pa 

Pi 

= 

78065  Pa 

Tt2 

= 

606  R 

Tj 

= 

298.15  A' 
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Method 

1 

2  . 

Calculation 

Procedure 

P4 

with  7 

P4 

without  7 

_ * _ 

Parameter 

Equation 

Source 

Equation 

Source 

Mi 

Eq.  6.25 

Method  1 

Pt2 

P«2  =  P2  (1  -I- 

Eq.  6.24 

Method  1 

Tm,  Pi4  1**  It. 

Pl4  —  P4,ezpf  ^ 

‘ 

Method  1 

■  ln(Pn/Ps.iK) 

■  ln(pii/p^.ih)-in(Ti,/Ts.i>,) 

Eq.  B.IO 

not  required 

Mi 

Eq.  6.3 

not  required 

P«4 

Eq.  6.4 

P“'= 

Eq,  6.5 

rerun  aerothermochemistry  code  using  pti|',  yields  new  theoretical  values 

^4,ezp  —  ■^4,th 

Method  1 

'  c* 

1  '-tip 

fjk  Pf4AaCDa 

‘'etp  -  m. 

Eq.  6.2 

Method  1 

Tti,izp 

T  r  f  ^ 

Eq.  6.16 

JM,«p  -  t  £*j  )th 

Eq,  6.18 

recalculate  7p,j,  M^  and  ptn  —  iterate  to  desired  accuracy 

1  ^  (^)  (1  +  Tp,.CI55) 

Eq.  6.11 

•voc.erp  “  ^txp  \  c’ 

Eq.  6,12 

/)i,  Pm  I  c*j.p  and  code 

Fig.  7.1 

Method  1 

(f/a)stcich  = 

Method  r 

1  (^/ /^*  1  r  )er|i 

Eq.  6.21 

Method  1 

r}c’ 

Eq.  6.1 

Method  1 

0 

d 

_  _  *wac.erp 

r?i*  =  7^-' — 

‘vc.lk 

Eq.  6.8 

Method  1 

r}£iT 

— T|a 

Eq.  6.15 

Method  1 

Eq.  6  20 

Method  1 

Pti/pa 

£m  ■ 

Pia 

Eq  6.22 

Method  1 

{pt2  -Pm)/pi2 

Pi3-Pm 

Pl2 _ 

Eq  6.23 

Method  1 

All  required  theoretical  values  are  obtained  from  the  NASA  CET89  code; 

Pti/Pi,  7?4,<A,  72.  <I>I  Tti^th,  Ps.ts.  Kac 

Table  7.1:  Summary  of  equations  for  performance  calculations 


7. 1.1. 4  Preliminary  Calculations 


The  following  preliminary  calculations  must  be  made  to 
provide  inputs  for  the  efficiency  parameter  equations  in 
Chapter  6: 


•  Fuel  enthalpy  Enthalpy  of  fuel 


Since  the  aerothermochemical  equilibrium  code  re¬ 
quires  enthalpy  units  in  cal /mole: 


>  1000  cal 

hi  =  -2016.6-;— X X 

^  kg  4.184  W 


1  kg 
1000  g 

=  -67607 


X  140.27 
cal 


mole 


mole 


•  Nozzle  stream  thrust  (Eqs.  C.25-C.26) 


^LC  ^tare 

.13400  -  5000 
8400  N 


book  value  at  measured  Tj 
—2016.6  kJ/kg 
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Method 

3 

4 

Calculation 

Procedure 

Fs 

with  7 

Fs 

without  7 

Parameter 

Equation 

Source 

Equation 

Source 

M2 

Method  1 

Method  1 

Pt2 

Method  1 

Method  1 

7)4  ■Pl4  1’*  ft. 

Method  1 

Method  1 

7fi,s 

Method  1 

not  required 

M4 

Method  1 

not  required 

Pt4  ■ 

’9.* 

_  r5+p.„i.4j  /v,,.+i\ 

Eq.  6.6 

(  c’  \  Fj+P.mivtj 

Eq.  6.7 

(1+7p.,CD5)j45  V  2  / 

rerun  aerothermochemistry  code  using  p^;  yields  new  theoretical  values 

^4, exp 

Method  1 

Method  1 

Method  1 

Method  1 

^(4, exp 

Method  1 

Method  2 

recalculate  7p  j,  Af^  and  p^  —  iterate  to  desired  accuracy 

t* 

‘vac.exp 

no  procedure 

2*  _ 

Sac, exp  rht. 

Eq.  6.13 

4>b 

Method  1 

Method  1 

Method  1 

Method  1 

Method  1 

Method  1 

Method  1 

Method  1 

’ll* 

'•voc 

Method  1 

Method  1 

TJaT 

Method  1 

Method  1 

Method  1 

Method  1 

1 

Pt4lPt2  ■ 

Method  1 

Method  1 

(P(2  —  Pl4)/P«2 

Method  1 

Method  1 

All  required  theoretical  values  are  obtained  from  the  NASA  CET89  code: 

P(4/P4,  ,thi  B21  72i  C*f,, 


Table  7  2'.  Summary  of  equations  for  performance  calculations  (cont’d) 


7^5  —  f^mcoj  Aj(Pj  Pamb) 

=  8400-0.004304(78065-  101300) 

=  8500' JV 

•  Station  2  conditions 

With  tables  or  an  aerothermochemical  equilibrium 
code  (Fig.  D.l),  determine  the  following  at  station 
2 


o  inputs:  P12  (assumed  equal  to  measured  p2),  7)2, 
air  composition 

o  output; 

M2  =  28.965  kff/kmole 
72  =  1  3750 

^2  =  ht2  =  hair 


■=  310.9  JcJ/Jcff 

converting  to  cal(molt: 

I  1000  cal 

hair  =  310.9t—  X  X 


kg 


1  kg 
1000  g 

=  2152.4 


4.184  kJ 
9 


X  28.965 
cal 


mole 


mole 


^2  = 


M2 
8314.51 
28,965 
=  287.05 


kg  K 
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•  Chamber  total  pressure  (pn) 

The  process  for  determining  pt4  requires  several  se¬ 
quential  runs  of  the  aerothermochemical  equilibrium 
code.  Method  2  (Table  7.1)  is  being  used  for  the 
sample  calculation.  However,  for  clarity,  the  required 
processes  to  obtain  pa  for  all  four  methods  are  shown. 

For  all  methods,  run  the  aerothermochemical  equilib¬ 
rium  code  using  a  measurement  or  estimate  of  (in 
this  example  p4  =  568800  Pa)  as  a  first  approxima¬ 
tion  to  generate  P(4- 


o  input  (Fig.  D.2):  p4,  Aa/IAsCds),  ma,r/m^ 


A4 

Ascps 

iTlflir ■ 
Thj 

o  output  (Fig.  D.2): 


P5,ih 


^4,1(1 


4 


0.022698 
0.012668  X  0.996 
1.7990 
SUBARCl) 

6.692 

0.311 

21.5177 

mx(i) 


314560  Pa 

1.0793 

2065.80  K 
1834.46  K 
1170  m/s 
1462.1  Nsjkg 
0.6877 


o  obtain  first  estimate  for  the  chamber-to-throat 
process  7  (tp,.): 
with  p£4  p4  and  3)4  ^  74 


Pm  _ 

568800 

P5,«h 

314560 

= 

1.8082 

Tt4 

2065.80 

1834.46 

= 

1,1261 

From  Eq,  B.IO; 

In  (pi4/p5,tft) 

“  ln{p£4/P5,M)-ln(r,4/T 

In  (1.8082) 

In  (1.8082) 

-In  (1.1261) 

=  1.2508 

Method  1;  using  mezisured  p4  and  7 


o  determine  the  combustor  Mach  number  using 
Eq.  6.3 


CD5A5 


4- 


VTp.j  +1  '  7p,»  +  1 
0.996  X  0.012668 


M4 


^  X’p.-'') 


0,022698 


(r 


2508  4-1 


4- 


1.2508-  1 
1.2508  4-1 
0.3504 


1  .aaofl-f  2 

1) 


o  compute  pt4  using  Eq.  6.4; 


P,4  =  P4  [^1  4- 

=  568800  X 

,  1.2508-1 

=  613838  Pa 


— 1 


0.3504- 


i.23oa_ 


o  rerun  the  aerothermochemical  equilibrium  code 
using  pi4  =  613838  Pa  (all  other  input  data  un¬ 
changed) 

output  (Fig.  D.3): 


PS.ih 

=  339460  Pa 

Tt4,ih 

=  2066.02  A" 

=  1834.49  K 

C*th 

=  1170m/s 

■J* 

^vac^th 

=  U&2.\Nsfkg 

Mi 

=  28.909  ij/fcmofe 

o  recalculate  7^,,,  using  Eq.  B.IO: 


7p,,  =  1.2510 


and  p(4  using  Eq.  6.3  and  6.4 


P!4  =  613845  Pa 


o  perform  additional  iterations  for  increased  accu¬ 
racy,  only  if  the  successive  values  in  pt4  differ 
significantly 

Method  2:  using  measured  P4  without  7 


o  calculate  combustor  total  pressure  using  Eq.  6.5 

( 

Pi'l  —  I  )  P'l.erp 
V  P'1  /  (4 

=  1,0793  x  568800 
=  613906  Pa 


o  rerun  aerothermochemical  equilibrium  code  at 
Pf4  =  613906F’a  (all  other  input  data  un¬ 
changed) 
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output  (Fig.  D.4); 

VtilPi  =  1.0793 

2^(4, <fc  —  2066.02  K 
Ts„th  .=  ISM  49  K 
c*)j  =  1170  m/s 

Cch  =  1462  1  m/s 
Mi  =  28.999  kg/ kmole 

Method  3:  using  measured  thrust  and  7 

o  calculate  chamber  total  pressure  using  Eq.  6.6 


_  F5  +,  Pamirs  / * 

(1  +  7p,3Cd5}A5  \  2  J 

8500+  101300  x  0  012668 
(1  +  1.2508  X  0.996)0.012668  ^ 

n.2508+  1\  ‘  «“•'> 

=  619848  Pa 


output  (Fig.  D.6); 

Ti4.ih  =  2066.05  1< 

Cff^  =:  1170  rn/s  , 

Cc.is  =  1462.1  lVs/l:s 

A44  =  28.906  kg/kmole 

o  reiterate  if  required 


♦  Stoichiometric  fuel/air  ratio 

using  Eqs.  4.2  and  4.3  and  <j)  from  Fig.  D.2 


rh  air 

0311 

6.692 

0.04647 


0.04647 


0.6877 
=  0.06757 


o  rerun  aerothermochemical  equilibrium  code  us¬ 
ing  pt4  =  619848  Pa  (all  other  input  unchanged) 
output  (Fig.  D.5); 


Ps.tft 

=  342780  Pa 

7i4,(h 

=  2066.05  K 

=  1834.50  A' 

C*h 

=  1170  m/s 

i* 

Soc.th 

=  1462.1  A^s/zts 

Mi 

=  28. 906  kg/kmole 

•  Mass  flows 

7^2  —  ^air 

=  6.692  kg/s 
rii^  =  rhs  =  m2  +  m/ 
=  6.692+0.311 
=  7.99Zkg/s 

7. 1.1. 5  Performance  Calculation 


o  calculate  7p,,  using  Eq.  B  10 
7p,i  =  1.2510 
and  p(4  using  Eq.  6.6 


The  performance  parameters  can  now  be  calculated  by 
employing  the  equations  of  Chapter  6  and  using  the  output 
shown  in  Fig.  D.l  for  Station  2  and  Figs.  D.2  and  D.4 
for  Stations  4  and  5.  All  required  equations  are  given  in 
column  2  of  Table  7.1.  Results  are  presented  in  Table  7.3. 


P(4  =  619884  Pa 
a  reiterate  if  required 

Method  4:  using  measured  thrust  without  7 


o  calculate  chamber  total  pressure  using  Eq.  6.7 

_  /  C*  \  Fs  +  pamtAs 

\*Jac/  f(i  A5CD5 


8500+  101300  X  0.012668 
0.012668  X  0.996 
=  620477  Pa 


o  rerun  aerothermochemical  equilibrium  code  us¬ 
ing  pfi  =  620477  Pa  (all  other  input  unchanged) 


•  Combustion  Efliciency 


1.  Efficiency  based  on  characteristic  velocity,  ijc* 
Using  Eq.  6.2 


c 


★ 

exp 


PmAsCjis 

m4 

613906  X  0  012668  x  0.996 
7.003 

1106.08  m/s 


From  Eq.  6.1 


1106.08 

1170 

0.9454 


Method 

1 

2 

3 

4 

Calculation  Procedure 

P4  using  7 

P4  without  using  7 

Fs  using  7 

Fz  without  using  7 

W,  1  - 

0.351 

Pt2 

Pa 

706824 

7p,,  U'  iter. 

- 

1.2508 

- 

1.2508 

: 

Af4 

- 

0.350 

- 

0.350 

- 

P<4 

Pa 

613838 

613906 

619848 

620477 

^4,<xp 

J/kgfK 

287.6 

287.6 

287.6 

287.6 

C* 

^exj> 

.  mfs 

1106 

no6 

1117 

1118 

TtA.exp 

K 

1843 

1846 

1879 

1886 

7p,,  2'*  iter. 

- 

1.2510 

' 

1.2510 

- 

1* 

‘vac.cxp 

Ns/kg 

1383 

1382 

1397 

<!>h 

- 

0.567 

0.568 

0.587 

0.587 

{f  /  ^sioich 

- 

0.0676 

0.0676 

0.0676 

0-0676 

- 

0.688 

0688 

0.688 

0.688 

'Ic* 

- 

0.945 

0.945 

0.955 

0,955 

ni- 

- 

0.945 

ffi945  ■ 

- 

0.955 

- 

0.847 

0.850 

0.872 

0.877 

'la 

- 

0.825 

0.825 

0.853 

0.856 

pulpa 

0.868 

0.869 

0.877 

0.878 

(pt2  -  Pt‘i)/Pi2 

0.132 

0.132 

0.123 

0.122 

Table  7.3:  LFRJ  performance  with  ideal  air  (Case  1) 


2.  Efficiency  based  on  vacuum  specific  impulse, 

^•:.c 

From  Eq.  6.12 


=  1382.22  m/s, 


From  Eq.  6.8 

I* 

*tioc,eip 
1* 

*vac,tfi 

1382.22 

1462.1 
0.9454 


3.  Efficiency  based  on  temperature  rise,  t}^x 
with  and  Eq.  6.18 

o  ^ 

—  1.  j 

Ma 

8314.51 

28.909 

=  29n.&\JikglK 


Tm, 


exp 


(T^\ 

\  /  ih  •R'l.erp 

/2066.02X  287.6n 

V  (1170)2  ) 

1846.44  K 


1106.082 


~  Tt2 

1846.44-  606 
2066  -  606 
0.8496 

4.  Efficiency  based  on  equivalence  ratio,  rj^ 

A  plot  of  theoretical  characteristic  velocity 
versus  equivalence  ratio  (1^)  has  to  be  generated 
using  the  aerothermochemical  equilibrium  code 
as  shown  in  Fig.  D.7  and  7.1.  This  is  achieved 
by  using  the  input  file  Fig.  D.7  while  varying 
the  ratio  between  rhaiv  and  rh/.  The  equivalence 
ratio  ((/ij)  necessary  to  theoretically  produce  the 
measured  c*^j,  is  determined  from  this  plot. 

<i>i  =  0.5675 

Using  Eq.  6.21 

\  1  . 
mairj„p  if/a)  aioich 

0.311\  1 

6.692/  0,06757 
=  0.6878 


Using  Eq.  6.15. 

bAT  = 


0s 


X 


287.61 


Using  Eq.  6.20 
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Figure  7.1:  Determination  of  (fib  from  c*  =  f[(f>) 

0.5675 

0.6878 
=  0.8251 


•  Combustor  total  pressure  loss 


o  calculate  the  Mach  number  at  station  2 
using  Eq.  6.25: 


2  “ 


,,  ,,  1  3750-  1  ,3 

M2\/1  + - - Mi 


2 

6.692 


650200  X  0.010325 

Ms 


7712  /  R2Tt2 


P1A2  v  72 


^287. C 


.05  X  606 


3750 


=  0.3505 


o  calculate  the  total  pressure  at  station  2  using  Eq. 
6.24: 


/  72  —  1  ,,2'' 

Pt2  =  P2  (  1  +  ^ - Mi 


=  650200  X 

t 

/,  1.3750  -1 

=  706824  Pa 


0.3505‘ 


o  calculate  the  relative  difference  of  the  total  pres¬ 
sure 


Pi2  ~  Pt4 
Pi  2 


706824  -  613906 
706824 
0.1315 


7.1.2  Discussion  of  Results 

Combustion  efficiencies  for  Case  1,  LFRJ  performance 
with  ideal  air,  are  tabulated  in  Table  7.3  and  the  influ¬ 
ence  coefficients  and  uncertainty  levels  are  summarized  at 
the  end  of  this  section  in  Table  7.12. 

Considering  first  the  results  in  Table  7.3,  one  can  see  that 
the  various  performance  parameter  calculations  yield  dif¬ 
ferent  numerical  values.  For  example,  the  and 
for  Case  1  are  jb  95  %,  however,  is  w  85  %  and  7}^  is 
w  83  %.  The  main  reason  for  the  differences  between  the 
higher  values  (tJc*  and  and  the  lower  values  and 

r]^)  is  the  range  over  which  the  parameters  can  vary.  This 
demonstrates  that  when  examining  the  combustion  effi¬ 
ciency  calculated  for  a  ramjet,  one  must  first  understand 
the  efficiency  parameter  that  is  being  used. 

The  influence  coefficients  information  presented  m  Table 
7.12  is  more  important  to  the  reader  than  the  uncertainty 
levels,  which  are  presented  merely  to  demonstrate  the  un¬ 
certainty  methodology.  The  influence  coefficients  reflect 
the  sensitivity  of  performance  calculation  methods  to  in¬ 
put  parameters.  Critical  input  parameters  can  be  iden¬ 
tified  which  are  independent  of  facilities  or  measurement 
systems.  For  example,  for  the  tjat  based  on  combustor 
pressure,  the  input  parameters  with  the  most  influence  on 
uncertainty  (greatest  influence  coefficients)  are  Pa, 

Cd5  and  A5,  whereas  for  the  t)at  based  on  thrust,  rhair 
and  Fs  have  the  largest  influence  coefficients. 

The  reader  should  be  careful  not  to  generalize  conclusions 
from  the  uncertainty  levels  in  Table  7.12.  Uncertainty  lev¬ 
els  are  highly  dependent  on  each  test  facility,  measurement 
system,  calculation  methods,  installation  effects  and  envi¬ 
ronmental  conditions.  The  error  sources  utilized  herein 
to  estimate  uncertainties  are  not  necessarily  typical  val¬ 
ues  for  the  entire  operating  range.  Table  7. 12, shows  a 
difference  between  the  uncertainty  levels  of  the  first  two 
parameters  {pe*  and  and  the  second  two  (pat  and 

p,f,).  However,  a  1  %  uncertainty  in  pc*  or  Pi^^^  is  2  %  of 
the  entire  range  for  these  performance  parameters.  With 
Pat  and  a  1  %  uncertainty  is  1  %  of  the  range  of  these 
parameters. 


o  calculate  the  total  pressure  ratio 

Pi4  _  613906 

^  ~  7068^ 

=  0.8685 


Error  analysis  is  a:n  essential  engineering  tool  for  design¬ 
ing  measurement  systems,  selecting  calculation  methods, 
identifying  critical  data  validation  requirements,  ensuring 
compliance  with  test  data  requirements,  interpreting  test 
results,  and  providing  a  bookkeeping  process  for  certifying 
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Type 

Error 

Description  of  Error  Sources 

(%) 

b 

5 

Standards 

Calibration 

0.12 

bi-  Standard  lab  calibration,  including  trace- 
ability  to  national  standards 

Hierarchy 

Calibration  System  Errors 

•  Reference  Source  . 

•  Sensor 

62-  Determination  of  reference  pressure 

-  Hysteresis  . 

(*7) 

S3-  Sensor  hysteresis  (combined  with  S7) 

Data 

-  Non-linearity  . 

(is) 

bs-  Sensor  non-linearity  (combined  with  bs) 

Acquisition 

—  Repeatability  . 

(S7) 

S5-  Sensor  repeatability  (combined  with  s?) 

Recording  System  Errors 

•  Sampling . 

* 

Sfi-  Data  variations  due  to  facility  or  engine 
instabilities 

•  Channel  . 

0.15 

0.15 

67  &  Sr-  Signal  conditioning,  electrical  calibra¬ 
tions  and  digital  systems 

Data 

Reduction 

Data  Processing  Errors . 

0.05 

bs-  Curve  fits  of  calibration  data 

. 

Installation  Effects  Errors 

•  Pressure  Probe  . . 

* 

br-  Design/fabrication  of  probes 

Other  Effects 

Environmental  Effects  Errors 

•  Temperature  . 

* 

sio-  Temperature  change  effect  on  sensor 

•  Vibration  . 

* 

611-  Vibration  effect  on  sensor 

Root  Sum  Square 

0.20 

0.15 

♦  Negligible  Error 


Table  7.4:  Pressure  measurement  error  sources  (Case  1) 


results,  Pre-test  error  analysis  allows  corrective  action  to 
be  taken  prior  to  the  test  to  reduce  uncertainties  that  ap¬ 
pear  too  high.  In  practice,  it  is  an  iterative  procedure  to 
tailor  the  entire  process  and  minimize  uncertainty  levels. 
Post-test  error  analyses,  which  are  based  on  actual  test 
results, 

•  permit  refinement  of  final  uncertainty  levels, 

•  check  for  consistency  of  redundant  measurements, 

•  identify  data  validation  problems. 


7.1.3  Uncertainty  Analysis 

Representative  uncertainties  were  determined  using  the 
methodology  defined  in  Section  3.7  and  Appendix  A,  Re¬ 
sults  are  presented  for  input  parameters  (both  measured 
and  calculated)  and  performance  parameters.  Represen¬ 
tative  uncertainties  for  measured  parameters  of  pressure, 
temperature,  force  (load  cell),  area  and  fuel  flow  rate  are 


presented  in  Tables  7.4  through  7.8.  Next,  the  uncer¬ 
tainties  of  calculated  input  parameters  were  determined 
by  error  propagation  using  appropriate  measured  parame¬ 
ter  uncertainties  and  Influence  Coefficients  for  airflow  rate 
(Table  7.9),  (Table  7.10)  and  equivalence  ratio  (Table 
7.11),  and  by  choosing  realistic  values  for  nozzle  discharge 
coefficient  and  stream  thrust.  Then,  the  uncertainties  of 
performance  parameters  were  determined  for  efficiencies 
based  on  characteristic  velocity  (c*),  vacuum  specific  im¬ 
pulse  (iJoc),  temperature  rise  (Tm  — T(2),  equivalence  ratio 
(<j))  and  combustion  chamber  pressure  loss  {pt‘i/pt2),  again 
by  error  propagation  of  appropriate  measured  and  calcu¬ 
lated  input  parameters  (Table  7.12). 


7.1.4  Influence  Coefficients  Comparisons 

A  sensitivity  analysis  was  performed  for  each  input  para¬ 
meter  (both  measured  and  calculated)  used  in  the  calcu¬ 
lation  of  combustion  efficiency.  Influence  coefficients  were 
established  for  each  input  parameter  by  numerically  per- 
turbating  the  perfomiance  equations  for  a  1%  difference 


Type 

Error 

Description  of  Error  Sources 

b 

i 

Standards 
Calibration  • 

o.rc 

6i-  Manufacturer  specification  of  wire  or  stan- 

Hierarchy 

dard  lab  calibration 

Calibration  System  Errors 

•  Level  . 

0.4’C 

62-  Reference  temperature  level 

Data 

•  Stability  . 

0.2‘’C 

Sa*  Reference  temperature  stability 

Acquisition 

Recording  System  Errors 

•  Sampling . 

o.rc 

64-  Data  variations  due  to  facility  or  engine 
instabilities 

•  Channel  . 

o.^c 

65-  Signal  conditioning,  electrical  calibrations 
and  digital  systems 

-  Data 
Reduction 

Data  Processing  Errors  . . 

0.28‘’C 

Sfi-  Curve  fits  of  calibration  data 

Installation  Effects  Errors 

•  Probe  Recovery  . 

♦ 

bj-  Probe  design  caused  by  radiation,  convec¬ 
tion,  etc. 

Other  Effects 

•  Conduction  Error  . 

* 

bs-  Heat  conduction 

•  Temperature  Gradients 

♦ 

69-  Temperature  gradients  along  nonhomoge- 
neous  thermocouple  wire 

Root  Sum  Square 

0.5‘‘C 

0.3“C’ 

*  Negligible  Error  with  proper  design/installation 


Table  7.5:  Temperature  measurement  error  sources  (Case  1) 


in  that  parameter,  while  keeping  all  other  parameter.s  con¬ 
stant  at  their  nominal  values  (using  the  procedure  speci¬ 
fied  in  Appendix  A,  Section  A. 2. 2).  Results  are  presented 
in  Table  7.12. 
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Type 

Error 

Description  of  Error  Sources 

(%) 

b 

5 

Standards 

Calibration 

0.12 

bx-  Standard  lab  calibration,  including  trace- 

Hierarchy 

ability  to  national  standards 

Calibration  System  Errors 

•  Off-Axis- Effects  . 

* 

62-  Force  measurement  on  axis  different  from 
centerline 

•  Tare  Correction 

-  Hysteresis  . 

0.1 

S3-  Force  measurement  system  hysteresis 

Data 

—  Non-repeatability  .... 

(sa) 

S4-  Sensor  non-repeatability  from  repeat  cali¬ 
brations  (combined  with  53) 

Acquisition 

-  Non-linearity  . 

0.05 

65-  System  calibration  non-linearity 

Recording  System  Errors 

•  Sampling  . 

« 

S6-  Data  variations  due  to  facility  or  engine 
instabilities 

•  Channel  . 

0.12 

S7-  Signal  conditioning,  electrical  calibrations 
and  digital  systems 

Data 

Reduction 

Data  Processing  Errors  . . . 

0.1 

tg-  Curve  fits  of  calibration  data 

Installation  Effects  Errors 

•  Stand  Alignment  . 

+ 

-  69-  Misalignment  of  engine  and  data  load  cell  I 
force  vectors  j 

•  Load  Cell  . 

* 

610-  Shift  in  load  cell  calibration  caused  by 
attachments 

Pressure  Effects  Errors 

•  Load  Cell  . 

* 

611-  Cell  pressure  change  on  load  cell 

Other  Effects 

•  Test  Cell  . . 

+ 

6j2-  Cell  pressure  change  on  test  cell  wall 
ground  ' 

•  Service  Lines  . 

0.05 

6i3-  Line  pressure  change  on  tare  forces 

Temperature  Effects  Errors 

•  Load  Cell  . 

614-  Temperature  change  on  load  cell 

•  Service  Lines  . 

(bie) 

615-  Line  temperature  change  on  tare  forces 
(combined  with  tie)  | 

•  Thrust  Stand  . 

0.1 

tie-  Thermal  growth  of  stand 

•  Vibration  Error  . . 

* 

si7-  Vibration  of  load  cell 

•  Scrub  Drag  Error . 

* 

sig-  Secondary  airflow  external  drag 

Root  Sum  Square 

0.2 

'0H51 

»  Negligible  Error 


Table  7.6:  Scale  force  measurement  error  sources  (Case  1) 
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Type 

Error 

(%) 

Description  of  Error  Sources 

Standards 

Calibration 

Hierarchy 

6 

s 

bi-  Standard,  lab  calibration  for  measurement 
instrument,  including  traceability  to  national 
standards 

63-  Determining  cross-sectional  area 

63-  Difference  in  measurement  and  test 
temperatures,  including  effect  of  temperature 
error 

■■ 

Lj 

(ither  Effects 

•  Diameter  Measurement  Error  .... 

•  Temperature  Compensation  Error 

Root  Sum  Square 

0.11 

Oj 

. 

Table  7.7;  Area  measurement  error  sources  (Case  1) 
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Type 

1 

Error 

Description  of  Error  Sources  ’ 

(%) 

b 

s 

Standards 

Calibration 

0.1 

fc]-  Standard  lab  calibration,  including  trace- 

Hierarchy 

ability  to  national  standards 

■ - - - * - 1 

Calibration  System  Errors 

•  Cal  Fluid  Properties  . 

0,15 

62-  Mismatch  between  calibration  and  test 
fluids 

•  Flowmeter  Repeatability 

0.12 

i>3-  Non-repeatability  from  repeat  flowmeter 
calibrations 

Data 

Recording  System  Errors 

Acquisition 

1 

•  Sampling . ■. . 

* 

S4-  Data  variations  due  to  facility  or  engine 
instabilities  1 

•  Channel  . 

0.17 

S5-  Signal  conditioning,  electrical  calibrations 
and  digital  systems 

Data 

Reduction 

Data  Processing  Errors  . 

0.1 

bs-  Curve  fits  of  calibration  data  j 

Installation  Effects  Errors 

•  Cavitation  . 

67-  Insufficient  static  pressure  in  flowmeter 

•  Turbulence  . 

* 

fcg-  Sharp  bends,  etc.,  upstream  flowmeter 

•  Meter  Orientation  . 

♦ 

fcg-  Orientation  difference  from  calibration  to  , 
test 

Other  Effects 

Environmental  Effects  Errors 

•  Temperature 

—  Flowmeter  . 

•  * 

frio-  Ambient  temperature  change  on 
flowmeter 

—  Viscosity  . 

0.12 

0.05 

5ii  fe  Sii-  Determination  of  test  fluid  viscosity 

-  Specific  Gravity  . 

0.1 

0.05 

b\2  &  *12'  Determination  of  test  fluid  specific 
gravity 

•  Vibration  . 

t 

Si3-  Vibration  on  flowmeter 

•  Pressure  . 

* 

bi4-  Ambient  pressure  change  on  flowmeter 

Root  Sum  Square 

0.28 

0.18 

*  Negligible  Error 


Table  7.8'.  Fuel  flow  measurement  error  sources  (Case  1) 


Basic  Measurements 

Air  Mass  Flow  Rate  (maiv) 

Input  1 

Bias 

Precision 

Influence 

1  Bias 

Precision 

Parameters 

Limits 

Index 

Coefficients 

Limits 

Index 

1  1 

flv 

5,' 

IC 

'  BilC 

Sk  =  SilC 

(%) 

(%) 

(%/%) 

1  (%) 

(%) 

1 

Pi  ' 

0.2 

0.15 

m 

T,  ! 

0.2 

0.15 

0.5 

0.1  1 

0.07 

B 

>1  ! 

0,05 

0 

1.0 

'  05  1 

0 

cd 

0,4 

0 

1.0 

1  ! 

0 

B  =  0,46  %  1 

5  =  0.17% 

U  =  1 

Table  7.9:  Error  propagation  for  air  flow  rate  (Case  1) 


Basic  Measurements 

Process  7p  , 

Input 

Parameters 

1 

Bias 

Limits 

Bi 

i%) 

Precision 

Index 

(%) 

Influence 

Coefficients 

IC 

(%/%) 

Bias 

Limits 

Fi  =  BilC 
{%)  ■ 

Precision 

Index 

5i  =  SilC 

i%) 

1 

Tt2 

0.20 

0,15 

0.024 

0.0048 

0.0036 

0.40 

0.25 

AifAi 

0.07  1 

0 

0 


38 


0.00056 


B  =  0.038 ' 


0 


0.24 


0 


5  =  0.24' 


{/  =0.52' 


Table  7.10:  Error  propagation  for  process  (Case  1) 


Basic  Measurements 


Input 

Parameters 

/■ 


Bias 

Limits 

Bi 

(%) 


0.29 


0.28 


Precision 

Index 

5. 

(%) 


0.17 


8 


Equivalence  Ratio  (i^) 


Influence 

Coefficients 

IC 

{%/%) 


1.0 


Dials 
Limits 
Bt  =  BilC 
{%) 


0.29 


8 


Precision 
Index 
5t  =  SilC 
{%) 


0.17 


8 


B  -  0.40  % 

II 

1 

5  =  0.25% 

L^  =  0.9% 


Table  7.11:  Error  propagation  for  equivalence  ratio  ^  (Case  1) 


•  P4.  yp,t 

•  p4  without  7 

•  Bi,  7p,j 

•  F,  without  7 


p 

Ai  1 

CDS  //a 

-0.2 

1 

1 

1.2  -2.6 

-0.2 

1 

1  -2.6 

<0.1 

0.1 

0  -2.7 

0 

1  0.1 

0  -2.9 

Influence  Coefficients 


-^5  Ps  Pamb 


Errors 


[%]  [%]  [%] 


•  P4.  Jp.i 

•  p4  without  7 

•  Fs)  7p,j 

•  F5  without  7 


0.9  0 


•  P4,  7p,j 

•  P4  without  7 

•  Bi,  7p  j 

•  F5  without  7 


•  P4.  7p,J 

•  p4  without  7 

•  ^5,  7p,i 

•  F5  without  7 


Table  7.12;  Sensitivity  analysis  and  uncertainty  for  performance  parameters  (Case  1) 
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7.2  LFRJ  Performance  with  Viti¬ 
ated  Air  Heater  (Case  2) 

Sample  calculations  are  provided  for  a  methane-fueled  air 
vitiator  with  make-up  oxygen,  The  geometrical  data  and 
ramjet  measured  test  data  are  the  same  as  for  Case  1  ex¬ 
cept  for  air  and  vitiator  flowrates. 


7.2.1  Calculation  Procedure 

7. 2. 1.1  General  Test  Information 

The  vitiated  air  consists  of  the  following  ingredients  with 
their  associated  composition  and  heat  of  formation  at 
298, 15A". 


•  Pressures  and  Temperatures 


Pamb 

= 

101300  Pa 

P2 

= 

650200  Fa 

P4 

= 

568800  Pa 

P5 

= 

341203  Pa 

Pb 

= 

78065  Pa 

Tn 

= 

606  K 

Tf 

298,15  K 

7.2. 1.4  Preliminary  Calculations 

The  following  preliminary  calculations  must  be  made  to 
provide  inputs  for  the  efficiency  parameter  equations  in 
Chapter  6; 


•  Ramjet  fuel',  €10^20  kerosene 

•  Vitiator  fuel;  CHu 

•  Air;  A^iss.^C^i.geAro  s34Co.o3i4 


•  Ramjet  fuel  enthalpy 

hj  =  book  value  at  measured  Tf 

hj  =  -2016  60  = -67607ca//mole 

•  Nozzle  stream  thrust  (Eqs.  C.26  -  C-25) 


7.2. 1.2  Geometric  Data 
•  Flow  areas 


Fmtcs  ~~  ^LC  ^tare  ■ 

=  13400  -  5000 
=  8400  N 


A2  — 

0.010325 

Ts  —  F meat  -^biPb  Pamb) 

A4  = 

0.022698  77)2 

■=  8400  -  0,004304  (78065  -  101300) 

=  8500  N 

•  Exhaust  nozzle  (Fig.  5.1 

and  C.l) 

•  Station  2  conditions 

As  = 

0.012668  m2 

The  properties  of  the  vitiated  air  are  determined  us¬ 

CDS  = 

0.996 

ing  one  of  the  methods  of  Section  4. 4. 5. 3.  Method  3 

Ah  = 

0.004304  77)2 

v/as  used  in  this  example.  The  aerothermochemical 
equilibrium  code  was  run  to  determine  the  properties 
of  the  vitiator  products  at  the  measured  Tij,  In  this 

7. 2. 1. 3  Measured  Test  Data 

case,  the  thermodynamic  state  of  the  vitiator  prod¬ 
ucts  is  specified  by  p(2  and  Tt2,  thus,  the  enthalpies 

•  Mass  flows 

used  in  the  input  file  are  irrelevant.  Samples  of  the 
code  input  and  output  are  shown  in  Appendix  D  (Fig 

rhv/  = 

QMlkg/s  ■ 

D.8). 

moa  = 

0.262  kg/s 

o  inputs:  pt2  (assumed  equal  to  the  measured  P2), 

6.383  kg/s 

Tt2 ,  vitiator  mass  flows 

TTlf  = 

0.311  kg/s 

0  outputs; 

•  Thrust  and  forces 

M2  =  2S.90i  kg/ kmol 

72  =  1.3687 

Plc  - 

=  13400  N 

/l2  =  ht2 

Fjars  - 

~  5000  N 

=  —69.041  kJ /kg  =  —477.0  cai/mole 
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composition  (mole  fractions): 

0  calculate  the  stoichiometric  fuel/air  ratio  using 

XN,  =  0.74339 

Eq.  4.2  . 

(L\  = 

Xo^  =  0.20946 

J  etp 

XHtO  =  0.02531 

my 

XcOa  ■  =  0.01295 

ffiflir  +  +  rhoj 

X-Ar  =  0.00889 

0.311 

6.383  +  0.047  +  0.262 

Vitiated  air  composition: 

- 

0.0465 

Nj.4a680o.470lCo.01295Ho.05062Aro. 00889 

(my  /mair)^xp 

V  ^  /  sfoi'cA 

<!> 

Ri  = 

Af  2 

- 

0.0465 

8314.51 

0.7201 

28.908 

0.06457 

=  287.62  7-^ 

•  Mass  flows 

kg  K 

m2  = 


•  Chamber  total  pressure  (pm) 


6.383  +  0.047  +  0.262 
6.692  kg/s 


o  run  aerothermochemical  equilibrium  code  using 
P4  as  a  first  approximation  to  pt4  (Fig.  D.9) 

*  inputs:  P4,  Ail{As  C£i5),  h-i,  hj,  m2,mj,  air 
composition 

*  outputs: 

7.2. 1.5 


1714  =  7775  =  rh2  + 

=  6.692  +  0,311 
=  7.003  kg/s 


Performance  Calculation 


The  performance  parameters  can  now  be  calculated  by 
employing  the  equations  of  Chapter  6  as  summarized  in 
Table  7.1.  Only  method  2  is  presented,  Results  are  pre¬ 
sented  in  Table  7,13. 


— =  1.0791 
J  th 

o  compute  pi4  using  (Eq.  6,5): 


=  568800  x  1.0791 
=  613792  Pa 

0  rerun  aerothermochemical  equilibrium  code  at 
Pt4  =  6l3792Pa  to  get  more  accurate  results 
(Fig.  D.IO) 

*  input;  pt4,  A^/i^As  Cds),  ^2,  hj,  m2,  m;, 
air  composition 

(all  unchanged  except  for  p!4) 

*  output; 


Pf4  \ 

/  (7i 

= 

1.0791 

Ps.tK 

= 

339950  Pa 

= 

2044.96  A' 

Ts.iA 

= 

1819.27  A' 

cTh 

1167  m/s 

‘vac,th 

= 

1457.1  Ns/kg 

M4 

= 

28.855  kg/kmole 

4> 

= 

0.7201 

1.  Efficiency  based  on  characteristic  velocity 
From  Eq.  6.2 

.*  _  PtiAscos 

^eip  —  ■ 

m4 

613792  X  0.012668  x  0.996 
7.003 

=  1105.87  m/s 

Using  Eq.  6.1 


1105.87 
1167 
=  0,9476 

2.  Efficiency  based  on  vacuum  specific  impulse 
From  Eq.  6.12 


=  1380.77  Ns/ kg 
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Method 

1  . 

2 

3 

4 

Calculation  Procedure 

P4  using  7 

P4  without  using  7 

Es  using  7 

Fs  without  using  7 

Mi 

z 

0.352 

Pt2 

Pa 

706963 

7p,,  1’‘  iter. 

- 

- 

M4 

- 

- 

• 

Pm 

Pa 

613792 

i?4,exp 

JfkglK 

288.2 

C* 

mfs 

1106 

^4, exp 

K 

1836 

7p,,  2“  iter. 

- 

i* 

‘'t/ae.erp 

Ns/ kg 

1381 

<Pi> 

- 

0.573 

{f  /^^sioieh 

- 

0.0646 

4>in3 

. 

0.720 

Vc* 

- 

0.948 

.  _ 

- 

0.948 

Vat 

0.855 

Vt 

- 

0.796 

Pm/P«2 

0.868 

(Pt2  ~Pi4)IPt2 

- 

032 

Table  7.13:  LFRJ  performance  with  vitiated  air  (Case  2) 


From  Eq.  6.8 


2” 

1380.77 

1457.1 
=  0.9476 


3.  Efficiency  based  on  temperature  rise 
Using  Eq.  6.18  with  It4^exp  =  -R-t.ifc 

M4 

8314.51 

28.855 

=  288.15  J/Fs/ZC 


^  <4, exp 


/  Tt^R^  ^  ^ezp 
\  J  th  ^‘i.txp 

/2044.96X  288.15'\  1105.87* 
V  1167=  )  288.15 

1836  33  K 


Using  Eq.  6.15 


tlAT 


'Rj'i.zxp  —  'R,2 

Tn^th  ~'Tt2 

1836.33-  606 
2044.96-  606 
0.8550 


4.  Efficiency  based  on  equivalence  ratio 

The  equivalence  ratio  (i^^)  necessary  to  theoretically 
produce  the  measured  =  1105. 87  m/s  is  deter¬ 
mined  from  a  plot  similar  to  Figure  7.1 ,  but  calculated 
for  vitiated  air. 

<1.1,  =  0.5727 

Using  Eq.  6.21. 


^inj 


=  (— ) 


txp  (//^)a<oich 

/0.311\  1 

V6-692/  0.06457 
=  0.7197 


Using  Eq.  6.20 


<i>i 

4>in} 

0.5727 

0.7197 

0.7957 


5.  Calculation  of  combustor  total  pressure  loss 


•  calculate  the  Mach  number  at  station  2; 
using  Eq.  6.25 

m2  I  RjT, 2 

P2A2  T2 


2\A 


Mi 


3687  -  1 


M|  = 
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7.3  Ducted  Rocket  Mass  Flow 
Rate  (Case  3) 

The  fuel  mass  flow  of  a  ducted  rocket  can  be  determined 
by  equation  5  4  for  a  choked  gas  generator,  a  constant 
injection  throat  area  and  a  constant  c*: 

Fig.  7.2  shows  a  typical  gas  generator  pressure-time  trace 
(two  transducer  outputs  are  shown  in  the  diagram)  for  an 
end  burning  grain.  The  integration  of  the  gas  generator 
pressure  over  time  is  done  between  the  points  where  the 
gas  generator  is  choked  relative  to  the  ramburner.  For  the 
case  presented,  the  integration  limits  were  at  0.3  AfPa  for 
the  start  and  at  0.5  MPa  for  the  end.  A  correction  of  the 
pressure  integral  may  be  needed  to  take  into  account  the 
mass  flow  injected  at  unchoked  conditions  during  burnout. 

The  integral  derived  from  the  pressure-time  trace  of  Fig. 
7.2  amounts  to  69.2994  MPa  x  s.  Together  with  the  used 
propellant  mass  (initial  propellant  mass  minus  residues  re¬ 
tained  in  the  gas  generator  after  burnout)  of  5. 053  ig  and 
using  Eq.  5.4,  a  mass  flow  given  by  Fig.  7.3  is  calculated. 

It  should  be  mentioned  that  a  period  of  increasing  pressure 
between  ignition  and  the  stationary  level  (e.g.,  during  the 
first  four  seconds  of  the  trace  in  Fig.  7.2)  can  be  attributed 
to  a  rising  c*  or  temperature.  This  leads  to  an  error  in  the 
mass  flow  determination  by  underestimating  the  mass  flow 
in  the  phase  of  rising  pressure  and  overestimating  the  mass 
flow  in  the  phase  of  constant  pressure. 


Pressure  progressivity  during  or  at  the  end  of  the  station¬ 
ary  burning  phase  may  be  attributed  to' 

•  intentional  increases,  of  the  propellant  burning  surface 
■  area  Induced  by  grain  geometry 

•  accidental  increases  of  the  propellant  burning  surface 
area  induced  by  nonuniform  (e.g.,  conical)  burning  or 
voids  in  the  grain 

•  clogging  of  the  gas  generator  nozzle  throat 

The  first  two  phenomena  do  not  affect  the  accuracy  of  the 
mass  flow  evaluation  according  to  Eq.  5.4  as  long  as  the  c* 
does  not  vary  within  the  range  of  the  pressure  progressiv¬ 
ity.  The  pressure  rise  at  the  end  of  burning  (Fig.  7.2)  was 
caused  by  voids  afthe  bottom  of  the  end  burning  grain 

Any  significant  variation  of  the  gas  generator  nozzle  throat 
area  denies  the  applicability  of  Eq.  5.4  for  mass  flow  evalu¬ 
ation.  Even  when  corrections  are  conceivable  to  take  into 
account  a  varying  gas  generator  throat  area,  the  uncer¬ 
tainty  of  the  mass  flow  evaluated  from  the  pressure  inte¬ 
gration  will  be  high. 

The  fuel  mass  flow  can  be  determined  alternatively  by  Eq. 
5.5  for  a  choked  or  unchoked  gas  generator. 

In  the  example  the  constant  burning  surface  area  of  the 
end  burning  grain  is  0.0201  and  the  propellant  density 
is  1510  kg/m^.  The  burning  rate  as  a  function  of  pressure 
is  given  by  Fig.  7.4.  Thus,  the  mass  flow  can  be  evaluated: 

e.g.  at  10s 

Pig  =  4.74MPa 
rj  =  ll,5mm/s 
rh}  =  0.349ip/s 

For  this  procedure  of  mass  flow  evaluation  it  is  essentia! 
that  the  ballistic  data  (like  given  in  Fig.  7.4)  are  valid 
for  the  full-scale  grain  and  its  geometry.  Burn  rate  data 
are  usually  evaluated  by  firing  small  motors  or  by  burning 
strands  in  a  combustion  bomb.  The  applicability  of  strand 
burn  rate  data  to  full-scale  motors  is  generally  poor.  Burn 
rates  determined  with  small  ballistic  test  motors  usually 
show  a  better  comparability  to  the  results  from  full-scale 
motors.  Nevertheless,  the  possibilities  that  differences  be¬ 
tween  model  and  full-scale  grains  may  affect  the  burn  rate 
(thermal  effects,  erosion,  etc.)  should  be  carefully  con¬ 
sidered  before  applying  the  above  method  for  mass  flow 
evaluation. 

Finally,  it  has  to  be  mentioned  that  thermal  erosion  of 
the  propellant  Insulation  (liner,  boot)  may  contribute  to 
the  gas  generator  mass  flow;  especially  for  end  burning 
grains  where  the  effect  may  be  significant.  Mass  flow  from 
the  insulation  consumed  is  included  when  the  mass  flow 
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Time  [si 


Figure  7.2;  Gas  generator  pressure  (Case  3) 


Time  [s] 


Figure  7,3:  Gas  generator  pressure  and  calculated  fuel  mass  flow  (Case  3) 


is  determined  by  the  pressure  integration  method.  How-  consumption  of  the  grain  insulation  be  checked  in  order 
ever,  the  contribution  of  insulation  to  the  gas  generator  to  be  aware  of  the  magnitude  of  the  possible  errors  of  the 
mass  flow  may  not  be  constant  during  the  burn  time  and  gas  generator  mass  flow  evaluation, 
the  heating  value  usually  differs  from  that  of  the  propel¬ 
lant,  The  method  for  determining  the  gas  generator  mass 
flow  based  on  the  ballistic  data  does  not  include  possible 
insulation  mass  flow. 


Procedures  to  correct  for  the  mass  flow  corresponding  to 
the  consumed  insulation  have  to  be  found,  depending  on 
the  individual  case.  It  is  strongly  recommended  that  the 


4S 


P  fMPeJ  7h»I«l 


Figure  7.4;  Ballistic  data  of  the  ducted  rocket  propellant  Figure  7-5r  SFRJ  burn  time  determination  (Case  4) 
(Ceise  3) 

7.4  Solid  Fuel  Ramjet  Mass  Flow 
Rate  (Case  4) 


The  fuel  flow  rate  for  the  SFRJ  is  usually  a  time  aver¬ 
aged  value  determined  from  the  mass  of  consumed  fuel 
and  combustion  burn  time  according  to  equation  5  3: 


m, 


■ 

h  . 


An  example  is  provided  for  clarity. 

m/,i  =  0.684^3 

myj  =  0.305  /cfl 

/j  =  9.55  s 

determined  from  Fig.  7.5. 

The  initial  and  final  fuel  masses  were  determined  with  a 
balance  or  scale.  The  combustion  burn  time  was  deter¬ 
mined  from  a  pre.ssure  versus  time  plot  cis  shown  in  Fig, 
7.5.  The  times  at  combustor  ignition  and  burnout  were 
determined  using  75%  of  the  difference  in  pressures  for 
combustion  and  no  combustion.  This  method  provides  a 
reasonable  alternative  to  an  integral  approach  and  permits 
quick  determination  of  the  burn  time.  Thus,  using  equa¬ 
tion  5  3,  the  value  for  fuel  mass  flow  rate  was  determined. 

0.684 /tfl- 0.305  fcp 
~  ■  9.55  s 

=  0.0397  kg(s 
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8  Air  Heaters 


Aerodynamic  compression  through  the  supersonic  inlet  of 
a  ramjet  causes  an  increase  in  the  static  temperature  of  the 
air  entering  the  combustor.  Since  combustor  inlet  Mach 
numbers  are  quite  low,  the  static  temperature  is  essen¬ 
tially  equal  to  the  stagnation  temperature.  At  a  flight 
Mach  number  of  4,  the  combustor  inlet  temperature  is  ap¬ 
proximately  lOOO/i,  at  Mach  6  it  is  near  2000/f,  and  at 
Mach  8  the  temperatures  are  in  the  3000A1  range.  The 
precise  amount  of  heating  required  to  simulate  true  flight 
conditions  also  varies  with  altitude.  In  order  to  examine 
combustion  performance  at  realistic  flight  conditions  these 
combustor  inlet  temperatures  must  be  reproduced  in  test 
facilities,  and  a  method  for  producing  this  heated  air  must 
be  found. 


8.1  Real  Gas  Effects 

Since  air  exhibits  non-ideal  gas  characteristics  at  elevated 
temperatures  (thermodynamic  properties  vary  with  tem¬ 
perature),  prediction  of  air  total  temperature  require¬ 
ments  for  a  test  facility  must  include  consideration  of  real 
gas  effects.  Figure  8.1  shows  how  combustor  inlet  total 
temperature  varies  with  flight  Mach  number  and  altitude. 
Figure  8.2  shows  the  overprediction  of  the  air  tempera¬ 
ture  that  results  when  assuming  constant  properties  for 
air  as  found  in  standard  tables  for  supersonic  flow  using 
a  specific  heat  ratio  of  1.4.  These  tables  should  not  be 
used  to  determine  the  correct  total  temperature  for  sim¬ 
ulation.  The  correct  values  of  Figure  8.1  were  obtained 
by  calculating  enthalpy  from  the  given  Mach  number  and 
altitude  conditions  and  using  an  aerothermochemical  equi¬ 
librium  code  to  solve  for  the  air  temperature  assuming  an 
isentropic  compression.  This  method  takes  into  account 
species  concentration  changes  as  well  as  specific  heat  vari¬ 
ations  to  predict  the  total  temperature  required  of  the  test 
facility  air  heater. 


8.2  Heater  Requirements 

The  ideal  heater  should  deliver  air  over  a  wide  operational 
envelope  of  mass  flow,  temperature  and  pressure  and  be 
able  to  support  a  wide  range  of  test  durations.  It  should 
also  deliver  good  air  flow  quality  (uniform  temperature 
profile  and  low  turbulence  levels)  and  be  free  of  air  con¬ 
taminants.  A  heater  should  be  easy,  safe  and  affordable  to 


Figure  8.1:  Influence  of  Mach  number  and  altitude  on  total 
temperature 


Altitude  [km] 

Figure  8.2:  Real  gas  and  Mach  number  effects  on  air  total 
temperature 
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use.  Ignition  and  control  of  the  output  flow  should  be  sim¬ 
ple  and  dependable  and  steady  state  conditions  achieved 
quickly.  The  heater  should  be  small  and  have  low  oper¬ 
ation  and  maintenance  costs.  The  fuel  or  energy  source 
should  also  be  economical,  safe  and  available. 

No  heater  meets  all  of  these  criteria,  and  each  test  facility 
has  its  own  priorities.  Therefore,  compromises  must  often 
be  made  to  match  the  requirements  of  a  particular  facility 
with  the  capabilities  of  available  heaters. 


8.3  Heater  Types 

Different  air  heating  methods  are  used  in  ramjet  test¬ 
ing.  They  may  be  grouped  into  three  broad  categories’, 
combusting  heaters  (vitiators),  non-combusting  heaters 
(including  heat  exchangers),  and  combinations  of  these 
heaters. 


8.3.1  Combusting  Heaters  (Vitiators) 

Combusting  heaters  heat  air  directly  with  a  fuel-oxidizer 
reaction.  FUel  is  burned  in  the  air  stream,  often  in  a  jet 
engine  style  combustor.  The  exhaust  gases  are  used  for 
ramjet  testing  after  the  consumed  oxygen  is  replenished. 

The  main  advantages  of  combusting  air  heaters  are  the 
low  cost,  because  of  the  low  fuel  flow  required,  and  their 
simplicity  of  operation  and  maintenance.  The  disadvan¬ 
tages  include  the  effects  of  air  heater  combustion  products 
on  subsequent  ramjet  combustion  and  the  change  of  air 
properties  like  molecular  weight. 


8.3.2  Non-Combusting  Heaters 

Non-combusting  heaters  avoid  contamination  of  the  air 
stream  with  combustion  products  and  deliver  clean  air  to 
the  ramjet  combustor. 

In  these  heaters,  air  flows  through  a  heat  exchanger.  The 
heat  source  may  be  electrical  resistance  or  combustion  of 
a  separate  fuel  and  oxidizer.  The  heat  exchangers  can 
be  quite  large  and  exit  temperatures  are  limited  by  the 
material  properties  of  the  heat  exchanger. 

Heat  storage  devices  are  heat  exchangers  of  high  thermal 
mass  that  are  gradually  heated  to  operating  temperature 
and  during  a  test  run  give  up  stored  thermal  energy  to 
air  passing  through  them.  Commonly,  large  vessels  filled 
with  ceramic  or  metal  pebbles  are  used,  heated  by  either 
electrical  resistance  or  hot  combustion  gases.  These  stor¬ 
age  heaters  can  heat  air  to  higher  temperatures  than  con¬ 
ventional  heat  exchangers  because  the  pebbles  are  more 


resistant  to  thermal  damage  than  fragile  heat  exchanger 
tubes.  Nevertheless,  repeated  thermal  expansion  and  con¬ 
traction  can  rub  the  pebbles  together  and  introduce  dust 
particles  into  the  air  stream.  Heater  life  span  is  lengthened 
by  reducing  the  number  and  magnitude  of  these  thermal 
cycles.  The  main  challenge  in  designing  heat  storage  de¬ 
vices  is  obtaining  a  constant  output  temperature  for  the 
required  run  time  and  range  of  test  conditions.  This  usu¬ 
ally  results  in  large  heaters  with  a  high  heat  capacity  so 
that  only  a  small  fraction  of  the  energy  is  extracted  during 
a  run. 

Electric  arc  heaters  heat  air  through  the  release  of  en¬ 
ergy  produced  by  an  electric  arc  between  two  electrodes. 
Arc  heaters  are  capable  of  producing  very  high  air  tem¬ 
peratures.  The  facility  capabilities  depend  primarily  on 
the  limits  of  the  anode  and  cathode  producing  the  arc. 
Due  to'  the  extremely  high  temperatures  produced  by  the 
arc,  ionized  species  are  created  that  react  to  form  NOx 
and  other  undesired  constituents  that  contaminate  the  air 
stream.  The  presence  of  these  contaminants,  and  the  fact 
that  oxygen  dissociation. begins  at  approximately  2500A', 
sets  the  upper  limits  of  combustion  testing  in  arc  facilities. 
Above  this  temperature,  care  must  be  taken  to  account  for 
the  chemical  effects  of  the  contaminants  and  dissociated 
species  on  combustion.  As  with  any  electric  heater  the 
high  power  requirements  of  arc  heaters  make  them  very 
expensive  to  operate. 


8.3.3  Combination  Heaters 


Combination  heaters  use  a  combination  of  the  previously 
discussed  methods  to  take  advantage  of  the  characteristic 
strengths  of  one  method  to  offset  the  weaknesses  of  the 
other.  For  example,  one  could  use  a  vitiator  to  boost  the 
temperature  from  an  electrically  powered  heat  exchanger. 
This  would  allow  higher  temperatures  without  damaging 
the  heating  elements  in  electric  heater,  and  would  deliver 
lower  levels  of  combustion  products  when  compared  to 
pure  combustion  heating.  In  addition  this  combination 
allows  temperature  variation  during  a  run  (transient  sim¬ 
ulations). 


8.4  Special  Considerations  for  Vi¬ 
tiators 

Since  vitiators  are  so  widely  used,  the  following  sections 
address  some  important  issues  concerning  their  operation. 
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8.4.1  Make-up  Oxygen 


Hence 


In  heating  the  air,  combustion  in  a  vitiator  depletes  part  of 
the  available  oxygen.  In  order  to  conduct  combustion  ex- 
periments  the  oxygen  content  of  the  flow  must  be  restored 
to  the  proper  percentage. 


Substituting: 


Calculation  of  make-up  oxygen  flow  rates  is  often  made 
to  yield  a  mole  fraction  of  oxygen  in  the  vitiator  products 
equal  to  that  in  atmospheric  air  (0.2095  for  the  first  air 
composition  in  Table  4.2).  To  get  this  mole  fraction  the 
stoichiometric  chemical  reaction  for  the  vitiator  fuel  and 
oxidizer  combination  must  be  examined.  For  a  hydrogen 
-  air  vitiator  the  reaction  is: 


y.=  0.1696 


l-t-2j/ 


1  -i-2  X  0.1696 


0.6696 


0.6696 


moles  O2 
mole  H2 
31.9988 


2.016 
10.628ity  Os/fcff 


H2  -I-  XO2  -I-  Air  =  H2O  4-  yOz  +  Air 

where  x  is  the  number  of  moles  of  make-up  oxygen  and  y 
is  the  number  of  moles  of  oxygen  in  the  vitiator  products 
per  mole  of  hydrogen  vitiator  fuel.  Notice  that  because  the 
air  in  the  reaction  already  has  the  correct  mole  fraction, 
calculation  of  the  make-up  oxygen  is  only  a  function  of 
the  fuel  flow.  To  get  a  0.2095  oxygen  mole  fraction  among 
the  products  other  than  air  the  following  equation  must 
be  satisfied: 

0.2095  =  y/(l  +  y) 


This  make-up  oxygen  mass  ratio  is  12.5%  lower  than 
the  12.142  ratio  that  was  computed  on  a  molar  basis. 
The  ratio  of  oxygen  in  the  ’non-air’  vitiator  products  is 
0.1696/(1  -I-  0.1696)  =  0.1450  instead  of  0  2095  as  in  air. 
This  lowers  the  partial  pressure  of  oxygen  in  the  entire 
vitiator  products. 

The  make-up  oxygen  may  be  mixed  into  the  flow  at  any 
point  upstream  of  the  ramburner.  However,  it  is  usually 
helpful  to  add  the  oxygen  upstream  of  the  vitiator  to  en¬ 
sure  good  mixing  and  help  raise  the  heater’s  efficiency  and 
broaden  its  operating  range. 


This  gives  y  =  0.265  for  the  hydrogen-air  vitiator. 

Using  an  oxygen  balance  to  define  r:  8.4.2  Air  Contaminants 


2r=  l-|-2y=  1.53 


thus,  X  =  0.7650. 
Converting  to  a  mass  ratio; 


/ 0.765  mole  O2  ^ 

\  .  1  mole  H2  J 

1  2  ) 

12.142 


kg  O7 

kg  H2 


Note  that  this  method  is  not  equivalent  to  adding  oxygen 
to  the  vitiator  exhaust  to  yield  an  oxygen  mass  fraction 
equivalent  to  atmospheric  air  (0  2315).  Maintaining  the 
mole  fraction  of  oxygen  preserves  its  partial  pressure,  and 
therefore  reaction  rates  which  may  be  important  at  low 
combustor  pressures  or  high  combustor  Mach  numbers. 
If  the  mass  fraction  of  oxygen  is  maintained  instead,  the 
oxygen  available  to  react  in  the  combustor  will  be  less. 


Using  the  above  oxygen  balance  equation  (2i  =  1  -|-  2y), 
and  the  oxygen  mass  fraction  in  atmospheric  air  (0.2315): 


0.2315 


31.9988y 

18.015 -i- 31. 9988y 


The  use  of  vitiated  air  instead  of  ideal  air  to  study  ramjet 
combustion  in  connected-pipe  testing  will  have  an  effect 
not  only  on  the  experimental  performance  but  also  on  the 
theoretical  performance  prediction. 

The  most  important  problem  associated  with  vitiated  air 
tests  is  the  influence  of  the  contaminants  on  the  combus¬ 
tion  process.  Unreacted  vitiator  fuel  will  add  energy  to 
the  ramjet  and  species  such  as  CO  and  H2O  can  change 
ignition,  combustion  efficiency  and  flame  holding  charac¬ 
teristics  measured  in  the  combustor. 

Concerning  the  theoretical  performance,  contaminants 
complicate  the  problem  of  calculating  properties  such  as 
the  molecular  weight,  specific  heat  ratio  and  the  enthalpy 
of  the  air  entering  the  combustor.  However,  if  the  air 
heater  operates  efficiently  these  thermodynamic  proper¬ 
ties  can  be  calculated  with  computer  codes,  assuming  the 
heater  products  are  at  equilibrium  by  the  time  they  enter 
the  combustor. 

In  the  following  the  effect  of  vitiated  air  on  both  theoreti¬ 
cal  as  well  as  experimental  performances  will  be  discussed 
in  more  detail. 
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Figure  8.3:  Effect  of  vitiator  fuel  on  Ms  when  comparing 
vitiated  air  with  ideal  hot  air  [34] 

8,4.2. 1  Effects  of  Vitiated  Air  and  Vitiator  Fuel 
Type  on  Ramjet  Theoretical  Performance 

The  effects  of  vitiated  air  on  ramjet  theoretical  perfor¬ 
mance  were  investigated  using  the  NASA  CET89  aero- 
thermochemical  equilibrium  code.  The  study  ([34])  quan¬ 
tified  the  effect  of  vitiated  air  on  solid  fuel  ramjet  per¬ 
formance  for  hydrogen,  methane  and  natural  gas  vitiation 
while  maintaining  the  oxygen  mole  fraction  of  0.2095.  The 
study  was  performed  with  an  initial  air  temperature  of 
298.15  K  and  the  fuel  temperature  was  varied  between 
50  K  and  200  K.  The  more  important  results  are  pre¬ 
sented  here  in  Figs.  8.3  to  8.6.  It  can  be  observed  that 
variations  in  ramjet  performance  increase  with  increasing 
combustor  inlet  temperature.  Fig.  8.6  shows  that  ramjet 
specific  impulse  decreases  by  as  much  as  10%  for  natu¬ 
ral  gas  (Groningen)  and  3%  for  methane,  but  increases 
up  to  14%  with  hydrogen  vitiation.  Similar  but  smaller 
variations  are  observed  for  c*  in  Fig  8  5.  Fig.  8.4  shows 
the  variation  of  c*  with  combustor  equivalence  ratio  for 
hydrogen  vitiation.  The  major  reason,  for  the  significant 
effect  of  hydrogen  vitiation  on  hydrocarbon  theoretical 
performance  is  the  decrease  in  molecular  weight  of  the 
vitiated  air  and  subsequently  in  the  ramjet  exhaust  (Fig. 
8.3).  These  results  show  that  vitiation  can  significantly  al¬ 
ter  ramjet  performance.  These  effects  must  be  accounted 
for  when  relating  the  data  generated  from  vitiated  air  fa¬ 
cilities  to  achievable  ramjet  performance  in  atmospheric 
air. 


8. 4. 2. 2  Example  Test  Results 

A  comparison  was  made  between  the  c*  efficiencies  ob¬ 
tained  from  data  of  some  previously  conducted  test  on  a 
ramcombustor  fed  by  vitiated  and  ideal  air.  The  ram¬ 
jet  fuel  was  kerosene.  Ideal  hot  air  was  provided  by  a 
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Figure  8.4:  Change  in  SFRJ  characteristic  velocity  when 
comparing  vitiated  air  for  a  hydrogen  fueled  vitiator  with 
ideal  hot  air[34j 
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Figure  8.5:  Effect  of  vitiator  fuel  on  SFRJ  characteristic 
velocity  when  comparing  vitiated  air  with  ideal  hot  air  [34] 
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Figure  8.6'.  Effect  of  vitiator  fuel  on  SFRJ  specific  impulse 
when  comparing  vitiated  air  with  ideal  hot  air  [34] 
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heat  exchanger.  In  case  of  vitiated  air  tests,  vitiator  fuel 
was  hydrogen  and  —  contrary  to  the  calculations  for  the 
above  figures  —  the  mass  percentage  of  oxygen  was  main¬ 
tained  at  23%  in  the  incoming  air  of  the  ramcombustor. 
The  highest  simulated  temperature  was  850A',  which  rep¬ 
resents  a  Mach  number  of  about  3.8  at  an  altitude  higher 
than  llirm. 

The  theoretical  performances  were  determined  using  one 
of  the  methods  of  Chapter  4  (method  3  of  Section  4. 4. 5. 3). 
The  computations  were  made  for  chemical  equilibrium  at 
station  2  and  4  with  frozen  flow  in  the  nozzle.  The  major 
results  are  shown  in  Figure  8.7. 

Tables  8.1  and  8.2  give  theoretical  and  experimental  re¬ 
sults  for  the  ideal  air  tests  and  the  vitiated  air  tests,  re¬ 
spectively.  Analysis  shows. 

*  an  increase  in  the  theoretical  c*  with  temperature 
similar  to  the  one  shown  in  Fig.  8.5  for  HTPB  in 
which  the  mole  fraction  of  O2  was  held  constant 

*  up  to  a  total  temperature  of  800/<’  at  station  2,  no 
significant  difference  between  the  c*  efficiency  with 
ideal  or  vitiated  air 

*  for  a  total  temperature  higher  than  800/^  at  station 
2  the  c'  efficiency  is  higher  with  ideal  air  than  with 
vitiated  air 

These  results  have  been  obtained  with  only  one  set  of  ex¬ 
perimental  data  and,  of  course,  it  is  necessary  to  make 
further,  tests  with  more  measurements  and  calculations. 
Nevertheless,  this  study  shows- that  caution  is  needed  when 
an  experimental  combustion  efficiency  is  determined  in 
connected-pipe  tests  with  vitiated  air. 


8.4.3  Fuel  and  Oxidizer  Choices 

The  choice  of  a  fuel  and  oxidizer  combination  for  heating 
air  in  a  particular  test  facility  is  complicated  by  all  the 
concerns  listed  in  Section  8.2.  Fig.  8.8  lists  some  common 
fuels  and  oxidizers  together  with  some  of  their  positive 
and  negative  attributes. 


8.5  Heater  Installation  and  Use 

The  performance  of  a  heater  is  significantly  affected  by 
the  way  it  is  installed  in  the  facility. 

The  flow  quality  of  air  to  the  ramjet  combustor  under  test 
can  be  modified  by  means  of  flow  straighteners,  screens 
or  plenums  to  reduce  turbulence  and  spatial  variations  in 
flow  pressures  and  temperatures  from  the  heater 


If  water  cooling  is  necessary  for  the  heater  and  the  ducting 
to  the  combustor,  heat  loss  may  induce  a  temperature 
gradient  in  the  flow.  Insulating  the  hardware  to  limit  heat 
loss  minimizes  this  problem. 

Acoustic  decoupling  of  the  air  heater  from  the  fuel  and 
oxidizer  feed  lines  and  the  ramjet  combustor  is  often  nec¬ 
essary  to  avoid  introducing  acoustic  oscillations  unique  to 
the  facility  installation  of  the  ramjet  under  test.  A  bas¬ 
ket  diffuser  having  many  small  holes  is  often  used  in  the 
air  heater  to  suppress  these  pressure  oscillations.  Sonic 
orifices  between  the  combustor  and  heater  prevent  feed 
back  from  the  ramburner  to  the  vitiator.  This  solution, 
however,  requires  higher  vitiator  supply  pressures  for  the 
air,  fuel  and  make-up  oxygen,  which  necessitates  higher 
pressure  supply  tanks  and/or  pumps. 

8.6  Heater  Performance  Determi¬ 
nation 

The  performance  of  a  heater  must  be  measured  during 
its  development  in  order  to  optimize  its  ability  to  prop¬ 
erly  condition  air  for  ramjet  testing.  The  methodology  for 
heater  performance  is  essentially  the  method  outlined  in 
Chapter  6  for  ramjet  combustors.  During  routine  testing, 
simple  monitoring  for  changes  in  operating  performance  is 
usually  sufficient  to  indicate  problems  that  may  affect  test 
results. 


8.6.1  Performance  Parameters 

Combustion  efficiency  is  the  most  important  heater  perfor¬ 
mance  parameter  since  unburned  heater  fuel  will  affect  the 
results  obtained  with  the  ramjet  tested.  The  temperature 
rise  efficiency  definition  as  given  in  Section  6.2.3  is  also 
appropriate  for  measuring  heater  performance,  especially 
as  long  as  the  exhaust  temperatures  are  below  the  lim¬ 
its  of  direct  temperature  measurement.  An  efficiency  that 
compares  actual  heater  fuel  flow  to  the  fuel  flow  theoreti¬ 
cally  needed  for  the  achieved  temperature  rise,  see  Section 
6.2.4,  also  clearly  characterizes  the  heater  performance. 

Uniform  temperature  and  pressure  profiles  in  the  vitiator 
exhaust  are  desirable. 

The  heater  should  have  a  smooth  combustion  behavior 
with  a  minimum  of  pressure  oscillations  and  be  free  of 
distinct  resonance  frequencies. 

Heat  loss  to  the  heater  structure  or  downstream  devices  is 
a  performance  related  parameter  since  it  increases  heater 
fuel  flow  and  subsequently  the  air  contaminants. 

The  envelope  of  heater  operation  should  be  wide  enough  to 
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Figure  8.7:  Change  in  LFRJ  characteristic  velocity  when  comparing  vitiated  air  with  ideal  air  [35] 


Test  case 

m 

ideal  air 
[tff/s] 

T,2 

exp 

[^] 

m 

C10H20 

c* 

th 

[m/s] 

c* 

exp 

[m/s' 

600K,  ssO.60 

1.623 

604 

0.066 

1119 

1084 

0.969 

'  600K,  4>  w  0.70 

1.602 

611 

0.076 

1172 

nil 

0,948 

700K,  <?S«0.60 

1.586 

686 

0.065 

1138 

1106 

0,972 

800K,  «S«0.50 

1.598 

815 

0.057 

1122 

1107 

0.987 

800K,  d«0.60 

1.632 

819 

0.067 

1164 

1137 

0.973 

850K,  <i) «  0.60 

1.544 

854 

0.063 

1169 

1163 

0.995 

850K,  4>  «  0.80 

1.563 

857  i  0.084 

1252 

1240 

0.990 

Table  8.1;  Ramcombustor  test  results  with  ideal  air 


Test  case 

1 

m 

ideal  air 

[icg/s] 

m 

O2 

\kg/s] 

m 

H2 

Tt2 

exp 

m 

Cio772o 

c* 

th 

[m/s] 

C* 

exp 

[m/s] 

'Jc* 

600K,  (^ssO.OO 

1,611 

52.3  X  10"^ 

4,95  X  10-"* 

605 

0,068 

1122 

1080 

0.962 

600K,  0  w  0.70 

1.618 

50,96  X  lO-^* 

4.83  X  10^^ 

601 

0.080 

1174 

1109 

0  945 

700K,  4>  w  0,60 

1.510 

67,20  X  10-^ 

6.32  X  10-^ 

703 

0.065 

1146 

1114 

0,972 

800K.  <;i«0.50 

1.433 

88.28  X  10-^ 

8,20  X  10"^ 

815 

0.056 

1137 

1106 

0,973 

'WOK,  0  s»0.60 

1.440 

87.23  X  lO''* 

8.10  X  10-=^ 

807 

0.063 

1168 

1129 

0,967 

850K,  <!>  w  0.60 

1,384 

94.56  X  10~3 

'"8,77  X  10-“ 

846 

0.061 

1177 

1146 

0.974 

,  850K,  w  0.80 

1.380 

94.44  X  10~^ 

8,76  X  10-^ 

849 

0.081 

1262 

1226 

0.972 

Thble  8.2;  Ramcombustor  test  results  with  vitiated  air 
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allow  for  the  temperature  and  mass  flow  variations  needed 
for  flight  trajectory  simulation. 


8.6.2  Experimental  Methods 

The  experimental  methods  applied  to- determine  the  per¬ 
formance  of  a  combusting  heater  correspond  to  standard 
rocket  or  ramjet  development  procedures. 

The  heater  can  be  equipped  with  a  choked  exhaust  noz¬ 
zle  to  provide  the  desired  Mach  number  in  its  combus¬ 
tion  chamber.  Subsequently  all  combustion  efficiencies  de¬ 
fined  in  Chapter  6  can  be  deduced  from  measured  mass 
flows  and  combustor  temperature  or  combustor  pressure 
(or  thrust,  but  this  is  rare), 

Sampling  and  analysis  of  the  heater  exhaust  products  can 
also  be  used  as  a  method  to  characterize  the  completeness 
of  heater  combustion  if  the  desired  heater  performance 
level  justifies  the  expense  of  this  more  sophisticated  per¬ 
formance  evaluation  technique.  In  addition  to  measuring 
concentrations  of  combustion  products,  this  method  can 
also  evaluate  levels  of  such  species  as  NOr  that  may  differ 
from  that  of  hot  atmospheric  air. 

The  uniformity  of  pressure  and  temperature  of  the  flow  at 
the  heater  exit  can  be  measured  directly  by  rakes  or  grids 
equipped  with  pitot  or  thermocouple  probes. 

The  steadiness  of  heater  combustion  can  be  checked  by 
unsteady  pressure  measurements. 

Specialized  non-intrusive  measuring  techniques  may  be  ap¬ 
plied  (at  considerable  expense)  if  comprehensive  charac¬ 
terization  of  the  heater  exhaust  is  desirable  or  if  high 
temperatures  prohibit  using  probes,  as  in  simulation  of 
hypersonic  flight  Mach  numbers.  Laser  velocimetry  can 
characterize  flow  velocity  profiles  in  lieu  of  measuring  pres¬ 
sure  profiles.  Nevertheless,  this  technique  requires  par¬ 
ticle  seeding  which  raises  the  problems  of  particle  injec¬ 
tion,  velocity  relaxation  and  particle  survival  depending 
on  the  flow  environment.  Spectroscopic  techniques  such 
as  Rayleigh  scattering,  Raman,  CARS,  LIF,  etc.,  can  be 
applied  to  characterize  temperature  profiles  and  combus¬ 
tion  products. 


8.6.3  Theoretical  Performance  Parame¬ 
ters 

Experimentally  derived  performance  parameters  need  to 
be  referenced  to  theoretical  values  calculated  for  complete 
or  equilibrium  combustion.  The  theoretical  values  most 
relevant  to  determining  a  combusting  heater  performance 
are  the  stagnation  temperature  and  the  mole  fraction  of 
species  at  the  heater  exit  (corresponding  to  a  combustor 


•  Hydrogen 

-t-  wide  flammability  and  ignition  ranges 
-I-  efficient  combustion  . 

-  water  vapor  and  low  molecular  weight  in 
exhaust 

•  Hydrocarbon  (liquid  or  gaseous) 

-b  minimal  safety  requirements 

-  carbon  dioxide  and  water  vapour  in  ex¬ 
haust 

•  MMHand.UDMH 

-  toxicity 

•  NHa 

-  water  vapor  in  exhaust 

-  toxicity 

•  NjO 

-I-  exhaust  is  similar  to  air 
+  heat  release  upon  decomposition 

-  special  handling  required 

•  N2O4 

-  toxicity 


Figure  8.8:  Fuels  and  oxidizers  for  combusting  heaters 

station  4  as  defined  in  Fig.  3.1).  Additional  parameters  as 
specified  in  Section  4.1  for  stations  4  and  5  are  required 
to  determine  the  stagnation  temperature  at  heater  exit 
indirectly  by  pressure  or  thrust  measurement  (Chapter  6). 

8.6.4  Performance  Monitoring 

Heater  performance  monitoring  at  the  test  facility  is 
mandatory  to  ensure  the  correctness  of  the  simulated  op¬ 
erating  conditions  for  the  ramjet  tested  and  for  test  safety. 
On-line  monitoring  during  the  test  is  desirable,  but  a 
heater  performance  check  by  post  test  data  reduction  may 
be  sufficient  if  the  overall  cost  of  the  test  and  test  article 
is  moderate, 

Heater  monitoring  should  look  at; 

•  the  injected  mass  flows  (heater  fuel,  make-up  oxygen, 
etc.) 

•  the  exit  temperature  (most  frequently  measured  by 
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thermocouples  installed  through  the  heater  wall, 
preferably  at  different  circumferential  positions). 

The  mass  flow  ratio  of  heater  fuel  to  make-up  oxygen  must 
be  kept  in  a  narrow  tolerance  band  to  maintain  the  correct 
oxygen  content  of  the  ait  fed  into  the  ramjet.  Combustion 
efficiency  of  the  heater  can  be  determined  on-line  from 
measured  data  with  modern  computers.  However,  on-line 
performance  monitoring  can  also  be  done  by  referencing 
the  measured  heater  mass  flow  to  a  min/max  tolerance 
table  as  a  function  of  temperature.  Depending  on  the 
possibility  of  tunning  heated  air  before  the  test  to  achieve 
steady  state  conditions  for  the  heat  losses  to  the  structure 
of  the  test  setup,  steady  state  or  transient  heat  transfer 
out  of  the  flow  must  be  determined  prior  to  testing  and 
be  taken  into  account  for  heater  performance  monitoring. 

The  uniformity  of  the  (circumferential)  heater  exit  tem¬ 
perature  must  be  monitored  within  a  predetermined  tol¬ 
erance  band.  The  width  of  the  band  should  take  into  ac¬ 
count  the  number  and  locations  of  the  thermocouples  and 
the  uncertainty  of  the  measurements. 

A  possible  flameout  of  the  heater  may  be  monitored  by 
using  thermocouples  to  measure  the  average  exhaust  tem¬ 
perature  or  a  higher  temperature  in  the  combustor  zone, 
or  by  using  optical  detectors. 
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9  Summary,  Conclusions  and  Recommendations 


In  the  initial  efforts  of  AGARD/PEP  Working  Group  22 
it  was  thought  that  perhaps  a  “recommended  method" 
could  be  established  within  NATO  for  the  determination 
of  connected-pipe  ramjet  and  ducted  rocket  performance. 
It  soon  became  apparent  that  there  wais  more  than  one 
acceptable  method.  The  result  was  a  detailed  description 
of  each  performance  parameter  and  the  interrelationships 
between  them,  and  an  expansion  of  the  Working  Group 
objectives  to  include  detailed  examples  in  order  to  provide 
users  with  a  working  document. 

Agreement  was  reached  (Chpt.  3)  on  the  methods  that 
should  be  used  for  reporting  test  results. 

Theoretical  performance  determination  (Chpt,  4)  depends 
upon  the  parameter  chosen  to  represent  performance  and 
the  specific  aerothermo chemical  equilibrium  code  and  op¬ 
tions  which  are  used.  It  wcis  found  that  the  various  codes 
(NASA  CET89,  PEP  . . .)  gave  substantially  the  same  re¬ 
sults,  providing  that  the  latest  properties  data  available 
were  used  as  input.  Specific  techniques  for  handling  viti¬ 
ated  air  heaters  have  been  recommended  and  illustrated. 

In  Chapter  5  the  experimental  parameters  that  are  re¬ 
quired  to  be  measured  are  given. 

Chapter  6  and  Appendix  C  present  the  fundamental  re¬ 
lationships  that  are  needed  to  calculate  experimental  per¬ 
formance  based  on  characteristic  velocity,  vacuum  specific 
impulse,  temperature  rise,  equivalence  ratio  and  pressure 
losses.  Various  “isentropic”  exponents  are  often  used  in 
performance  calculations.  An  explanation  of  these  are  pre¬ 
sented  in  Appendix  B,  together  with  recommendations  as 
to  which  values  are  most  appropriate  in  the  various  regions 
of  the  combustor  and  exhaust  nozzle. 

In  Chapter  7  detailed  sample  calculations  arc  presented, 
including  uncertainty  analysis.  This  chapter  presents  step- 
by-step  examples  for  the  use  of  the  methods  and  recom¬ 
mendations. 

Since  most  facilities  which  experimentally  evaluate  ram¬ 
jets  and  ducted  rockets  utilize  air  heaters,  some  detailed 
explanations  of  their  operation  and  recommended  tech¬ 
niques  for  their  use  have  been  included  in  Chapter  8,  es¬ 
pecially  for  vitiated  ait  heaters  where  a  variety  of  tech¬ 
niques  have  been  used  throughout  the  NATO  community. 
It  should  be  emphcisized  that  the  influence  of  air  vitiation 
on  ramjet  performance  is  dnly  generally  understood,  while 
there  is  still  a  lack  of  detailed  investigations  to  quantify 


the  resulting  effects. 

A  very  important  result  from  the  experimental  determi¬ 
nation  of  performance  is  the  uncertainty  of  the  final  re¬ 
ported  performance  parameter.  The  extensive  work  of 
AGARD/PEP  Working  Group  15  and  AGARD  Lecture 
Series  169  on  Comparative  Engine  Performance  Measure¬ 
ments  for  gas  turbines  was  utilized  in  the  present  effort 
(Appendix  A),  with  specific  examples  presented  for  ram¬ 
jets  and  ducted  rockets. 

Special  efforts  have  been  made  by  Working  Group  22  to 
emphasize  the  working  document  aspect,  in  order  to  facil¬ 
itate  the  practical  use  of  the  performance  determination 
procedures  given  in  this  report. 

In  conclusion,  it  is  the  hope  of  Working  Group  22  that 
this  document  will  prove  to  be  valuable  to  both  new  and 
experienced  investigators  in  the  area  of  ramjet  and  ducted 
rocket  performance  determination.  If  the  recommended 
procedures  are  followed,' the  document  should  permit  all 
investigators  to  be  able  to  communicate  their  results  in  a 
common  “performance  language”  and  to  readily  convert 
performance  values  obtained  in  one  facility  .to  those  used 
in. other  facilities. 

No  document 'is  without  some  shortcomings  and  very  few 
are  fortunate  to  be  without  a  few  errors.  It  is  hoped  that 
users  of  this  document  will  forward  their  comments  and 
recommendations  to  AGARD/PEP  Standing  Committee 
02  so  that  the  material  can  be  revised  and  improved  in  the 
future. 

A  recommended  next  step  is  to  perform  a  similar  effort  for 
supersonic  and  dual-mode  combustion  ramjets. 
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A  Uncertainty  Analysis  Methodology 


The  methodology  and  terminology  used  by  AGARD-PEP 
WG  15  in  their  interfacility  comparisons  of  turbine  engine 
performance  ([2])  has  been  adopted  in  this  report. 


A.l  Error  Types 


It  is  an  acceptable  fact  that  measurements  have  errors 
and,  therefore,  error  sources.  Elemental  error  sources  are 
classified  as  either  precision  {random  errors)  or  bias  (fixed 
errors).  The  following  sections  define  precision,  bias  and 
uncertainty  interval  (combined  error). 


A.  1.1  Precision  (or  Random  Error) 


Random  error  is  seen  in  repeated  measurements  of  a  single 
parameter.  Measurements  do  not  and  are  not  expected  to 
agree  exactly.  There  are  always  numerous  small  effects 
which  cause  disagreements.  The  variations  between  re¬ 
peated  measurements  can  be  quantified  by  the  precision 
index  (S). 


N  -  1 


(A.l) 


where  x  is  the  average  value  of  N  individual  measurements 
Xi,  in  the  sample. 


A. 1.2  Bias  (or  Fixed  Error) 

The  second  error  component  is  the  systematic  error,  which 
is  constant  for  repeated  measurements  and  can  only  be  de¬ 
termined  by  comparison  with  the  true  value  of  the  quan¬ 
tity  measured.  A  true  value  comparison  is  normally  im¬ 
possible  within  a  single  measurement  process,  but  tests 
can  be  arranged  to  provide  some  bias  information.  Exam¬ 
ples  are: 


1.  Interlab  and  interfacility  test  comparisons  on  mea¬ 
surement  devices,  test  rigs,  and  full  scale  engines. 

2.  Calibration  of  the  measuring  instruments  against  lab 
standards  during  the  test,  e.g.,  incorporating  a  stan¬ 
dard  in  the  scanning  cycle. 


3.  Comparisons  employing  redundant  instruments  or 
measuring  techniques. 


Large-  differences  in  measurements  can  usually  be  at¬ 
tributed  to  a  mistake,  but  this  progressively  gets  more 
difficult  as  the  size  of  the  difference  reduces.  Hence,  one 
tends  to  be  left  with  small  unexplained  differences,  which 
constitute  part  of  the  bias  limit. 


A. 1.3  Uncertainty  (Combined  Error) 

For  comparison  of  measurement  results,  a  single  value  is 
desirable  to  express  a  reasonable  error  limit  or  uncertainty 
interval.  This  value  must  be  a  relevant  combination  of 
bias  and  precision.  Precision  is  a  statistic,  which  lends 
itself  to  the  calculation  of  confidence  limits,  within  which 
the  actual  measurement  can  be  reasonably  expected  to  lie 
in  the  absence  of  bias  error.  It  is,  however,  impossible  to 
define  a  single  rigorous  statistic  for  the  total  error,  because 
bias  is  an  upper  limit,  which  has  unknown  characteristics, 
and  is  to  some  extent  dependent  on  engineering  judgment. 

Usually,  the  bias  (B)  plus  a  multiple  of  the  precision  index 
IS  used  to  estimate  the  total  error  or  uncertainty  interval 

(C/). 

t/  =  ±{B-t-t93S)  (A.2) 


in  which  tgg  is  the  95th  percentile  point  for  the  two-sided 
Student’s  ”t”  Distribution  ([3]). 


A.2  Error  Analysis  Process 


A  single  measuring  chain  stretches  from,  the  physical 
phenomenon  being  measured  (e.g.  pressure,  temperature, 
thrust),  via  probe  and  connecting  line,  to  the  transducer, 
and  from  there  usually  via  an  electric  line  -  sometimes 
preamplified  -  to  the  multiplexer  amplifier  and  signal 
conditioner  and  then  to  the  recorder.  Afterwards  the  sig¬ 
nal  is  played  back,  and  instrumental  calibration  applied, 
and  a  number  of  measurements  are  combined  to  determine 
a  value  representative  of  the  physical  phenomenon  be¬ 
ing  measured,  usually  by  averaging  in  space  and/or  time. 
Such  Basic  Measurements  are  then  used  to  calculate  Per¬ 
formance  Parameters,  (e  g.,  thrust,  combustion  efficiency. 
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Type 

Error 

Description  of  Error  Sources 

Standards 

Calibration 

bi 

■  ■Sl 

Standard  lab  calibration,  including  traceabil- 

Hierarchy 

ity  to  national  standards 

1 

Calibration  System  Errors 

•  Reference  Source  . 

•  Sensor 

h 

*2 

Determination  of  reference  pressure  . 

-  Hysteresis  . 

f>3 

S3 

Sensor  hysteresis 

Data 

-  Non-linearity  . 

S4 

Sensor  non-linearity 

Acquisition 

-  Repeatability  . 

frs 

S5 

Sensor  repeatability 

Recording  System  Errors 

•  Sampling . 

be 

S6 

Data  variations  due  to  facility  or  engine 
instabilities 

1 

•  Channel  . 

b? 

S7 

Signal  conditioning,  electrical  calibrations  and 
digital  systems 

Data 

Reduction 

Data  Processing  Errors  . 

be 

Sg 

Curve  fits  of  calibration  data 

Installation  Effects  Errors 

•  Pressure  Probe  . 

^9 

S9 

Design/fabrication  of  probes 

Other  Effects 

Environmental  Effects  Errors 

•  Temperature  .  . 

bio 

sio 

Temperature  change  effect  on  sensor 

•  Vibration  . 

bn 

S|1 

Vibration  effect  on  sensor  ' 

Root  Sum  Square 

Table  A.l:  Example  of  pressure  measurement  error  sources 


total  pressure  loss),  which  constitute  the  end  product  of 
the  measurement. 

Each  step  in  the  above-mentioned  measuring  chain  con¬ 
tributes  to  the  overall  data  error  in  its  own  specific  way 
and  is  treated  in  the  error  analysis  process  below, 

1.  Define  ' Elemental  Errors  (bias  and  precision)  for 
the  Basic  Measurements;  pressure,  temperature, 
force,  length  and  time. 

2.  Perform  Sensitivity  Analysis  to  determine  Influ¬ 
ence  Coefficients  and  the  combined  effect  (attendant 
bias  and  precision)  for  the  Performance  Parameters. 

3.  Estimate  Uncertainty  Interval  by  combining  total 
bias  and  precision  values  for  each  Performance  Para¬ 
meter. 


A. 2.1  Elemental  Error  Sources 

The  first  step  is  to  assess  and  categorize  the  elemental  er¬ 
rors  for  both  bias  and  precision,  in  a  separate  table  (e.g., 
Table  A.l  for  pressure  measurement)  for  a  single  point  of 
each  Basic  Measurement,  keeping  bias  limits  B  and  preci¬ 
sion  indices  S  strictly  apart.  Each  elemantal  error  source 


may  be  composed  of  bias  and/or  precision  error.  These  el¬ 
emental  errors  are  combined  by  Root-Sum-Square  (RSS) 
addition  to  give  the  total  B  and  S  values  for  each  Basic 
Measurement,  An  important  condition  required  to  justify 
RSS  combination  is  that  each  item  must  be  independent. 

The  Abernethy /Thompson  methodology  described  in  [3] 
details  the  evaluation  of  the  elemental  errors.  The  elemen¬ 
tal  error  of  a  single  measuring  chain  can  be  categorized 
into  four  groups  as  follows: 

1.  Calibration  Hierarchy 

2.  Data  Acquisition 

3.  Data  Reduction 

4.  Other  Effects,  e.g.,  non-instrument  effects,  errors  of 
method,  sensor  system  errors,  spatial  profile  sam¬ 
pling,  etc. 

For  the  purpose  of  conducting  a  detailed  assessment  of  the 
facility  measurement  uncertainties,  it  may  be  necessary 
to  define  error  subgroups  for  each  measurement  system 
(e.g.  Table  A.l  through  A. 4).  The  Uniform  Engine  Test 
Programme  elemental  error  groups  arc  documented  in  [2], 
The  general  definitions  of  the  elemental  error  groups  are 
given  below. 
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Type 

Error 

Description  of  Error  Sources 

Standards 

Calibration 

b: 

S] 

Manufacturer  specification  of  wire  or  standard 

Hierarchy 

lab  calibration 

Calibration  System  Errors 

•  Level  . 

h 

«2 

Reference  temperature  level 

Data 

•  Stability  . 

bs 

S3 

Reference  temperature  stability 

Acquisition 

Recording  System  Errors 

•  Sampling . 

b, 

S4 

Data  variations  due  to  facility  or  engine 
instabilities 

•  Channel  . 

b. 

S5 

Signal  conditioning,  electrical  calibrations  and 
digital  systems 

Data 

Reduction 

Data  Processing  Errors  . . 

be 

S6 

Curve  fits  of  calibration  data 

Installation  Effects  Errors 

•  Probe  Recovery  . 

b? 

S7 

Probe  design  caused  by  radiation,  convection, 

etc. 

Other  Effects 

•  Conduction  Error  . 

6a 

Ss 

Heat  conduction 

•  Temperature  Gradients 

69 

*9 

Temperature  gradients  along  nonhomoge- 
neous  thermocouple  wire 

Root  Sum  Square 

Table  A. 2:  Example  of  temperature  measurement  error  sources 


1.  Calibration  Hierarchy  traces  the  possible  instru¬ 
ment  error  back  to  the  National  Standard,  usually  in 
steps  via  a  Working  Standard,  a  Laboratory  Standard 
and  a  Transfer  Standard.  In  each  step  the  original 
bias  of  the  instrument  is  removed  by  the  calibration 
and  replaced  by  the  (smaller)  combination  of  system¬ 
atic  error  of  the  reference  instrument  and  the  random 
error  of  the  comparison.  (Additional  details  are  pro¬ 
vided  in  [2]). 

2.  Data  Acquisition  errors  can  be  caused  by  slight 
variations  in  exciter  voltage,  outside  influences  on 
data  transmission  and  on  the  transducer,  signal  con¬ 
ditioning  and  recording.  The  first  three  items  cause 
non-repeatability  (precision  error).  Another  factor  is 
sensor  hysteresis;  this  usually  depends  on  the  mea¬ 
suring  range  and  could  be  reduced  if  the  sensor  is 
only  calibrated  over  the  minimum  range  and  if  the 
measuring  history  is  known.  In  this  case,  hysteresis 
is  classified  as  bias.  Usually  this  is  not  a  practical 
proposition;  however,  with  modern  instruments  hys¬ 
teresis  is  small. 

3.  Data  Reduction  errors  consist  of  resolution  error 
and  calibration  curve  fit  errors  and  can  usually  be 
made  negligible,  compared  with  the  other  groups.  An 
error  of  half  the  biggest  error  elsewhere  only  con¬ 
tributes  10%  to  the  overall  error  when  added  RSS; 
therefore,  it  is  not  effective  to  use  extreme  resolution 


in  the  computational  hardware  and  software.  Cali¬ 
bration  curve  fit  errors  can  be  minimized  by  choosing 
the  appropriate  functional  relationship,  qualified  by 
visual  and  numerical  inspection. 

When  a  higher  than  second  order  curve  fit  is  used  it 
is  important  that  the  calibration  points  are  spaced 
evenly,  otherwise  the  densely  populated  part  may  in¬ 
troduce  a  calibration  bias  in  the  sparsely  populated 
part. 

4.’ Other  Effects  are  difficult  to  separate  and  as  such 
are  open  to  different  interpretations.  In  general 
they  are  concerned  with  the  interaction  between  the 
medium  and  the  measuring  chain.  This  is  the  case 
for  design  and  fabrication  of  probes  and  hole  patterns, 
which  renders  the  measured  pressure  sensitive  to  flow. 

Internal  flow  is  nearly  always  non-uniform,  both  in 
space  and  in  time,  and  not  necessarily  the  same  in 
different  installations.  This  nonuniformity  can  give  a 
bias  error  even  when  using  the  same  instrumentation, 
both  for  pressure  and  temperature.  Another  possi¬ 
ble  error  is  constituted  by  the  assumption  that  static 
pressure  is  constant  over  the  flow  area  of  the  parallel 
section  of  a  duct,  where  total  pressure  is  measured. 

The  mechanics  of  the  thrust  stand  can  introduce  bias 
and/or  precision  errors  —  notably  in  the  thrust  stand 
zero  —  which  can  not  be  determined  exactly,  not  even 


63 


Type 

Error 

Description  of  Error  Sources 

Standards 

r 

Calibration 

fei 

Si 

Standard  lab  calibration,  including  traceabil- 

Hierarchy 

1 

ity  to  national  standards 

Calibration  System  Errors 

0  Off-Axis  Effects  . 

62 

^2 

Force  measurement  on  axis  different  from 
centerline 

•  Tare  Correction 

—  Hysteresis  . 

S3 

Force  measurement' system  hysteresis 

Data 

-  Non-repeatability  .... 

Si 

Sensor  non-repeatability  from  repeat 
calibrations 

Acquisition 

—  Non-linearity  ........ 

bs 

S5 

System  calibration  non-linearity 

Recording  System  Errors 

•  Sampling . ; _ 

be 

S6 

Data  variations  due  to  facility  or  engine 
instabilities 

•  Channel  . 

b?' 

S7 

Signal  conditioning,  electrical  calibrations  and 
digital  systems 

Data 

Reduction 

Data  Processing  Errors  ... 

f>8 

ss 

Curve  fits  of  calibration  data 

Installation  Effects  Errors 

•  Stand  Alignment  . 

*9 

Misalignment  of  engine  and  data  load  cell 
force  vectors 

•  Load  Cell  . 

bio 

sio 

Shift  in  load  cell  calibration  caused  by 
attachments 

Pressure  Effects  Errors 

•  Load  Cell  . 

Cell  pressure  change  on  load  cell 

Other  Effects 

•  Test  Cell  . 

bi2 

S12 

Cell  pressure  change  on  test  cell  wall  ground 

•  Service  Lines  . 

bi3 

Sl3 

Line  pressure  change  on  tare  forces 

Temperature  Effects  Errors 

•  Load  Cell  . 

bi4 

Sl4 

Temperature  change  on  load  cell 

•  Test  Cell  . 

biS 

Sis 

Line  temperature  change  on  tare  forces 

•  Thrust  Stand . 

bie 

S16 

Thermal  growth  of  stand 

•  Vibration  Error  . 

bi7 

Sir 

Vibration  of  load  cell 

•  Scrub  Drag  Error . 

bi8 

•^18 

Secondary  airflow  external  drag 

Root  Sum  Square 

Table  A. 3:  Example  of  scale  force  measurement  error  sources 


in  an  end-to-end  calibration  (environmental  condi¬ 
tions  with  an  operating  ramjet  are  different  from  the 
calibration  environment).  Pre-test  and  post-test  ze¬ 
roes  are  different,  and  it  is  usually  assumed  —  but 
without  true  justification  —  that  the  test  zero  lies  in 
between. 

Length  and  time  can  generally  be  measured  very  ac¬ 
curately,  but  when  determining  flow  area  the  metal 
temperature  must  be  known  as  well  to  compensate 
for  growth.  Fuel  flow  depends  on  time  measurement, 
but  can  be  influenced  by  pulse  shape  and  by  residual 
swirl  with  less  than  10-20  diameters  of  pipe  straight 
section  upstream  and  downstream  of  the  flow  meter. 
In-situ  fuel  flow  calibration  is  preferred  but  discrepan¬ 


cies  still  exist.  Determination  of  fuel  properties  (lower 
heating  value  and  specific  gravity)  can  introduce  er¬ 
rors  of  0.3%  to  1%  because  of  reproducibility  and  re¬ 
peatability  of  evaluation  methods. 

Inadequate  sampling  or  averaging  is  another  error 
which  must  be  considered.  The  uncertainty  of  ef¬ 
fective  pressure  or  temperature  values  obviously  de¬ 
creases  with  the  number  of  probes  ([2]). 

When  comparing  the  value  of  a  Performance  Parame¬ 
ter,  which  is  a  dependent  variable,  it  is  necessary  to 
read  it  from  curves  at  a  chosen  value  of  an  indepen¬ 
dent  variable  (e.g.,  //a).  Any  uncertainty  in  the  cho¬ 
sen  independent  variable  translates  into  an  additional 
bias  error  in  the  performance  curve  (even  though  it 
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has  no  effect  on  the  individual  Performance  Parame¬ 
ter  values.)  A  procedure  for  combining  the  biases  to 
account  for  the  curve  shift  effect  is  described  in  [2], 


A. 2. 2  Sensitivity  Analysis 


The  second  step  is  to  propagate  B  and  S  separately  for 
each  Basic  Measurement  (or  Input  Parameter)  to  the  Per¬ 
formance  Parameter.  This  is  done  by  (1)  performing  a  sen¬ 
sitivity  analysis  to  determine  Influence  Coefficients  (7C) 
for  all  Input  Parameters  that  appear  in  the  Performance 
Parameter  calculation  and  (2)  multiplying  the  B  and  S 
values  by  the  appropriate  Influence  Coefficient  to  deter¬ 
mine  B  and  S  for  each  Performance  Parameter. 


S(Perf.  Param.)  = 


{A.3) 


S(Perf.  Param.)  = 


(A.4) 


The  standard  equations  used  to  calculate  Performance  Pa¬ 
rameters  from  the  Basic  Measurements  ate  presented  in 
chapter  6  of  this  report.  The  influence  of  an  error  in  any 
Basic  Measurement  on  the  outcome  can  be  determined 
either  by  Taylor  series  expansion  or  numerically  by  per¬ 
turbing  the  equation  for  a  difference  in  that  parameter, 
keeping  all  other  .parameters  constant  at  their  nominal 
values.  The  latter  method  is  preferred  when  used  with 
data  reduction'  software  because  it  accounts  for  implicit 
as  well  as  explicit  functional  relationships.  The  resulting 
Influence  Coefficient  is  usually  expressed  as  a  percentage 
variation  of  the  calculated  Performance  Parameter  (P)  for 
a  one  percent  deviation  of  a  single  Input  Parameter  (7). 


Influence  Coefficient  1C  = 


^P/P 

M/I 


(A.5) 


If  the  perturbation  is  small,  non-linearity  effects  will  be 
insignificant  —  but  of  course  the  value  of  the  Influence 
Coefficient  will  vary  over  the  operating  range  of  the  Per¬ 
formance  Parameter  and  is,  therefore,  a  function  of  ramjet 
operating  and  test  conditions. 


where  B  is  the  total  bias  error,  S  the  total  precision  er¬ 
ror  and  fgs  is  the  95th  percentile  point  for  a  two-sided 
Student’s  "t”  distribution  (a  function  of  the  number  of 
points  used  to  calculate  S)  for  the  respective  measure¬ 
ment/parameter.  If  the  predicted  5  is  determined  from 
a  large  number  of  points  {N  >  30)  the  value  <95  =  2.0 
can  be  taken;  Monte  Carlo  simulations  have  shown  that 
the  coverage  of  U  is  about  99  percent  ([36]).  This  means 
that  the  comparable  Performance  Parameter  results  from 
all  test  conditions  must  be  within  a  band  of  U.  If  this  is 
not  the  case,  either  a  data  error  exists  or  an  important 
aspect  of  the  uncertainty  estimate  has  been  overlooked. 

A  form  to  compute  estimated  uncertainty  for  a  single  Per¬ 
formance  Parameter  at  a  selected  test  condition  is  pre¬ 
sented  in  Table  A.5. 


A.3  Test  Data  Assessment 


When  the  pre-test  uncertainty  analysis  allows  corrective 
action  to  be  taken  prior  to  the  test  to  reduce  uncertainties 
which  appear  too  large,  the  post-test  assessment,  which  is 
based  on  the  actual  test  data,  is  required  to  refine  the  fi¬ 
nal  uncertainty  intervals.  Test  data  assessment  is  also  used 
to  confirm  the  pre-test  estimates  and/or  to  identify  data 
validity  problems.  It  can  also  be  made  to  check  for  consis¬ 
tency  if  redundant  instrumentation  or  calculation  methods 
have  been  used  in  the  data  collection  system. 


A.4  Glossary  for  Uncertainty  An¬ 
alysis  Methodology 


This  glossary  was  extracted  from  Ref  [2]. 


Accuracy  —  The  closeness  or  agreement  between  a 
measured  value  and  a  standard  or  true  value;  uncer¬ 
tainty  as  used  herein,  is  the  maximum  inaccuracy 
or  error  that  may  be  expected  (see  measurement 
error). 

Average  Value  (x)  —  The  arithmetic  mean  of  N 

readings.  The  average  value  is  calculated  as: 


Chapter  7  in  this  report  presents  results  of  a  sample  Sen-  ^ 

sitivity  Analysis.  *  =  average  value  =  ^ 

i  =  l 


A.  2. 3  Estimated  Uncertainty 

The  total  uncertainty  interval  for  both  Basic  Measure¬ 
ments  and  Performance  Parameters  is  estimated  by  (A. 2): 

U  =  i(77  -f  <9s5) 


Bias  (B)  —  The  difference  between  the  average  of  all 

possible  measured  values  and  the  true  value.  The 
systematic  error  or  fixed  error  which  characterizes 
every  member  of  a  set  of  measurements  (Fig.  A.l). 

Calibration  The  process  of  comparing  and  correct¬ 
ing  the  response  of  an  instrument  to  agree  with  a 
standard  instrument  over  the  measurement  range. 
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Type 

Error 

Description  of  Error  Sources 

r 

1  Standards 

Calibration 

T  : 

Si 

Standard  lab  calibration,  including  traceabil¬ 
ity  to  national  standards 

Hierarchy 

Calibration  System  Errors 

•  Cal  Fluid  Properties  .  ... 

62 

S2 

Mismatch  between  calibration  and  test  fluids 

•  Flowmeter  Repeatability  . 

^■3 

S3 

Non-repeatability  from  repeat  flowmeter 
calibrations 

Data 

Recording  System  Errors 

. 

Acquisition 

•  Sampling . 

64 

S4 

Data  variations  due  to  facility  or  engine 
instabilities 

•  Channel . 

65 

Ss 

Signal  conditioning,  electrical  calibrations  and 
digital  systems 

Data 

Reduction 

Data  Processing  Errors . 

k 

S6 

Curve  fits  of  calibration  data 

Installation  Effects  Errors 

•  Cavitation  . 

67 

St 

Insufficient  static  pressure  in  flowmeter 

•  Turbulence  . 

is 

Ss 

Sharp  bends,  etc.,  upstream  flowmeter 

•  Meter  Orientation . 

69 

Sg 

Orientation  difference  from  calibration  to  test 

Other  Effects 

Environmental  Effects  Errors 

e  Temperature 

—  Flowmeter  . 

iiQ 

sio 

Ambient  temperature  change  on  flowmeter 

—  Viscosity  . 

ill 

Sll 

Determination  of  test  fluid  viscosity 

—  Specific  Gravity . 

ii2 

S12 

Determination  of  test  fluid  specific  gravity 

•  Vibration  . 

ii3 

S13 

Vibration  on  flowmeter 

•  Pressure  .  . 

ii4 

Sl4 

Ambient  pressure  change  on  flowmeter 

Root  Sum  Square 

Table  A. 4:  Example  of  fuel  flow  measurement  error  sources 


Calibration  Hierarchy  —  The  chain  of  calibrations 
which  link  or  trace  a  measuring  instrument  to  a 
national  bureau  of  standards. 

Coverage  —  A  property  of  confidence  intervals  with 
the  connotation  of  including  or  containing  within 
the  interval  with  a  specified  relative  frequency. 
Ninety-five-percent  confidence  intervals  provide  95- 
percent  coverage  of  the  true  value.  That  is,  in  re¬ 
peated  sampling  when  a  95-percent  confidence  in¬ 
terval  is  constructed  for  each  sample,  over  the  long 
run  the  intervals  will  contain  the  true  value  95- 
percent  of  the  time. 

Cycle  —  A  whole  period  of  any  multiplexer. 

Data  Point  '  —  Can  be  made  up  from  a  number  of 
scans,  resulting  in  an  average  in  time  and/or  place 
(i.e.,  number  of  pick-ups). 

Defitred  Measurement  Process  (DMP)  —  encom¬ 
passes  the  overall  procedure,  including  calibration, 
etc.,  to  arrive  at  a  desired  test  result  using  a  spec¬ 
ified  installation  or  installations.  This  may  be  a 


single  test  point,  a  least  squares  curve  fit  to  a  num¬ 
ber  of  test  points,  or  a  collection  of  such  fits  for 
different  test  conditions.  Any  error  that  propagates 
to  the  result  as  a  fixed  error  is  classified  as  bias, 

.  otherwise  it  is  precision.  What  is  bias  for  a  single 
point  of  a  curve  becomes  precision  overall,  with  a 
remnant  test  bias  and  —  of  course  —  the  possibility 
of  an  installation  bias. 

Degree  of  Freedom  (df)  —  A  sample  of  N  values  is 
said  to  have  N  degrees  of  freedom,  and  a  statistic 
calculated  from  it  is  also  said  to  have  N  degrees  of 
freedom.  But  if  k  functions  of  the  sample  values  are 
held  constant,  the  number  of  degrees  of  freedom  is 
reduced  by  k.  For  example,  the  statistic 

»=i 

where  x  is  the  sample  mean,  is  said  to  have  N  — 
1  degrees  of  freedom.  The  justification  for  this  is 
that  (a)  the  sample  mean  is  regarded  as  fixed  or 
(b)  in  normal  variation  the  N  quantities  {x;  —  x) 
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are  distributed  independently  of  x  and  hence  may 
be  regarded  as  Af  —  1  independent  variates  or  N 
variates  connected  by  the  linear  relation  — 

x)  =  0. 

Dwell  —  Time  during  which  a  transducer  is  connected 
to  a  pick-up;  includes  settling  (line  or  filter  stabi¬ 
lization]  and  reading, 

Elemental  Error  — ■  The  bias  and/or  precision  error 

associated  with  a  single  component  or  process  in  a 
chain  of  components  or  processes, 

Fossilization  —  Random  (live)  errors  in  a  single  cali¬ 
bration  run  give  rise  to  an  uncertainty  in  the  value 
of  the  calibration  constants,  which  becomes  a  fixed 
“fossilized"  bizis  when  this  calibration  is  applied  to 
measurement  results. 

Laboratory  (Lab)  Standard  —  An  instrument 
which  is  calibrated  periodically  at  a  national  bu¬ 
reau  of  standards.-  The  laboratory  (lab)  standard 
may  also  be  called  an  interlab  standard. 

Mathematical  Model  —  A  mathematical  description 
of  a  system,  ft  may  be  a  formula,  a  computer  pro¬ 
gram,  or  a  statistical  model. 

Measurement  Error  —  The  collective  term  meaning 
the  difference  between  the  true  value  and  the  mea¬ 
sured  value.  Includes  both  bias  and  precision  er¬ 
ror;  see  accuracy  and  uncertainty.  Accuracy  implies 
small  measurement  error  and  small  uncertainty. 

Multiple  Measurement  —  More  than  a  single  con¬ 
current  measurement  of  the  same  parameter. 

Multiplexer  —  A  unit  which  connects  a  number  of 
pick-ups  sequentially  to  a  transducer,  or  a  number 
of  transducers  to  a  recorder. 

Parameter  —  An  unknown  quantity  which  may  vary 
over  a  certain  set- of  values.  In  statistics,  it  oc¬ 
curs  in  expressions  defining  frequency  distributions 
(population  parameters).  Examples:  the  mean  of  a 
normal  distribution,  the  expected  value  of  a  Poisson 
variable. 

Precision  Error  —  The  random  error  observed  in  a 
set  of  repeated  measurements.  This  error  is  the 
result  of  a  large  number  of  small  effects,  each  of 
which  is  negligible  alone. 

Precision  Index  (5)  —  The  precision  index  is  defined 

herein  as  the  computed  standard  deviation  of  the 
mezisurements. 

Usually,  _ 

but  sometimes 


Random  Error  Limit  of  Curve  Fit  (RELCF)  — 
The  limits  on  both  sides  of  a  fitted  curve  within 
which  the  true  curve  is  expected  to  lie,  with  95% 
probability;  apart  from  a  possible  bias  error  of  the 
DMP.  It  is  calculated  from  observed  random  statis- 
.tical  data,  including  the  Residual  Standard  Devia¬ 
tion. 

Reading  —  A  number  of  samples  or  an  averaged  value 
taken  during  a  dwell. 

Sample  —  A  single  value  giving  the  momentary  output 
of  a  transducer,  possibly  via  a  (low  pass)  filter. 

Sample  Size  (AT)  —  The  number  of  sampling  units 
which  are  to  be  included  in  the  sample. 

Scan  —  A  period  during  which  all  pick-ups  have  been 
read  at  least  once. 

Standard  Deviation  (cr)  —  The  most  widely  used 
measure  of  dispersion  of  a  frequency  distribution. 
It  is  the  precision  index  and  is  the  square  root  of 
the  variance,  5  is  an  estimate  for  tr  calculated  from 
a  sample  of  data. 

Standard  Error  of  Estimate  (See)  —  (also  known 
as  Residual  Standard  Deviation  (RSD))  The  mea¬ 
sure  of  dispersion  of  the  observed  dependent  vari¬ 
able  (Yobs)  about  the  calculated  least-squares  line 
(ycAi)  in  curve  fitting  or  regression  analysis.  It  is 
the  precision  index  of  the  output  for  any  fixed  level 
of  the  independent  variable  input.  The  formula  for 
calculating  this  is 

-  V - - 

for  a  curve  of  N  data  points  in  which  A'  constants 
are  estimated  for  the  curve. 

Standard  Error  of  the  Mean  —  An  estimate  of  the 
scatter  in  a  set  of  sample  means  based  on  a  given 
sample  of  size  N,  The  sample  standard  deviation 
(S)  is  estimated  as 

Then  the  standard  error  of  the  mean  is  S/'/N.  In 
the  limit,  as  N  becomes  large,  the  standard  error 
of  the  mean  converges  to  zero,  while  the  standard 
deviation  converges  to  a  fixed  non-zero  value. 

Statistic  —  A  parameter  value  baised  on  data,  x  and 
S  are  statistics.  The  bias  limit,  a  judgement,  is  not 
a  statistic. 

Statistical  Confidence  Interval  —  An  interval  esti¬ 
mate  of  a  population  parameter  based  on  data.  The 
confidence  level  establi.shes  the  coverage  of  the  in¬ 
terval.  That  is,  a  95-percent  confidence  interval 
would  cover  or  include  the  true  value  of  the  para¬ 
meter  95-percent  of  the  time  in  repeated  sampling. 
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Basic  Measurements 

Performance  Parameter 

Input 

Parameters 

/ 

Bias 

Limits 

Bi 

{%) 

Precision 

Index 

(%) 

Influence 

Coefficients 

1C 

(%/%) 

Bias 

Limits 

Bt  =  BiIC 
(%) 

Precision  | 
Index 

Sk  =  SilC 
{%) 

1 

2 

3 

n 

B=  % 

s=  % 

U=  %  1 

Table  A, 5:  Error  propagation  procedure  for  a  specific  performance  parameter  at  a  selected  test  condition 


Student’s  ”t”  Distribution  {<)  —  The  ratio  of  the 

difference  between  the  population  mean  and  the 
sample  mean  to  a  sample  standard  deviation  (multi¬ 
plied  by  a  constant)  in  samples  from  a  normal  pop¬ 
ulation.  It  is  used  to  set  confidence  limits  for  the 
population  mean. 

Traceability  — •  The  ability  to  trace  the  calibration  of 

a  measuring  device  through  a  chain  of  calibrations 
to  a  national  bureau  of  standards. 


Transducer  —  A  device  for  converting  mechanical  or 
other  stimulation  into  an  electrical  signal.  It  is  used 
to  measure  quantities  like  pressure,  temperature, 
and  force. 

Transfer  Standard  —  A  laboratory  instrument  which 
is  used  to  calibrate  working  standards  and  which  is 
periodically  calibrated  against  the  laboratory  stan¬ 
dard. 

True  Value  —  The  reference  value  defined  by  a  Na¬ 
tional  Bureau  of  Standards  which  is  assumed  to  be 
the  true  value  of  any  measured  quantity 

Uncertainty  (£/)  —  The  mtiximum  error  reasonably 

expected  for  the  defined  measurement  process. 

U  =  ±{B  +  tssS) 

Variance  (cr^)  —  A  measure  of  scatter  or  spread  of  a 

distribution.  It  is  estimated  by 

c2  _  ELifLill! 

*  "  V  -  1 


Measurement 


Figure  A.  1:  Sampling  systems 


from  a  sample  of  data.  The  variance  is  the  square 
of  the  standard  deviation.' 


Working  Standard  —  An  instrument  which  is  cali¬ 
brated  in  a  laboratory  against  an  interlab  or  trans¬ 
fer  standard  and  is  used  as  a  standard  in  calibrating 
measuring  instruments. 
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B  Isentropic  Exponents 


Simple  equations  relating  static  and  stagnation  properties 
can  be  derived  for  a  flowing  gas  if  the  following  assump¬ 
tions  are  made: 

•  nonreacting,  thermally  and  calorically  perfect  gas 

•  adiabatic  process 


combustor  and  exhaust  nozzle.  Two  limiting  flow  cases  are 
frozen  flow  and  local  equilibrium  flow.  Reaction  rates  are 
zero  for  frozen  flow,  whereas  local  equilibrium  flow  implies 
infinite  chemical  and  vibrational  rates.  The  specific  heats 
at  constant  pressure  and  constant  volume  are  defined  for 
local  equilibrium  and  frozen  flows  by  Eqs.  B.3  through  B.6 
[37). 


With  these  assumptions  the  equations  are  generally  in  the 
form  of  the  static  to  stagnation  ratios  for  a  property  (i.e. 
temperature,  pressure,  etc.)  as  a  function  of  Mach  number 
and  the  isentropic  exponent  7  (i.e.  Poisson  relations).  The 
process  is  easily  observed  with  the  help  of  a  temperature- 
entropy  diagram  (Fig.  B.l).  The  isentropic  exponent  de¬ 
rives  its  name  as  the  exponent  from  Eq.  B.l. 

=  constant  (B  1) 

7  =  —  (B.2) 


Specific  heats  at  constant  pressure: 


I 


Specific  heats  at  constant  volume; 

c„,,  =  c„.;-l-^ei 


This  value  of  7  does  not  change  for  a  fixed  property  flow 
and,  therefore,  it  is  constant  when  going  from  static  to 
stagnation  conditions. 


t 

where  -  indicates  a  molar  basis. 


(B.3) 

(B.4) 


(B.5) 

(B.6) 


The  difficulty  arises  when  trying  to  use  these  relationships 
for  chemically  reacting  flows  such  as  in  the  rocket  or  ramjet 


Figure  B.l:  General  temperature-entropy  diagram  for  a 
p  —  T  process 


One  can  see  that  for  local  equilibrium  flow  there  is  a  con¬ 
tribution  to  the  specific  heats  from  the  chemical  reaction. 
This  contribution  can  be  large  and  is  affected  by  a  change 
in  species  as  the  temperature  changes.  Therefore,  7  also 
changes  as  a  function  of  temperature  and  strictly  speaking, 
the  Poisson  relations  do  not  hold  any  longer  in  this  case. 
Nevertheless,  they  conveniently  approximate  the  process. 
Since  the  static  temperature  increases  when  going  from 
the  initial  static  condition  to  the  stagnation  condition  (re¬ 
fer  to  Fig.  B.2),  7  does  not  have  the  same  value  at  the 
stagnation  conditions  as  it  had  at  the  initial  static  condi¬ 
tions.  This  creates  the  problem  of  which  7  to  use  in  the 
isentropic  relations. 

The  isentropic  exponent  is  calculated  in  three  ways. 


Frozen  flow: 


7;  =  ^ 
c„(T) 


Local  equilibrium  flow: 


91np 

51nu 


7/ 

^_ln  V 
d  In  p  j. 


Process  7p: 


(B.7) 


(B.8) 
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isentropic  exponent  changes  for  a  chemically  reacting  equi¬ 
librium  flow  that  is  isentropically  expanded  from  a  stag¬ 
nation  condition  to  astatic  one.  Dissociation  or  recombi¬ 
nation  of  species  as  the  static  temperature  changes  causes 
a  change  in  thermochemical  properties,  such  as  specific 
heat  and  molecular  weight.  Even  for  frozen  flow  the  yy 
changes  as  the  static  temperature  changes.  A  process  isen¬ 
tropic  exponent  is  required  that  accounts  for  the  change  in 
composition  and  temperature  when  going  from  the  static 
condition  to  the  stagnation  condition. 

The  process  isentropic  exponent  (yp),  given  in  Eq,  B.IO, 
exactly  relates  the  properties  of  the  specified  end  states  of 
an  isentropic  process.  It  also  quite  accurately  relates  the 
properties  of  Intermediate  states  to  the  end  states.  For 
these  reasons  it  is  the  recommended  isentropic  exponent 
to  be  used  in  equations  of  Chapter  6. 


Figure  B.2:  Temperature-entropy  diagram  for  a  p-T  pro¬ 
cess 


In  general  the  process  isentropic  exponent,  yp,  can  be  de¬ 
fined  in  terms  of  pressures  and  temperatures  between  any 
two  stations  i  and  j  by 


.y  -.V  .  HPi/Pj) 

ln(p./p;)-ln(Ti/ry) 


(B.9) 


More  specifically  the  pressure  and  temperature  at  station 
i  can  be  either  static  or  stagnation  values,  leading  in  the 
second  case  to  the  following  form  of  yp  as 


^  _ _ HptifPj) 

ln(p„7P;)-ln(r,i/T;) 


(B.IO) 


For  the  process  between  stations  4  and  5  the  difference  in 
yp,4_5  obtained  using  the  above  two  definitions  has  been 
found  to  be  small  since  is  small. 


An  aerothermochemical  equilibrium  code,  such  as  the 
NASA  CET89  code  [7,  8,  9],  is  used  to  calculate  values 
for  yy  and  y,  as  defined  by  Eqs,  B.7  and  B.8,  respectively. 
Values  at  the  combustor,  the  nozzle  throat. and  the  nozzle 
exit  are  given  by  the  code  yp  as  defined  by  Eq.  B.IO  is  not 
calculated  directly  by  all  codes,  but  it  can  be  determined 
from  the  program  output. 

Both,  frozen  (yy)  and  local  equilibrium  (y, )  isentropic  ex¬ 
ponents  are  static  point  properties  of  the  flow  as  defined  by 
Eqs.  B.7  and  B.8,  respectively.  .The  combustor  yy  and  y, 
are  static  point  properties  when  a  combustor  static  pres¬ 
sure  is  entered  into  the  code  and  stagnation  point  prop¬ 
erties  when  a  combustor  stagnation  pressure  is  used.  A 
combustor  stagnation  pressure  must  be  entered  in  order 
for  the  code  to  correctly  calculate  the  static  point  prop¬ 
erty  Isentropic  exponents  for  the  nozzle  throat  and  nozzle 
exit. 


A  ramjet  engine  has  four  major  station  locations  where 
thermochemical  properties  are  needed  to  determine  theo¬ 
retical  and  experimental  performance.  They  are  the  inlet 
(station  2),  the  combustor  (station  4),  the  nozzle  throat 
(station  5)  and  the  .nozzle  exit  (station  6).  Characteristics 
of  the  flows  are  given  below. 

1.  Flow  in  the  inlet  is  in  chemical  equilibrium.  In  ad¬ 
dition,  a  yp  can  not  be  easily  determined.  For  low- 
values  of  A/2  (<  0.5)  y,  or  yy  can  .be  used  to  relate 
stagnation  to  static  properties. 

2.  Flow  in  the  combustor  is  chemically  reacting  and  is 
usually  in  local  equilibrium,  The  Mach  number  is  rel¬ 
atively  low  (<  0.5)  and  therefore,  the  static  properties 
are  not  very  different  from  the  stagnation  properties, 
y,  4  relates  these  states  with  good  accuracy,  but  it  is 
recommended  that  the  local  equilibrium  process  isen¬ 
tropic  exponent,  yp^,  =  yp,4_5,  be  used  for  compati¬ 
bility  with  nozzle  throat  calculations.  • 

3.  Flow  in  the  converging  nozzle  is  generally  in  local 
equilibrium.  Therefore,  it  is  recommended  that  the 
combustor-to-nozzle  throat  local  equilibrium  process 
isentropic  exponent,  yp,,  =  yp,4-5,  be  used. 

4,.  Flow  in  the  diverging  part  of  the  nozzle  generally  has 
fixed  composition.  Therefore,  it  is  recommended  that 
the  nozzle. throat-to-nozzle  exit  frozen  flow  process 
isentropic  exponent,  ypj  =  yp,5-6>  be  used. 


None  of  these  point  property  isentropic  exponents  are  the 
correct  ones  to  use  in  the  isentropic  relationships.  The 
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C  Compilation  of  Equations  for  Performance  Evaluation 


The  calculation  of  performance  is  based  on  a  simplified 
gasdynamical  model  of  the  engine  process; 

•  The  flow  through  the  engine  is  one-dimensional.  Over 
the  cross-section  of  the  stream  tube,  variable  values 
are  replaced  by  average  values. 

•  There  is  no  loss  of  air  or  fuel  within  the  engine  chan¬ 
nel, 

•  The  change  of  state' of  the  average  values  follows  the 
aerothermodynamical  rules  of  uniform  flow.  In  par¬ 
ticular,  the  average  values  of  the  stagnation  flow  prop¬ 
erties  can  be  derived  from  the  average  values  of  local 
static  flow  properties. 

•  The  geometrical  cross-sectional  areas  of  the  engine 
channel  (e.g.  the  nozzle  throat)  are  adjusted  by  flow- 
coefficients  in  order  to  obtain  conformity  between  the 
gasdynamical  functions  of  density  and  velocity  and 
the  law  of  mass  conservation, 

•  The  flow  conditions  in  the  convergent  part  of  the  noz¬ 
zle  are  particularly  important  for  the  performance 
evaluation  process.  In  this  nozzle  section,  isentropic 
flow  is  assumed  (no  change  of  total  temperature  or 
total  pressure).  Moreover,  for  the  sake  of  clarity,  the 
assumption  of  choked  nozzle  operation  is  made.  This 
condition  exists  in  the  majority  of  practical  tests. 

In  the  following,  the  equations  will  be  simplified  for  a  ther¬ 
mally  and  calorically  perfect  gas  (p  =  gRT,  Cp  ^  /(T)). 
In  practice  these  equations  are  also  applied  to  real  com¬ 
bustion  gases,  even  if  there,  are  comprehensive  deviations 
from  the  idealisation  of  thermal  and  calorical  perfection, 
by  modifying  the  coefficients  and  exponents  of  the  equa¬ 
tions. 

The  geometry  used  in  the  following  sections  is  shown  in 
Figure  C.l.  It  should  be  noted  that  is  commonly 

called  in  many  textbooks,  and  is  the  effective  flow  area. 
In  addition,  the  ^-location  is  often  used  for  station  5.  The 
7  used  in  this  Appendix  refers  to  the  process  7  (7,,)  in 
Chpt.  6  and  Appendix  B 


Figure  C.l:  Definitions  for  the  derivation  of  the  momen¬ 
tum  equation 


At  the  throat,  this  becomes 


Fss  =  rhjCs  4- P5/I5  (C-l) 

By  introducing  the  formulas  of  mass  flow 

rns  =  gsC^AsCoi  (C.2) 

and  dynamic  pressure  (with  M5  =  1,0) 

?5  =  ^Qscl  =  |p5M|  =  |p5  (C.3) 

the  stream  thrust  is  given  by  the  following  equation 

^55  =  esCsAscus -1- P5A5 

=  7P5A5CD5  +  P5A5 

Fsi  =  (C.4) 

For  some  deductions  it  is  more  convenient  to  avoid  the 
pressure  term  in  the  equation.  This  is  obtained  as  follows 


WsCs  +  P5A5 

'  P5A5 


rriscs  I  1  + 


(C.5) 


C.l  stream  Thrust 


By  combination  of  equations  C.2  and  C.3  a  useful  relation 
of  the  term  mjcs  is  derived. 


The  stream  thrust  is  defined  as 

Fs  =  rhc  -h  pA 


I 


mjcs  =  QiclA^CDi 

=  7PS.A  5  CD  5 


(C.6) 
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Using  equation  C.6  the  second  term  in  the  bracket  {C.5)  In  the  nomenclature  of  this  report 
reduces  to 

ViM  _  1  >  _  Pts^sCQs 

771505  7C£)5 

With  equation  C.7  a  different  formulation  of  stream  thrust 
is  obtained 


7715 


^  1+7CD5.,.  . 

fsi  =  - 7775^5 

7CD5 


Since  M  =  1  (nozzle  throat)  equation  C.8  specializes  to 

(C,9) 


Equation  C,16  offers  a  relation  between  and  c*,  which 
(C.8)  proceeds  to  a  useful  formula  by  application  of  equation 
C.ll. 


„  7CP5  +  1  .  * 

=  - TTJSO 


7CD5 


with  a*  being  the  critical  speed  of  sound. 

C.2  Local  and  Total  Parameters 


(r«m  +  l)(^) 


CDS 
7CD5  +  1 
CDS 
CDS 


7CD5  +  1 


(C.17) 

(C.18) 


By  Bernoulli’s  equation  for  compressible  flow  the  local  and 
total  parameters  are  related  in  the  following  way: 

^  =  (l  +  l^M=)^  (CIO) 


M  =  1 


Pts  _  _  / 7  +  1^ 

P’^  ~  Ps  ~  V  2  ) 


M  =  1  :  = 


T* 


Tg  _  7  +  1 
Ts 


—  -  -  i 

a  “  V  T  ~  V 
_  h  +  1 
a*  “  V  2 


C.3  Vacuum  Specific  Impulse  and 
Characteristic  Velocity 


C.4  Combustion  Temperature 


(C.ll) 

(C.12) 

(C.13) 


from  or  c* 


The  combination  of  equations  C. 9  and  C.15  gives  a  relation 
between  vacuum  specific  impulse  and  sonic  speed  at  total 
temperature  (combustion  temperature). 

•*  — 

*iiac  ■ 

7715 

7CD5  +  1  , 


1  +  (C.14) 


(C.15) 


7CD5 
7CDS  + 
7CD5 
7CD5  + 


7CDS 


a«5 


(C.19) 


With  help  of  the  equation  for  sonic  speed 
fliS  =  \/jRiTg 

a  relation  between  the  total  temperature  and  1*,,^  follows. 


is  defined  as  the  thrust  per  unit  mass  flow  of  a  conver¬ 
gent  nozzle  discharging  into  vacuum.  This  value  is  identi¬ 
cal  with  the  stream  thrust  per  unit  mass  flow  in  the  sonic 
throat. 

Having  in  mind  this  definition  and  with  the  help  of  equa¬ 
tion  C.4  the  following  relation  can  be  derived; 

^  _  ^55 

^vac  ■ 

7775 

=  + 

77I5 

=  +  ^  (C16) 

\PtsJ  7715 

The  definition  of  the  characteristic  velocity  is 

c*  =  £i^ 

m 


T  :=-{ ^  7+1 

2V7CD5+1/  7  Rs 


(C.20) 


The  combination  of  equations  C.20  and  C.18  leads  to  the 
relation  between  total  temperature  and  characteristic  ve¬ 
locity. 

.  Tls  =  7 

The  practical  evaluation  always  has  to  deal  with  the  prob¬ 
lem  of  choosing  the  proper  value  of  7.  It  is  recommended 
to  take  the  so  called  process-7,  which  norrnally  is  calcu¬ 
lated  by  the  aerothermochemical  equilibrium  code. 


The  more  exact  method  (which  accounts  for  variable  gas 
properties)  is  to  take  directly  from  the  aerothermochemi¬ 
cal  equilibrium  code  the  ratios  between  the  different  para¬ 
meters,  or  substitutions,  and  to  use  them  as  proportion¬ 
ality  factors.  For  example,  equations  C.20  and  C.21  turn 


over  to 


Because  of  the  relation 


or 


C.5 


(C.22) 

(C.23) 


Determination  of  the  Stream 
Thrust  by  Thrust  Measure¬ 
ment  with  a  Convergent  Noz¬ 
zle 


hts  _  m 
A5CD5  c* 

a  formula  for  total  pressure  is  formed. 


/  c*  \  ^5  +  PambAf, 


(C.28) 


By  introduction  of  equation  C.18  an  explicit  formulation 
can  be  geae;ated. 


P<5  = 


7  +  J5  +  VamhA^ 

1  )  (1  +  7C£,3)yl5 


(C.29) 


From  the  requirement  of  momentum  conservation  follows 
the  definition  of  ramjet  nozzle  thrust  (see  Fig.  C.l): 

P’s  =  FS5-P<.m6>l5  (C.24) 

This  can  also  be  written  as 

F's  =  -  (P6  -  Pamb)At  (C.25) 


Regarding  the  choice  of  7,  the  same  problem  exists  as 
mentioned  in  section  C.4,  Again  the  preferable  alterna¬ 
tive  is  to  use  equation  C.28,  taking  the  ratio  of  c*  and 
fJaj  as  the  theoretical  value  from  the  aerothermochemical 
equilibrium  code. 


Vtac/iA  A5CD5 


(C.30) 


where  \ 

Fmeai  =  measured  thrust  =  Flc  -  ^tare  (C.26) 


Flc  load  cell  force 
Fiart  thrust  stand  preload 

As  geometrical  cross-sectional  area  of  nozzle  throat 
pb  pressure  on  the  nozzle  base  area 


In  some  cases  the  effects  of  base  pressure  are  negligible. 
Then  the  difference  between  Fs,  and  Fmeat  is  very  small 
and  the  usual  form  of  the  momentum  relationship  as  seen 
in  textbooks  becomes 

F55  =  Fmeaz  +  PambA^  (C.27) 


C.6  Determination  of  Total  Pres¬ 
sure  in  the  Combustion 
Chamber  from  Thrust  Mea¬ 
surement  (Convergent  Noz¬ 
zle) 


From  the  definition  of  vacuum  specific  impulse  Eq.  6.9  and 
Eq.  C.24: 

Fs5  =  =  Fs  -b  Pam6  As 


Division  by  Ascds  gives 


ms  ^  : 

Asccs*'"'' 


Fs  -b  PflTTitAs 

AsCds 
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D 


Input  and  Output  Files  of  Sample 


Calculations 


The  following  pages  contain  the  input  files  used  with  the  NASA  CET89  aerothermochemical  equilibrium  code  ([9])  to 
perform  the  calculations  of  Chapter  7.  All  input  files  are  included  with  the  appropriate  sections  of  the  output  for  the 
Case  \  chamber  stagnation  conditions  and  for  the  Case  2  vitiated  air  inlet  conditions. 


For  better  understanding  the  spaces  or  blanks  in  input  files  are  written  as  readable  characters. 


fnputfile. 


REiCTASTS 

1  562uu0u  ■  4196^j(jAR  .  00934uCij  ■  000314i  m  it  U^n  nun  mlOO  .  OOt  u  I^I  n  n  u  u  e  iGl  II  II  p  t  H  lO 
IAHSL15TS 

u*IIPT2uRASE=l ,TP=T,Sr0irT=T,BSQl!=T,T=608.0,P=6S0200.0u/ 

Outputfile: 


TBERriQOYIlHIC  EQUILIBRIUM  PROPERTIES  AT  ASSICIEO 

( 

TEMPERATURE  AID  PRESSURE 


CASE  10.  1 

CHEMICAL  FORMULA 

FUEL  H  1.S6200  0  .41960  AR  .00934  C  .00031 


VT  FRACTIOI 
(SEE  lOTE) 
1.000000 


EIERGY  STATE  TEMP 
KJ/KC-MCL  OEO  R 

.  000  G  .  00 


0/F*  .0000  PERCEIT  FUEL*  100.0000 


EQUIVALEICE  RATIO*  .OOlS  PHI*  .0000 


TBERMOOYIAMIC  PROPERTIES 


P,  MPA 

. 6S020 

T,  OEG  K 

606.00 

RBO,  KG/CU  H 

3.7377  0 

H,  KJ/KG 

310 . 87 

U,  KJ/KG 

136.91 

G,  KJ/KG 

-3962 . 01 

S,  KJ/{KG>(K) 

7.0S10 

M,  MOL  WT 

28.96S 

(I)LV/DLP)T 

-1 . 00000 

<DLV/OLT)P 

1.0000 

CP.  KJ/(KG)(K) 

1.0525 

GAMMA  <SJ 

1 . 3750 

SOB  VEL,K/SEC 

489.1 

HOLE  FRACTI055 

AR 

.00934 

CD2 

.00031 

S2 

.78089 

02 

.  20946 

Figure  D.l:  Input  and  output  file  1  (Case  1)  for  NASA  CET89;  station  2  conditions 
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Inputfile: 

REACTAITS 

Hul .  562y[jDy ,  4196jju&R  .  00934l^u  4  000314i  h  n « n  ■  n  w  1 » ilOO.Oft  n « i21S2  4  ^  jGi  14  n  1 44  n  *  h  iD 
CylO  4 i_nxju®u2® ■  I  n « 11  ■  H II  ■  n  H II  w « II  ■  i>  *1 "  "  M  *1  ■ » II  ■  II  “  *100 4 OOtj"67G07  4 ft  »  h  ■  11 « >  h  iF 


B4NELI5TS 

u*HPT2uKA5E=2,RXT=T,5HraiI=T,15Ql!=I,0F«T,«IX«2l4S17685, 
uuP=S 68800. Ou/ 

utRKTISPuFR0Z=F,SU81R=l 47989Su/ 

Outputfile: 


TBEDRETICIL  rocket  PERFORMiBCE  ASSUKIBG  EqUILI8RIUH  CDKPDSITIDI  DinilBG  EXPABSIOB 
FROH  IIFIIITE  AREA  COMOUSTOR 


PIIF  »  8245  P5IA 

CASE  10 4  2 

CSEItlCAL  FORltULA 

OIIDAIT  1  1.56200  0  .41960  IR  .00934  C  .00031 

FUEL  C  10.00000  B  20.00000 


UT  FRACTIOB 
(SEE  BOTE) 

1 4  000000 

t. 000000 


EBERGY  STATE 
KJ/KG-MOL 
900s. 642  G 
-282867.700  L 


TEMP 
DEG  K 
.00 
.00 


0/F=  21.5177  PERCEIT  FUEL=  4.4410  EQUIVALEBCE  RATIO*  .6877  FBI*  .6872 


CBAHBER 

T8R0AT 

EXIT 

PIBF/P 

1 4  0000 

1.8082 

1 .0793 

P,  HPA 

.56880 

.314S6 

. 52701 

T,  DEG  K 

206S .80 

1834 . 46 

2035.14 

RBO,  XG/CU  H 

9.5724-1 

5.9639-1 

9.0035-1 

B,  KJ/KG 

207. SI 

-124. 4S 

162.50 

U,  KJ/KG 

-386.70 

-6S1.89 

-422.84 

G,  KJ/KG 

-17901.3 

-16205.3 

-17677.5 

5,  KJ/(KO)(K) 

8 . 7660 

8.7660 

8.7660 

M,  MOL  VT 

28.906 

28.918 

28.908 

(DLV/DLP)T 

-1.00022 

-1.00006 

-1 .00019 

(DLV/DLT)P 

1 .0073 

1.0021 

1.0063 

CP,  KJ/(KG)(K) 

1.4845 

1.4041 

1.4715 

GAMMA  (5) 

1.2444 

1.2588 

1.2465 

SOB  VEL,M/5EC 

859.9 

814.8 

854.2 

MACS  BUH8ER 

.000 

1.000 

.351 

PERFORHABCE  PARAMETERS 

AE/AT 

1.0000 

1 . 7989 

CSTAR,  M/SEC 

.1170 

1170 

CF 

.696 

.2S6 

IVAC.  M/SEC 

1462.1 

2251.0 

ISP,  M/5EC 

614  .8 

300.0 

MOLE  FRACTI0B5 

AR 

.00891 

.00891 

.00891 

CO 

.0003S 

.00006 

. 00029 

CD2 

.09147 

.09180 

.09154 

B 

.00001 

.00000 

.00001 

B2 

.00007 

.  00002 

.00006 

820 

.09079 

.09131 

.09089 

BO 

.00474 

.00245 

.00438 

102 

.00001 

.00001 

.00001 

52 

.74232 

.74379 

.74257 

0 

.00011 

.00002 

. 00009 

08 

.00129 

.00046 

.00114 

02 

.05993 

.06118 

.06012 

Figure  D.2:  Input  and  output  file  2  (Case  1)  for  NASA  CET89;  first  run  with  p4 
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Inputfild; 

REiCTiSTS 

■  S62i  II  lOi^i  1 41 9Ci  i!  lAR .  00934i  ^  i  .00031  4l  miimm«iibwii  ilOO  OOyuij21S2  1 iGi  » ii  ■  n  » « n  ;0 
C(jlO .  1 M  ■  II  iHj  i20 .  lll■llJl■l■#^l■lilllllll^^^llMll<llllll^l1llli>»i  ilOO .  00u“67 607  .  Q  iL  im  n  n  «  ■ « iF 


HINELISTS 

utIIPT2uRlSE=2,MT=T.SIUSIT=T.ISB«=T,0F=T,MH=21  .S1768S, 

uuP =6 13838.0 

uIERD 

utRRTISPuFR0Z=F,SU81R=l  7989S 
utEIO 

Outputfile: 


T8E0aETICAL  ROCKET  PERFORJIAECE  A55UHIIG  EQUILI8RIUH  CQHPaSITIOI  OURIIQ  EXPAISIOE 
FROH  ISFISITE  AREA  COHSUSTOR 


PIEF  =  89.0  PSIA 

CASE  Ea.  2 

C8EKICAL  FORMULA 

OXIDAET  E  1.S6200  0  .41960  AR  .00934  C 

FUEL  C  10.00000  8  20,00000 


WT  FRACTIOl 

EBERGY 

STATE 

(SEE  EOTE) 

XJ/KG-MOL 

,00031 

1.000000 

900S.642 

G 

1 .000000 

-282867 . 700 

L 

0/F=21.S177  PERCEET  FUEL=  4.4410  EQUIVALEECE  RATIO*  .6877  P8I*  .6872 


C8AHBER 

T8R0AT 

EXIT 

PIEF/P 

1 .0000 

1.8083 

1 . 0793 

P.  MPA 

.61384 

. 33946 

. 56873 

T,  OEG  K 

2066 . 02 

1834.49 

2035.32 

R80 ,  XG/CU  K 

1.0329  0 

6.43S9-1 

9.715S-1 

8,  XJ/XG 

207  I  51 

-124.49 

162.49 

U,  XJ/XG 

-386. 7S 

-6S1 , 93 

-422  89 

G,  XJ/XG 

-17858.0 

-16165.4 

-17634.5 

S,  XJ/(XG)(X) 

8.7441 

8.7441 

8.7441 

H,  MOL  WT 

28 . 906 

28.918 

28.909 

(OLV/OLP)T 

-1.00022 

-1.00005 

-1.00018 

(OLV/DLT)P 

1.0071 

1 , 0020 

1.0061 

CP,  KJ/(KG)(X) 

1 .4825 

1.4035 

1 . 4697 

GAMMA  (5) 

1 . 2447 

1 . 2S89 

1.2467 

50E  VEL,H/SEC 

860.0 

814.9 

854.3 

MACH  BUHBER 

.000 

1.000 

.351 

PERFORHAECE  PARAMETERS 

AE/AT 

1 .0000 

1.7989 

CSTAR,  M/SEC 

1170 

1170 

CF 

.698 

.256 

I VAC,  H/SEC 

1462.1 

22S1.0 

ISP,  H/SEC 

814.9 

300.1 

HOLE  FRACTIOES 

AR 

, 00891 

. 00B91 

.00891 

CO 

. 00034 

. OOOOG 

.00027 

C02 

. 09148 

.09180 

. 091 SS 

8 

.00001 

. 00000 

. 00001 

82 

.00007 

.00001 

.00006 

820 

. 09081 

.09132 

.09090 

10 

. 00474 

.0024S 

.00438 

B02 

. 00001 

.00001 

. 00001 

E2 

.74233 

.74379 

.74258 

□ 

.00011 

.00002 

.00009 

□8 

,00127 

.00045 

.00112 

02 

.05993 

.06118 

.06012 

TEMP 
DEG  X 
.00 
,00 


Figure  D.3;  Input  and  output  file  3  (Case  1)  for  NASA  CET89;  for  measured  using  7 
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Inputfile; 

RE&CT&ITS 

B|jl .  562ijlj0u  ■  4196uijAR ,  00934uCu ,  000314i  n  d  ]i » n » u  ■  n  ilOO ,  OOljiju21  52  ■  4i  iGi » w  ■ « n  o « lO 
CljI  0  ■  LIULiU^U^O  ■  I  II  a  II  II  n  n  ■  ii  ii  «  ii  ii  n  ii  n  ■  ii  ii  n  ii  i  ii  ii  n  ii  it  n  ii  (i  i  ti  ilOQ  ■  0Olj“67607  ■  Oi  ila  ii  n  a  n  n  a  n  iF 

B4HELISTS 

Li»iePT2uKiSE=2  ,RKT=T  ,SHniIT=T  ,«SQ!!=T  ,OF=»T  ,HII=21 , 517685 . 

uuP=61 3906.0 

utEIO 

LjtRETIHPuFR0Z=F.SU81R=l  ,79895 
uREID 

Outputfile: 


TBECRETICAL  ROCKET  PERF0RH4BCE  tS5(miSG  EQUILIBRIUH  CDHPOSITIQI  DURIIG  ElPllSIOB 
FROM  IIFIIITE  IREl  C0HBU5T0R 


PIIF  =  89.0  PSIl 

CASE  IG.  2 

C8EHIC1L  FORKULA 

OXIDIBT  I  1.56200  0  .41960  IR  .00934  C  .00031 

FUEL  C  10.00000  8,20.00000 


VT  FRACTIOI 
(SEE  BOTE) 

1 . 000000 
1 . 000000 


EBERGY  5TATE 
KJ/KG-MDL 
9005.642  G 
-282867 . 700  L 


TEKP 
DEG  X 
.00 
.00 


0/F-  21.5177  PERCEBT  FUEL=  4.4410  EQUIVALEBCE  RATIO*  .6877  P8I=»  .6872 


CBAMBER 

T8R0AT 

ElIT 

PIBF/P 

1 . 0000 

1.6063 

1 . 0793 

P,  KPi 

.61391 

. 33949 

. 56879 

T,  DEG  X 

2066 . 02 

1834  49 

2035 . 32 

RIO,  KG/CU  K 

1.0331  0 

6.4366-1 

9.7166-1 

B,  KJ/KG 

207 . 51 

-124  49 

162.49 

U,  KJ/KG 

-386.75 

-651.93 

-422  89 

G,  KJ/KG 

-17857'.  9 

-16165.4 

-17634,5 

S,  KJ/(XG)(K) 

8.7441 

8.7441 

8.7441 

B,  HOL  RT 

28.906 

26  918 

28.909 

(DLV/DLP)T 

-1 .00022 

-1 , 00005 

-1 , 00018 

(DLV/DLT)P 

1.0071 

1.0020 

1 . 0061 

CP,  KJ/(Ka)(K) 

1.4825 

1 . 4035 

1 . 4697 

GAMMA  (S) 

1.2447 

1.2569 

1.2487 

SOB  VEL, M/SEC 

860.0 

814.9 

654.3 

MACS  BUMBER 

.000 

1.000 

.351 

PERFORHABCE  PARAMETERS 

AE/AT 

1 , 0000 

1 , 7989 

CSTAR,  M/SEC 

1170 

1170 

CF 

696 

.256 

IVAC,  M/SEC 

1462.1 

2251.0 

ISP.  M/SEC 

.814,9 

300.1 

HOLE  FRACTIOBS 

AR 

. 00691 

. 00891 

.00891 

CO 

. 00034 

. 00006 

.00027 

ca2 

.09148 

.09180 

.09155 

6 

.00001 

. 00000 

. 00001 

62 

.00007 

.00001 

. 00006 

820 

.09061 

.09132 

. 09090 

10 

.00474 

.00245 

.00438 

102 

. 00001 

.00001 

. 00001 

B2 

.74233 

.74379 

,74258 

0 

.00011 

.00002 

. 00009 

08 

.00127 

.00045 

.00112 

02 

.05993 

.06116 

.06012 

Figure  D,4;  Input  and  output  file  4  (Case  1)  for  NASA  CET89;  for  measured  p4  without  using  7 
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Inputfile; 

REICTIBTS 

BljI  ■  562(jlj®li  ■  4196ijij^®-  'OO^S^yCij  ■  0003 14i  n  n » it  i  n « n »  jIOO ■  OOijlj(j2152  ■  4uGljuujuulajO 

CulO  •  ULIUuHu20  ■  I II  B  n  ]>  fi  II  B  II  LI  u  II  II  1  »  II  n  II  n  B  II  a  ii  ii  p  ii  ii  »  ii  n  n  ii  1 100  ,  OOu^  6760?  .  OutuULUULf  JuF 

liHEUSTS 

j4ISPT2u*ASE=2,MT=T.SIU1IIT=T,HSQM=T,0F=T,MI1=21, 517685, 

uuP=619848.0 

u»EBO 

uIKETIBPuFK0Z=F,SUHAR=1 .79895 
uAEND 

Outputfile: 


THEOHETICAL  ROCKET  PEKFORMABCE  ASSUHIBG  EqUILIHRIUH  COHPOSITIOI  DUAIBG  EIPABSIOl 
FKDM  IBFIBITE  AREA  COMBUSTOR 


PIBF  =  89.9  PSIA 

CASE  BO .  2 

CHEMICAL  FORMULA 

OIIDAHT  B  1.56200  0  41960  AX  .00934  C  ,00031 

FUEL  C  10.00000  H  20.00000 


WT  FXACTIOI  EBERGY  STATE  TEMP 

(SEE  HDTE)  .  KJ/KG-MOL  DEG  X 

1.000000  9005.642  G  .00 

1,000000  -282867,700  L  ,00 


D/F=  21.5177  PEXCEHT  F0EL=  4.4410 


EqurVALEHCE  RATIOS  .8877.  PHI=  .6872 


CEAHEEX 

TBXDAT 

EXIT 

PIFF/P 

1.0000 

1 , 8083 

1 . 0793 

P,  MPA 

.61985 

. 34278 

.57430 

T,  DEG  X 

2066.05 

1834.50 

2035 . 35 

RHD,  XG/CU  M 

1.0430  0 

6.4988-1 

9. 8106-1 

8,  KJ/KG 

207 . 51 

-124.49 

162.49 

U,  XJ/XG 

-386,76 

-651 . 94 

-422,89 

G,  XJ/XG 

-178S2.4 

-16160.4 

-17629.0 

S,  XJ/<XG)(X) 

8.7413 

8.7413 

8 .7413 

M.  MOL  UT 

28 . 906 

28.918 

28.909 

(DLV/DLPlT 

-1 .00021 

-1.00005 

-1.00018 

(DLV/DLT)P 

1 . 0071 

1.0020 

1,0061 

CP,  XJ/(XG)(X) 

1.4823 

1 . 4034 

1 . 4695 

GAMMA  (S) 

1 . 2447 

1.2589 

1.2468 

SOB  VEL, M/SEC 

860.1 

814.9 

854.3 

MACH  HU118EX 

.000 

1.000 

.351 

PERFQRHAHCE  PARAHETEXS 


AE/AT 

1 . 0000 

1 . 7989 

CSTAR,  H/SEC 

1170 

1170 

CF 

.696 

.256 

IVAC,  H/SEC 

1462.1 

2251.0 

ISP,  M/SEC 

814.9 

300,1 

MOLE  FXACTIOBS 


AR 

.00891 

.  00891 

.00891 

CO 

. 00034 

. 00006 

.00027 

CD2 

.09148 

.09180 

.09155 

H’ 

. 00001 

.00000 

.00001 

H2 

.00007 

.  00001 

.00006 

E2a 

.09081 

.09132 

. 09090 

HO 

.00474 

.00245 

.00438 

H02 

.00001 

.00001 

.00001 

B2 

.74234 

. 74379 

.74258 

0 

.00011 

. 00002 

.00009 

□H 

.00126 

, 00045 

.00112 

02 

.05993 

.06118 

06012 

Figure  D.5'.  Input  and  output  file  5  (Case  1)  for  NASA  CET89;  Pt4  using  measured  thrust  and  7 
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Inputfile: 

REACTiSTS 

•  562uuQu  ■  4196uuA8 , 00934|jCu .  000314i » n »  h  ■ «  n  i  n  ilOO .  00yij^21S2  ■  4i  iGi  n  n  n  n « ■  n  lO 
C[jlO  1 1 II  n  n  iHi  i20 1 1 II I  ti  II  ^  n  B II II H II II II II II  ■  ji  II H II  a  ii  ii  ii  ii  ii  ii  ii  ii  ■  ii  ilOO .  0^j“67 SOT  1  Oi  iL  11 II  It  II II II II  iF 


BAKELISTS 

u*IHPT2ultASE=2,RI[T=T.51UBIT-T,HSqH»T.0F"T.NII*21.5176H5, 

uuP= 620477,0 

uAEEO 

u»RKTISPuFR0Z=FiSU8AR=1 .79H95 
ukEBO 

Outputfile; 


T8E0RETICAL  ROCKET  PERFORXAECE  A5SUH1IG  EQUILISRIUH  COHPOSITIOB  OURIBS  EIPABSXDE 


FRQN  lEFlBITE  AREA  C0HBU5T0R 


PIBF  *  90,0  PSIi 

CASE  EO.  2 

C8EKICAI,  FORJniLA 

OIIOABT  B  1.S6200  0  .41960  AR  .00934  C  .00031 

FUEL  C  10.00000  8  20.00000 


VT  FRACTIOB  EIERGY  STATE 
(SEE  BOTE)  KJ/KG-KOL 
1.000000  900S.642  G 

1 . 000000  - 2B  2B67 , 700  L 


TEKP 
OEG  K 
,00 
.00 


0/F=  21.S177  PERCEBT  FUEL=  4,4410  EqUIVALEBCE  RATIO*  .6B77  PBI=*  .6872 


CBAKBER 

TBROAT 

EXIT 

PIBF/P 

1.0000 

1  .B083 

1 . 0793 

P,  KPA 

.62048 

.34313 

. S7488 

T,  DEG  K 

2066. OS 

1B34.S0 

203S  .3S 

RED  I  RG/CU  H 

1,0441  0 

6.S0S4-1 

9,B20S-1 

E,  KJ/KG 

207 , 51 

-124.49 

162.49 

U.  KJ/KG 

-3B6.76 

-6S1.94 

-422.90 

0,  KJ/KG 

-178S1,8 

-16159. B 

-1762B. S 

S,  KJ/(KG>(K> 

B.7410 

8.7410 

B.7410 

K,  KOL  VT 

2B  906 

28.918 

28.909 

(DLV/DLP)T 

-1.00021 

-1 .OOOOS 

-I.OOOIB 

(OLV/OLT)P 

1 . 0071 

1.0020 

1 . 0061 

CP,  KJ/(KG)(K) 

1 .4B22 

1  4034 

1 . 469S 

GAHHA  (S) 

1 , 2447 

1  2SB9 

1 . 246B 

SOB  VEL,H/5EC 

B60.1 

B14.9 

B54  3 

KAC8  IUHEER 

.000 

1.000 

.351 

PERFOHIUBCE  PARAHETERS 

AE/AT 

1 . 0000 

1 .7989 

CSTAR,  M/SEC 

1170 

1170 

CF 

.696 

.2S6 

I VAC,  H/SEC 

1462.1 

22S1.0 

ISPi  K/SEC 

814.9 

300, 1 

HOLE  FRACTIOIS 

AR 

. 00E91 

.00891 

.00891 

CO 

, .  00034 

, 00006 

.00027 

C02 

.09148 

091BO 

.09155 

E 

, 00001 

. 00000 

.00001 

E2 

,  00007 

.00001 

.00006 

E20 

.  090E1 

, 09132 

.09090 

10 

.00474 

. 0024S 

. 0043B 

102 

,00001 

,  00001 

,00001 

B2 

,74234 

,74379 

,74258 

0 

.00011 

.00002 

, 00009 

08 

.00126 

.00045 

.00113 

02 

.0S993 

,06118 

, 06012 

Figure  D,6.'  Input  and  output  file  6  (Case  1)  for  NASA  CET89;  pti  using  measured  thrust  without  7 
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Inputfile: 

RE1CTAIT5 

Bijl  •  SGSuijOiji  4136i_i;jAR  ■  00934|jC(j  ■  0003 14i»L>tiit  tlOO .  00|_fytj2152 . 4^_jG^j^jijuuLgj^jO 

CulO  .  UULAJ®u20  ■  UUUUULIUUUUJULJJUJJUUULl-'ULajULUUULIulOO.  00li“67607  .  OuLuuubjijlXjuF 

H1HELIST5 

o»ISPT2uRASE=2,RRT=T,5IUSIT=T,e55H>T,ERATI0=T,MII-HII=0. 55, 0.555, 0.56, 0.565. 0.67, 0.575, 0.50, 0,585, 0.59,0. 595, 0,60. 

uuP=613906.0 

uAEIO 

utRKTIIPuFR0Z>F,SU8AR=l .79895 
uAEIS 

Outputfile; 

T8E0XETICAL  ROCKET  PERFORHAICE  ASSUMIIG  EQUILIBRIUM  COHPOSITIOI  DURIIG  EIPAISIOI 
FROH  IIFIIITE  AREA  COHBUSTOR 

PIHF  =  09.0  P5IA 

CASE  SO.  2 

WT  FRACTIOI  EIERGY  STATE  TEMP 
C8EHICAL  FORMULA  (SEE  HDTE)  KJ/KG-HOL  OEG  K 

OH  DA  ST  H  1.56200  0  .41960  AR  .00934  C  .00031  1.000000  9005  642  G  OO 

FUEL  C  lO  .OOOOO  8  20.00000  1.000000.-282067.700  L  .00 

0/F=  26,9200  PERCEHT  FUEL=  3.5817  EQUIVALENCE  RATIO*  .5500  PBI=  .5493 


C8AM8ER 

T8RDAT 

EXIT 

PISF/P 

1 , 0000 

1.0183 

1,0801 

P,  MPA 

.61391 

.33763 

.56836 

T,  OEG  K 

1825.52 

1607 . 73 

1796,30 

R80,  KG/CU  M 

1.1700  0 

7.3072-1 

1.1009  0 

8,  KJ/KG 

227,51 

-67 , 174 

187 . 39 

U,  KJ/KG 

-297.18 

-529  23 

-328.09 

G,  KJ/KG 

-15356,7 

-13792.2 

-15147  5 

S,  KJ/(KG)(K) 

8  5369 

8.5369 

8.5369 

H,  MOL  UT 

28.928 

28  930 

28.928 

COLV/DLP)T 

-1.00004 

-1.00001 

-1.00003 

(OLV/DLTJP 

1  0014 

1,0004 

1 .0012 

CP,  KJ/(KG)(K) 

1.3785 

1.3313 

1.3716 

GAMMA  (5> 

1 . 2643 

1.2755 

1 . 2659 

SOS  VEL, M/SEC 

814.5 

767  ,7 

808.4 

HAC8  BUM8ER 

,000 

1.000 

.350 

PERFORMABCE  PARAHETER5 


AE/AT 

1  0000 

1,7909 

C5TAR,  M/5 EC 

1094 

1094 

CF 

,702 

.  259 

I VAC,  M/5EC 

1369.6 

2105.9 

ISP,  M/5EC 

767,7 

283,3 

HOLE  FRACTIOHS 


AR 

, 00899 

. 00699 

,00899 

CO 

. 00003 

. 00000 

.00002 

CD2 

, 07414 

. 07417 

,07415 

B2 

. 00001 

. 00000 

.00001 

820 

. 07360 

.07302 

, 07371 

SO 

.00200 

.00120 

.00262 

B02 

. 00001 

.00001 

,00001 

B2 

,75052 

.75138 

,75066 

D 

. 00002 

. 00000 

. OOOOl 

08 

00036 

.00010 

.00031 

02 

. 08936 

. 09024 

.00951 

Figure  D.7;  Input  and  output  file  7  (Case  1)  for  NASA  CET89,  equivalence  ratio  versus  characteristic  velocity 
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Inputfile: 


XEtCTiBTS 

B|jl  •  562u;j0|j  ■  4196|jyAR .  00934 yCy .  000314i  h  t  w  ■  ■  ■  m  h  ]^j6 . 383y^jiji » t  n  i  n  lO 


Cul  -UUiXluBu^  .1  a  ■  1  HMJI  ai  n  m  a  >  m  ■  m  lu  nintJuimmaiiiiiH  X)  .047i  mi  ■  ii  ii  ■  ii  »L*-pQLF_»-fLJUL*-‘U^ 
0y2  .  tjuLftJUUJUUJUUUUl  a  H  I  ILK  B.B  ^iJUUUUUUlJUUUJULijUL^UuQ  •  262^LJjyijjljUjG(juuuyLKJuO 


BiKELISTS 

u*IBPT2yUSE=3.TP=T.Srum=T,SSqH*T.T*=«06.0pP=650200.0 

uREBD 

Outputfile: 

TaEKHOOrimC  EgUILIBXIUH  PXOPEXTIES  IT  iSSIQIED 


TEHPEXATURE  AID  PKESSUXE 


CASE  10.  3 


¥T  FRACTIOI 

EIERGY 

STATE 

TEMP 

CHEHICAL  FORMJIA 

(SEE  lOTE) 

KJ/KG-HOL 

DEG  X 

OIIOAIT 

I  1.56200  0  .41960  AR  .00934  0 

.00031 

. 960572 

.000 

G 

.00 

FUEL 

C  1.00000  H  4.00000 

1.000000 

.000 

G 

.00 

OIIOAIT 

0  2.00000 

.039428 

.000 

G 

.00 

0/F=141.3830  PERCEIT  niEL=  .7023  EQUIVALEICE  RATIO"  .1089  PHI"  .1078 
THERHODriAXIC  PROPERTIES 


P.  MPA 

.66020 

T,  DEG  K 

606 .00 

RHO,  KG/CU  H 

3.7305  0 

H,  KJ/KG 

-69.041 

U,  KI/KG 

-243.34 

G,  KJ/KG 

-4389.49 

S,  KJ/(KG)(KI 

7.1294 

M,  MOL  HT 

28 . 908 

(DLV/DLP)T 

-1.00000 

(DLV/DLT)P 

1.0000 

CP,  KJ/(KG){K) 

1.0677 

GAMMA  (S) 

1 .3687 

SOI  VEL, M/SEC 

488.4 

MOLE  FRACTIOIS 

AR 

.00889 

CO  2 

.01295 

H20 

.02531 

12 

. 74339 

02 

.  20946 

Figure  D.8:  Input  and  output  file  1  (Case  2)  for  NASA  CET89;  station  2  conditions 
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Inputfile: 

REACTilTS 

I,  il  ■  486780, , .  47013,  .AR  .00889,  C,, .  01 295,  ,8. .  .05062.  ,.,6.692  .  „  .-476 . 7,  ,0,  „ . . .  „  „ ,  ,Q 
CulO  .  [  II  H  It  H[  |20 . 1 II  n  mi  1 11  n  mi  t  II  M  ■  ii  ii  i  ii  ii  w  ]i  i  ]i  ii  n  u  ii  ■  ii  ii  »  ii  n  lO  .  31  lij“67607  .  Oi  iL  »  ii  i  ii  ii  n  ii  iF 

lAMELISTS 

uAIIPT2uKA5E=4  ,  RICT=T ,  SIUFIT=T  ,SSI}II=T ,  P=S68800 . 0 
uAEED 

uARKTIIPuFR0Z=F,SU8AR=l .79895 
uAEBD 

Outputfile: 

T8E0RETICAL  ROCKET  PERFORHAICE  ASSUDIBG  EqUILI8RIUn  COHPOSITIOS  DURIRG  EIPABSIOB 


FROM  IRFIRITE  AREA  C0H8U5T0R 


PIRF  >  83.5  PSIA 

CASE  HO.  4 

C8EHICAL  FORMULA 

OIIDAHT  H  1.48678  0  .47013  AR  .00889  C  .01295  8  .05062 

FUEL  C  10.00000  8  20.00000 


9T  FRACTIOI  EHERGY 
(SEE  lOTEl  KJ/KG-MOL 
1.000000  -1994.513 

1.000000.  -282867.700 


STATE 

G 

L 


TEMP 
OEG  X 
.00 
.00 


0/F=  21.5177  PERCEIT  FUEL=  4.4410  EQUIVALEHCE  RATIO=  .7201  P8I=  .6859 


CIAH8ER 

T8R0AT 

El  IT 

PIHF/P 

1  0000 

1 . 8055 

1 .0791 

P,  MPA 

. 56880 

.31504 

.52711 

T,  OEG  R 

2044.75 

1819.23 

2014.87 

R80,  KG/CU  H 

9  6531-1 

6.0118-1 

9 . 0789-1 

8,  KJ/KG 

-155.49 

-484.14 

-200.02 

U,  KJ/KG 

-744 . 73 

-1008.17 

-780.60 

G,  KJ/KG 

-18196.8 

-16535.7 

-17977.8 

S,  KJ/(KG)(K) 

8.8233 

8.8233 

8.8233 

H,  MOL  UT 

28 . 853 

28.865 

28.855 

(DLV/DLP)T 

-1.00022 

-1 . 00006 

-1 .00018 

CDLV/DLT)P 

1 . 0073 

1.0021 

1.0063 

CP,  KJ/(KG)(K) 

1  5071 

1 .4283 

1 . 4940 

GAMMA  (5) 

1 . 2404 

1.2543 

1 . 2424 

SOJ  VEL, M/SEC 

854-9 

810.7 

849.3 

HAC8  HUM8ER 

.000 

1.000 

.351 

PERFORMAHCE  PARAMETERS 


AE/AT 

1 .0000 

1.7989 

C5TAR,  M/5EC 

1167 

1167 

CF 

.  695 

.256 

IVACi  M/SEC 

1457. 1 

2243 .9 

I5P,  H/5EC 

810.7 

.  298.4 

HOLE  FRACTI0H5 


AH 

.00848 

.00848 

.00848 

CO 

.  00033 

1 00006 

.00027 

CQ2 

.10336 

.10368 

.10344 

8 

. 00001 

.00000 

.00001 

82 

I  00008 

.00002 

.00007 

820 

.11474 

.11528 

.11484 

la 

00439 

.00228 

. 00406 

102 

1 00001 

.00001 

. 00001 

82 

.  70680 

.70816 

.70703 

0 

.00010 

. 00002 

.00008 

08 

.00132 

.00047 

.00117 

02 

.06037 

.06164 

.06055 

Figure  D.9:  Input  and  output  file  2  (Case  2)  for  NASA  CET89;  first  run  with  pn 
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Inputfile: 

REiCTilTS 

Hyl  •  4fl67 SOy  1470 13y &R  •  00  SdSyCy  •  01 235(j6y  •  0  5D62yyy6  <  692y(jij“476  ■  7 u^uuliulilliu^ 

CylO  .  I  ]|  n  It  iHi  i20 . 1  rl  It  H  It  It  It  c  it  ti  n  ii  ii  n  it  ii  i  it  ii  ii  ii  a  it  it  n  it  il  n  ii  it  ■  II  tt  |0 . 3l  lu“67007  ■  O  jLi  ti  it  it  it  it  r  n  iF 


SAHELISTS 

ytrHPT2yKiSE=4  ,RKT=T  ,SIU]IIT=T  ,1ISQJ!=T  ,P=613792 , 0 
yREID 

ukRKTIHPuFR0Z=F,SU84R=l .79895 
uRESD 

Outputfile: 


T8E0RETIC4L  ROCKET  PERFORJliSCE  ASSOKIBG  EOUILTBRIUK  COHPOSITIOB  DURIBG  EIPiJSIOB 
FROH  IIFIIITE  AREA  COMBUSTOR 


PIJF  =  89.0  PSIA 

CASE  10.  4 

C8EHICAL  FORMULA 

OIIDABT  B  1.48678.  0  .47013  AR  .00889  C  .01295  8  .05062 

FUEL  C  10.00000  8  20.00000 


WT  FRACTIOB 

EBERGY 

STATE 

(SEE  BOTE) 

EJ/KG-HOL 

1.000000 

-1994.513 

C 

1 . 000000 

-282867 . 700 

L 

TE)0> 
DEG  K 
.00 
.00 


0/F=  21.5177  PERCEBT  FUEL=  4.4410  EQUIVALEBCE  RATIO"  .7201  P8I=  .6859 


C8AH8ER 

T8R0AT 

EIIT 

PIBF/P 

1.0000 

1.8056 

1 . 0791 

P,  HPA 

.61379 

. 3399S 

. 56680 

T,  DEG  K 

2044 . 96 

1819.27 

2015.05 

RBO,  KG/CU  H 

1 . 0416  0 

6.4870-1 

9.7962-1 

8.  KJ/KG 

-155.49 

-484.17 

-200.02 

U,  KJ/KG 

-744 .78 

-1008 . 21 

-780.65 

G,  KJ/KG 

-18153.9 

-16496 . 2 

-17935.1 

S,  KJ/(KG)<K) 

8.8013 

8.8013 

8.8D13 

H,  HOL  UT 

28.8S3 

28.865 

28 . 855 

(DLV/DLP)T 

-1 .00021 

-1.00005 

-1 .00018 

(DLV/DLT)P 

1 .0071 

1 .0021 

1.0061 

CP.  KJ/(KG)(K) 

1 . 5051 

1.4257 

1.4923 

GAHHA  (5) 

1 . 2406 

1  2544 

1.2426 

SOB  VEL.H/SEC 

855.0 

810.8 

849.4 

HAC8  BUH8ER 

.000 

1.000 

.351 

PERFORHABCE  PARAMETERS 

AE/AT 

1 .0000 

1 .7989 

CSTAR,  H/SEC 

1167 

1167 

CF 

.695 

.  256 

IVAC,  H/SEC 

1457.1 

2243.9 

ISP,  H/SEC 

810.8 

298.4 

HOLE  FRACTIOBS 


AR 

. 00848 

. 00648 

. 00848 

CO 

.00032 

. 00006 

.00026 

C02 

.10338 

.10369 

. 10345 

E 

.00001 

. 00000 

.00001 

62 

.00008 

,00002 

. 00006 

620 

.11476 

.11529 

.11486 

BO 

.00439 

.00228 

.00406 

G02 

.00001 

.00001 

.00001 

S2 

.70681 

.70816 

,70703 

a 

.00009 

.00002 

.00008 

06 

.00129 

.00046 

.00115 

02 

.06037 

.06154 

.06055 

Figure  D.IO:  Input  and  output  file  3  (Case  2)  for  NASA  CET89;  pti  for  measured  without  using  7 
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Index 


Area 

error  measurement,  13,  40,  42 
geometric,  11,  22,  26-31,  43,  45,  46,  70-72 
propellant  burning,  23,  46 

Burning  area,  23,  46 
Burning  rate,  6,  12,  23,  46 
Burning  time,  12,  23,  46-48 
Burnt  thickness,  23 

Characteristic  velocity,  12,  26-27,  33-35,  44-45,  71-72 
Combustion  efficiency  based  on 

characteristic  .velocity,  26,  34,  36,  44,  52 
equivalence  ratio,  27,  35,  36,  45 
temperature  rise,  27,  35,  36,  45,  53 
vacuum  specific  impulse,  26,  35,  36,  44 
Composition 

of  air,  11,  12,  16-19,  30,  32,  43-44 
of  fuel,  2,  11,  12,  16,  30,  43 
of  vitiated  air,  11,  12,  17,  43-44 

Efficiency 

combustion,  11,  12,  25-27,  34-36,  44-45,  51,  53, 
55,  56 

error  measurement,  37 
expulsion,  12,  25,  28 
nozzle  expansion,  12,  25,  29 
Enthalpy 

of  air,  17,  20,  32,  43,  51 
of  fuel /propellant,  11,  16-20,  31,  43 
of  vitiated  air,  17 
Equivalence  ratio,  12 
burned,  26-28,  45 
error  propagation,  42 
injected,  27,  28,  35 
stoichiometric,  11,  21,  27,  34,  44 
Error  measurement,  12,  60-67 

Force 

error  measurement,  12,  39,  63 
load  cell,  13,  22,  26,  30,  31,  43 

Gas  constant,  27,  32,  35,  36,  44,  45 

Heat  of  formation,  xvi,  11,  18,  43 

Influence  coefficient,  12-13,  36,  37 
Isentropic  exponent,  12,  14-15,  25-29,  32-36,  43-46, 
68-72 

error  propagation,  42 


error  measurement,  13,  41 
of  air,  12,  16,  27,  30,  41 
offuel,30,  41 

of  fuel/propellant,  11,  12,  16,  22,  23,  27,  46 
of  vitiated  air,  16,  22 

Molecular  weight,  14,  23,  32-35,  43-45,  50,  51 

Nozzle  discharge  coefficient,  11,  23,  26,  30,  37,  44,  70- 
72 

Pressure,  12-14,  16,  22,  25-28,  30-36,  71-72 
error  measurement,  12,  37,  41,  61 
gas  generator,  23  ,  46-47 
losses,  25,  28,  36,  37,  45 

Specific  impulse,  12,  14,  21,  52 
Station  (vehicle),  10-11 

Temperature,  12-14,  16-20,  22,  25-27,  30-36,  43-45, 
49,  56,  71-72 

error  measurement,  12,  38,  41,  62 
Thrust 

coefficient,  15 
net,  7 

nozzle,  8,  9,  32,  34,  43 
stream,  xvi,  31,  37,  43,  70-72 

Vacuum  specific  impulse,  26-27,  33-35,  44-45,  71-72 


Mach  number,  6-9,  12,  20,  25,  26,  29,  33,  36,  45,  49 
Mass  flow  rate 
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